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Abstract 



The description of baryons as chiral solitons of the Nambu-Jona-Lasinio (NJL) model is 
reviewed. A motivation for the soliton description of baryons is provided from large Nc QCD. 
Rigorous results on the spontaneous breaking of chiral symmetry in QCD are discussed. It 
is then argued that the NJL model provides a fair description of low-energy hadron physics. 
The NJL model is therefore employed to mimic the low-energy chiral flavor dynamics of 
QCD. The model is bosonized by functional integral techniques and the physical content of 
the emerging effective meson theory is discussed. In particular, its relation to the Skyrme 
model is established. 

The static soliton solutions of the bosonized NJL model are found, their properties dis- 
cussed, and the influence of various meson fields studied. These considerations provide strong 
support of Witten's conjecture that baryons can be understood as soliton solutions of effective 
meson theories. The chiral soliton of the NJL model is then quantized in a semiclassical fash- 
ion and various static properties of the nucleon are studied. The dominating 1/ 'Nc corrections 
to the semiclassically quantized soliton are investigated. Time-dependent meson fluctuations 
off the chiral soliton are explored and employed to estimate the quantum corrections to the 
soliton mass. Finally, hyperons are described as chiral solitons of the NJL model. This is done 
in both, the collective rotational approach of Yabu and Ando as well as in the bound state 
approach of Callan and Klebanov. 
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1 Introduction 



It is generally accepted that Quantum Chromo Dynamics (QCD) is the theory of strong 
interactions (for an introduction see e.g. refs. [1, 2]). Besides coming in three colors the 
fermion fields of QCD, the quarks, also carry flavor. The interactions of QCD are flavor blind 
but sensitive to color. QCD is an asymptotically free theory which means that the forces 
between quarks become weak for small quark separations, or equivalently, large momentum 
transfers. This allows one to quantitatively calculate observables of strong interaction physics, 
which are sensitive to the short distance behavior of QCD, by perturbative techniques. As a 
matter of fact, the predicted scaling violations have been verified to a high accuracy at existing 
accelerators [1]. According to present knowledge QCD is the only renormalizable theory that 
can account for these scaling violations. 

The same self-interactions of gluons which give rise to asymptotic freedom lead to a strong 
quark-quark interaction for medium and small energies. Using additionally the empirical fact 
that neither quarks nor gluons have been detected as "free" particles has lead to the confine- 
ment hypothesis: Only singlets of the gauge group appear as physical particles. Especially, the 
quarks belonging to the fundamental representation and the gluons being in the adjoint repre- 
sentation can never be observed directly. The perturbative result" that the coupling constant 
increases as the momentum transfer becomes small, or the distance large, is in accordance 
with the confinement hypothesis. Unfortunately, this behavior of the coupling constant ex- 
cludes perturbative calculations for low energies. Thus it is still unproven that QCD is really 
a confining theory. Even worse, until today properties of hadrons have not been calculated 
from QCD without making use of severe assumptions or simplifications. 

There is probably one exception to this statement: With the help of Monte-Carlo tech- 
niques attempts have been made to calculate hadron properties directly from QCD using a 
discrete lattice (for an introduction to lattice gauge theories see e.g. refs. [3] or [4]). Despite 
the fact that these non-perturbative calculations should enable the determination of every 
physical quantity the results often deviate very strongly from the experimental values. There 
are several reasons for this. Even by using "high-performance" computers the possible lattice 
sizes are still modest, especially, if one wants to include dynamical fermions. Furthermore, a 
calculation using massless quarks is impossible. And there are still open conceptual questions 
concerning the continuum limit of lattice theory. 

Given this state of affairs it is natural to resort to effective models of strong interactions. 
These are intended to mimic the low-energy behavior of QCD as closely as possible. For 
this purpose the approximate chiral symmetry of QCD [5] provides a very useful guideline. 
Hadron phenomenology has left no doubt that this symmetry is broken dynamically by strong 
interactions. Requiring the known pattern of explicit, spontaneous and anomalous chiral 
symmetry breaking puts significant constraints on possible models for the strong interactions 
of quarks. Additional guidance can be obtained if one generalizes QCD to a gauge theory 
with an arbitrary number of colors Nc- This is because for large Nc QCD reduces to an 
effective theory of infinitely many weakly interacting mesons and glueballs [6] . Unfortunately, 
this effective meson theory cannot be constructed explicitly. Nevertheless, Witten was able to 
give arguments that within this effective theory baryons emerge as soliton solutions [7]. 

Although Witten's conjecture has never been proven rigorously the soliton picture of 
baryons has turned out quite successful in recent years. The starting points have been phe- 
nomenological effective meson theories, which possess soliton solutions. The most popular 

a See e.g. Chapter 3.4.2 in ref.[2]. 
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ones perhaps are the Skyrme model [8, 9] and the gauged cr-model [10, 11]. Investigations 
within these models have satisfactorily explained the wealth of spectroscopic baryon data, see 
e.g. ref. [12] for a recent compilation of references on soliton models for baryons. It is worth- 
while to mention that many of the early difficulties encountered for the Skyrme model, like e.g. 
form factors being too soft [13, 10], the missing intermediate range attraction in the nucleon- 
nucleon force [14] or the linearly rising phase shifts in pion nucleon scattering [11], have found 
satisfactory solutions. Partly, this was achieved by using effective meson Lagrangians [15, 16] 
describing the meson physics better than the original simple Skyrme Lagrangian. On the other 
hand, this increased complexity also generated more ambiguity in the effective action. 

At that point naturally the question arises whether a more microscopic picture of the soliton 
can provide some guide for choosing the effective meson theory. Even more, one may aim at 
a justification of the soliton picture of baryons in general. Therefore not only a microscopic 
realization of the soliton picture in terms of quark degrees of freedom is wanted but the 
dynamics of the model should be able, at least in principal, to determine its favored picture of 
the baryon. In this sense the Nambu-Jona-Lasinio (NJL) model [17] is unique. First of all, it 
is simple enough that such complicated field configuration like solitons can be determined self- 
consistently, see Chapter 6. Second, it contains the correct chiral symmetry breaking pattern 
and reproduces a lot of meson properties like masses, decay constants, scattering lengths etc., 
using only a few input parameters, see e.g. [18, 19] and references therein. Third, and more 
important, it allows for two complementary pictures of baryons, namely either as ordinary 
three-valence-quark bound state or as chiral soliton. The latter will be the subject of this 
review. In the other approach the baryon wave function is obtained as a solution of a Faddeev 
equation using diquarks as intermediate building blocks [20, 21]. Using functional integral 
techniques one can derive a generating functional which allows one to treat both pictures in 
one quark theory without any ambiguities or double counting [20]. Such a theory contains the 
relativistic quark model and the Skyrmion as asymptotic limits. Preliminary results indicate 
that such a hybrid model can display unexpected features, e.g. the diquark mass is drastically 
reduced in a soliton background [22]. 

Obviously, the simplicity of the NJL model which makes it suitable for such complicated 
investigations is also its most severe drawback. The NJL model is non-renormalizable and 
only uniquely defined if the necessary regularization prescription is specified. This introduces 
an ultraviolet cut-off which indicates the range of applicability of the model. The results 
for several observables differ in various regularization schemes [18]. Even worse, sometimes 
the qualitative behavior changes when altering the regularization prescription. Fortunately, 
the situation is not as dramatic for most of the observables, or can be understood from 
the deficiencies of the regularization procedure like e.g. missing gauge invariance. The other 
serious disadvantage of the NJL model is the absence of confinement, or more precisely, the 
appearance of two-quark (or quark-antiquark) thresholds. These cause unphysical imaginary 
parts in correlation functions for large (time-like) momenta. Therefore the results of the NJL 
model are restricted to low energies not only by the ultraviolet cut-off introduced via the 
regularization but also by the two-quark thresholds. 

This review is devoted to the soliton description of baryons within the NJL model. It is 
not the primary goal of such investigations to reproduce all the phenomenological successes 
of Skyrme type models. The basic motivation has rather been, and is still, to increase our 
understanding how the soliton emerges, and to better understand its relation to the underlying 
quark dynamics. Especially, the conceptional easy and direct access to the quark degrees of 
freedom allows one to study questions which escape our considerations when starting from 
a purely mesonic model. On the other hand, in order to arrive at physical baryons one has 
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to use techniques which are well known from the Skyrmion: semiclassical quantization [9] 
(cranking [23]), calculation of quantum corrections [24, 25], generalization to three flavors by 
either collective quantization [26, 27, 28] or bound state approach [29, 30], and so on. These 
subjects will be discussed in detail in this review, see Chapter 7. They should be considered 
as a basis for model calculations where one wants to describe baryons as chiral solitons and 
additionally wants to have access to the quark degrees of freedom in one consistent frame. 

The organization of this review is as follows: In Chapter 2 we will present some rigorous 
results from QCD which are relevant for the soliton description of baryons. After a short 
summary of the large Nc arguments for baryons being solitons we will outline some low- 
energy theorems based on chiral symmetry and the chiral anomaly. In Chapter 3 we will 
describe the NJL model and its bosonization [31]. In Section 3.3 a very important property 
of the NJL model, the dynamical breaking of chiral symmetry, is described. We will also use 
a local chiral rotation to display the hidden gauge symmetry [16] of the bosonized model [32]. 
Chapter 4 is devoted to the effective meson theory of the NJL model. Hereby, the gradient 
expansion does not only yield approximate meson masses and coupling constants but also 
reveals the relation to the Skyrme model. The determination of meson masses with the help 
of Bethe-Salpeter equations is explained in Section 4.3. Section 4.4 contains some comments 
on the chiral anomaly in the effective meson theory. In Chapter 5 the generic aspects of 
chiral solitons are discussed: the topological properties, the emergence of the soliton and its 
semiclassical quantization. Chapter 6 is the central piece of this review. After giving the 
energy functional of the static NJL soliton we discuss the self-consistent solutions for different 
meson fields included (or neglected). Chapter 7 embodies the description of baryons as NJL 
solitons. This also includes a comprehensive explanation of the cranking method for two and 
three flavors as well as the description of time-dependent meson fluctuations off the static 
soliton. In Chapter 8 we give a short summary. Some lengthy formulas which are nevertheless 
necessary to make this review reasonably self-contained are given in five appendices. 
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2 Rigorous results from QCD 



In this chapter we will present some rigorous QCD results which are relevant to hadron 
physics and in particular to the soliton description of baryons. As discussed in the introduction 
QCD cannot be treated perturbatively at low energies. As a consequence there are very few 
such results. They are either based on considerations assuming the number of colors Nc to 
be large or on the chiral symmetry of massless QCD and its spontaneous breaking. As we will 
see both provide substantial arguments for the chiral soliton picture of baryons. 

2.1 Large N c QCD 

In the low-energy regime there is no obvious expansion parameter to treat QCD pertur- 
batively. However, in the seventies 't Hooft [6] and Witten [7] demonstrated that generalizing 
QCD from the gauge group SU(3) to Sl^Ac) with Nc being large, 1/Nc might be consid- 
ered as an implicit expansion parameter. Naturally, one might wonder whether this provides 
a sound basis for a perturbative analysis since Nc = 3 in the real world. Nevertheless, we 
will see that this idea is very fruitful, especially it is the essential motivation for identifying 
baryons with solitons of meson fields. 

Assuming confinement in the large Nc world it can be shown that QCD with Nc — > oo 
and g 2 Nc fixed has a limit which can be described as keeping only planar Feynman diagrams 
[6]. The reason is that non-planar diagrams or diagrams with gluon handles are suppressed 
by Nq 2 , and the ones with quark loops by A^ 1 . Simple power counting arguments using 't 
Hooft 's diagrammatic analysis reveal that the correlation function of a color singlet operator 
is dominated by planar diagrams with exclusive gluon insertions and a quark loop on the 
exterior. Furthermore, in leading order in the 1/Nc expansion this correlation function cannot 
factorize, i.e. it is saturated only by color singlet intermediate states. This then leads to 't 
Hooft's most important result on large Nc QCD: In this limit QCD reduces to a theory of 
(infinitely many) weakly interacting mesons and glueballs. Their masses are of order N c , their 
mutial interaction is suppressed by powers of 1/Nc- 

The mesons are stable objects in the limit Nc — > oo. Their decay rates are of order 
1/Nc, their cross sections of order A^ 2 . This is very gratifying because it provides a natural 
explanation for the existence of narrow meson resonances in nature. Indeed, it is not obvious 
at all that e.g. the decay — > 27r is narrow enough to detect the p meson experimentally. 
Even more, the ratio of this decay width to the p meson mass very roughly follows the 1/Nc 
counting. To be more specific we consider an n-point correlation function of a color singlet 
operator. The summation of planar diagrams leads, among others, to an effective meson 
self-coupling which at tree level takes the form [7] 

C n = Nc* +1 \ il ... in * il ...* in . (2.1) 

Rescaling the meson fields $ 8 according to $ 8 = ^/Nc'-fi leads to 

C n = N c \ n ... tn tpn ■ ■ ■ (2-2) 

From eq. (2.2) one can then conclude that the generating functional of QCD becomes in 
leading order of the 1 /Nc expansion 

Z QCD = Jv[A,q,q]e iA *°° J V<pe iN ° A M (2.3) 
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where A[(fi] is an effective meson action that involves infinitely many mesons and is of the 
generic form (summation over repeated indices is assumed) 

A[<p] = j d 4 x fydtfi) 2 - ^rriijipupj + ^\ %3k ^^ 3 <^ k + . . )j . (2.4) 

1/Nc corrections are then given in terms of meson loop terms. 

Using the large Nc diagrammatics Witten conjectured [7] that baryons behave as soliton 
solutions of the effective theory (2.4). Assuming again that quarks are confined and that g 2 Nc 
is a constant the mass of a baryon composed of Nc quarks in their antisymmetric ground state 
is given by 

M B = N c (m q + T) + l -N c (N c - l)g 2 V « Ncf(g 2 N c , m q ) (2.5) 

to leading order in the 1/Nc expansion. Besides the masses m q and the kinetic energy T of 
the individual quarks also the one-gluon-exchange is included. By combinatoric arguments it 
can be shown that the n-gluon-exchange behaves like Nc(Ncg 2 ) n for Nc 3> n. Therefore the 
function J r (g 2 Nc } m) whose leading order is N c parametrizes the 1/Nc corrections in a smooth 
fashion. Since the strength of the meson coupling is proportional to 1/Nc the baryon mass 
Mb is proportional to the inverse of this meson coupling. As the baryon radius is proportional 
to the inverse of the quark kinetic energy up to suppressed binding effects it is independent 
of Nc in leading order. This is the typical behavior of a soliton field! These results (as 
well as the fact that the meson-baryon and baryon-baryon interactions are of order N c and 
N c , respectively) are the basis for Witten's conjecture that baryons emerge as solitons of an 
effective meson theory. 

Obviously, such a behavior can never be obtained in perturbation theory. Unfortunately, 
the effective meson theory (2.4) cannot be constructed explicitly." Nevertheless, this kind of 
reasoning provides some insight in the wealth of experimental results. The large Nc arguments 
account for suppression of exotics, confirms Zweig's rule and Regge phenomenology. Further- 
more, it is consistent with phenomenological models like e.g. the meson exchange models. 

On the other hand, to use the 1/Nc expansion also quantitatively one has to resort to 
effective models. The most prominent example for soliton models of baryons is the Skyrme 
model [8, 9], for a recent review see ref. [12]. In the following chapters we will present the 
description of baryons as solitons within the NJL model. Besides others it has the virtue that 
Witten's conjecture may be tested by the dynamics of the model. Nevertheless, one might 
also raise the question whether it is not more natural to describe baryons as three-quark 
bound states if one starts from an effective quark theory. Also here the 1/Nc expansion 
might give some hints. Starting from an effective quark interaction described by a current- 
current coupling of color octet quark currents (see eq. (3.2)) fierzing into attractive channels 
and subsequent functional integral hadronization leads to an effective theory which contains 
besides mesons also baryon fields composed of a diquark and a quark [20]. The interesting 
fact now is that in this effective theory the quark-quark, which is responsible for diquark 
formation, interaction is suppressed by 1/Nc as compared to the antiquark-quark interaction, 
i.e. in the limit Nc — > oo these explicit baryon fields are removed from the theory [35] and one 
has to find the baryons in the left-over effective meson theory which, by the way, is formally 
identical to the effective meson theory obtained by bosonizing the NJL model. This effective 
theory will be discussed in some detail in chapter 4, and its soliton solutions are exactly the 
NJL solitons. 



a In contrast to the situation in four space-time dimensions QCD2 can be bosonized exactly [33, 34]. 
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2.2 



Chiral symmetry and low— energy theorems 



Besides the exact color symmetry, the QCD Lagrangian possesses an additional symmetry 
in the limit of vanishing quark current masses, the chiral symmetry. In order to display this 
symmetry it is convenient to split the quark fields in right- and left-handed components. 
These are defined by 

qn = Pnq, qL = Piq (2.6) 

where 

Pr,l = \(1±1s) (2-7) 

are the corresponding projectors. Then the QCD Lagrangian taken at zero current masses 
decouples into a sum of two Lagrangians containing only right-(left-) handed fields, respec- 
tively. These two Lagrangians are invariant under global unitary flavor transformations of 
the corresponding right-(left-) handed fields, i.e. the QCD Lagrangian is invariant under 
Ul(Nj) X Ur(Nj) where Nj is the number of flavors under consideration. The decomposition 
into semi-simple subgroups 

U L (N f ) x U R (N f ) = U L+R (1) x U L - R (1) x SU L (N f ) x SU R (N f ) 

allows one to discuss the associated conserved charges. The invariance under Ul+r(1) is 
responsible for the conservation of baryon number whereas Ul-r(1) is subject to an anomaly 
[36]. The implications of this anomaly will be discussed in the subsequent section in more 
detail. 

The invariance under unitary transformations belonging to the subgroup SUl(Nj) xSUr(Nj) 
implies the conservation of the iV| — 1 currents 

A a 

Jl»( x ) = qL{ x )^i^qL{x) 

A a 

Jr»( x ) = qR( x )^i^qR{ x ) (2-8) 

where the matrices A a are the generators of SU(Nj). The corresponding conserved charges 

X 

J d 3 xJ a R0 (x) (2.9) 

fulfill commutation relations which are identical to the one of the generators of the group 
SUl(Nj) X SUr(Nj). However, they mix under parity, VQlP~ x = Qr. Parity eigenstates are 
given by the currents which are related to the diagonal subgroup SUl+r(Nj) and the coset 
(SU L (N f ) x SU R (N f ))/SU L+R (Nf): 

V;(x) = J a R ,(x) + J a L ,(x) = q(x)^q(x) 

A;(x) = J a R,(x)-J a L ,(x) = q(x)^ 75 q(x). (2.10) 
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Under parity these currents transform like vectors or axial- vectors, respectively. Furthermore, 
their time components fulfill the equal-time commutation relations 

[V a (x,t),V b (y,t)} = if abc V c (x,t)S^(x-y) 

[V a (^t) } A b (y } t)] = zr bc A c (x } t)6^(x-y) 

[A a (x,t),A b (y,t)} = if abc V c (x,t)S( 3 \x-y) (2.11) 

where f abc are the structure constants of the Lie algebra SU (Nj). These commutation relations 
are known as current algebra [37]. As the left hand side is quadratic in the currents whereas 
the right hand side is linear these relations also fix the normalization of the currents. This is 
the basis of the phenomenologically successful current algebra sum rules. 

Chiral symmetry is explicitly broken by the current masses. Nevertheless it is a good 
approximation for the two light flavors (up and down) since their current masses are of the 
order of a few MeV. For the strange quark whose current mass is of the same order of magnitude 
as the fundamental QCD scale, Aqcd, chiral symmetry is certainly not as good a symmetry. 
In spite of this, assuming approximate chiral symmetry for strange quarks is useful for a 
wide range of applications. For the heavy quarks the opposite type of expansion, namely 
in powers of 1/m, has proven to be successful [38]. Heavy quarks can be neglected in low 
energy processes as a consequence of the Appelquist-Carrazzone theorem [39] which states 
that particles decouple whose mass is much larger than the typical energy scale for the process 
under consideration. We will therefore treat only the cases of two or three flavors in the 
following. 

Chiral symmetry is spontaneously broken implying the existence of pseudoscalar (would- 
be) Goldstone bosons. These are identified with the isotriplett of pions or for three flavors 
with the octet of pseudoscalar mesons 7r, K and r\. A further consequence is the dynamical 
generation of a non-perturbative quark mass. It is expected that this quark mass, which is 
often labeled as constituent quark mass, is of the order of Aqcd- Furthermore, a finite non- 
perturbative quark mass implies a non-zero quark condensate which is defined in terms of the 
full quark propagator 

(qq) = —i lim trSi?(:c, y). (2.12) 



y- 



Note that the quark condensate (qq) is a gauge invariant quantity. Therefore it can be evalu- 
ated in any gauge. In a covariant gauge the quark propagator is of the form 

i 

SFiq) = 4A(q*) - B(q*) + ie ^ 

In perturbation theory in the chiral limit one has B(q 2 ) = 0. Therefore the quark condensate 
would vanish due to the vanishing Dirac trace of 7^. On the other hand, a dynamically 
generated B(q 2 ) ^ implies a non-vanishing condensate. 

Using (qq) ^ and the Goldstone theorem one can prove that the axial currents A^ (2.10) 
couple the Goldstone bosons to the vacuum. Denoting the one particle states of Goldstone 
bosons with momentum p by \Tr a (p)) one obtains 

(0\A;(x)\7r b (p))=zr b p,e-^. (2.14) 

The f ab are non-vanishing constants. In the isospin symmetric case they are proportional to 
the unit matrix, f ab = S ab f a . One usually determines their numerical values from weak pion 
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(kaon) decay, e.g. it — > \iv. f v (fx) lS therefore called pion (kaon) decay constant. For a 
compilation of recent experimental data see page 1443 of ref. [40]. For the purpose of this 
review it suffices to know that /„. = 93MeV and fx/fir = 1-22. 

Acting with the derivative operator on eq. (2.14) and using the Klein-Gordon equation for 
the pion state (p 2 = m 2 ) yields 

(0\d»A;(x)\Ap)) = 6 ah Umle-*** . (2.15) 

The conservation of the axial current implies either / T = or m v = 0. These two possibilities 
correspond to the Wigner-Weyl or Nambu-Goldstone realization of chiral symmetry, respec- 
tively. As both quantities are finite chiral symmetry has to be explicitly broken, i.e. one has 
to use finite current quark masses in order to describe nature. To further proceed eq. (2.15) 
is elevated to an operator identity. Let us first introduce the pion field operator <f>^(x) and 
choose its normalization with respect to the one pion state to be 

(0mx)\Ap)) =S ab e~^. (2.16) 

Then the operator identity 

d"Al(x) = f«mMx) (2.17) 

implies the relation (2. 15). 6 Eq. (2.17) is known as the PCAC (Partially Conserved Axialvector 
Current) hypothesis. Assuming additionally that meson and baryon form factors can be ex- 
trapolated smoothly from the corresponding mass shells allows one to relate different hadronic 
observables depending on both, weak and strong interaction parameters. One famous example 
is the Goldberger-Treiman relation 

f^g^NN = m N g A (2.18) 

which connects the axial coupling of the nucleon §a with the pion nucleon coupling constant 
9-kNN (i^in is the nucleon mass). This relation is experimentally fulfilled within ten percent. 
This (in)accuracy sheds some light on the usefulness of the PCAC hypothesis. 
Using the PCAC hypothesis one furthermore deduces [5] 

6 ah mlfi = i I J 4 x(0|%o)[A»(x),^(0)]|0) (2.19) 



in the limit of vanishing meson momentum. This may be rewritten in terms of the (partially 
conserved) axial charges Ql and the Hamiltonian 7i(x) as 

S ab mlf n = (0[Q a 5 ,[Ql,H(0)]}\0). (2.20) 

As these charges generate infinitesimal axial transformations it is obvious that only chiral 
symmetry breaking terms in the Hamiltonian contribute to the double commutator. In the 
case of three flavors the symmetry breaking mass term is given by (A = \/2/3 1) 



— 7 7 0— _A^ _A^ , , 

m u uu + m d dd + m s ss = c q—q + c 3 q—q + c s q—q (2.21 J 



'Note that eq. (2.17) is not derived from (2.15). One has also to assume that the matrix elements 
(Old 1 * A^(x)\X) for all states X except the one pion state vanish. 



12 



where the symmetry breaking parameters c 4 - are functions of the current quark masses 



1 / | | 0\ 

c = -^{m u + m d + m s ) 

1/0 0\ 

c 3 = ^yn u -m d ) 

c 8 = -l=( m o + m °-2m°). (2.22) 

The commutators (2.20) may be computed with the help of the explicit expressions (2.9) and 
(2.10) of the axial charges. Using the algebra of the Gell-Mann matrices and especially the 
relation 

A a A^ A c 

[<Zy7o759,<?y<?]= -id abc q—-f 5 q (2.23) 

(the d abc are the symmetric structure constants of SU(3)) one obtains the famous Gell-Mann- 
Renner-Oakes relations [41] 

flml = ±(m° u + m o d )(0\uu + dd\0) 

fW K = \(m O u + m O s )(0\uu + ss\0) 

f 2 v m 2 v = ^(m° + m° d )(0\uu + dd\0) + ^m° s (0\ss\0). (2.24) 

For the case of equal quark condensates 

(0\uu\0) = (0\dd\0) = (0|ss|0) 

the definition of the decay constants (2.14) implies that all decay constants are equal. For this 
special case eqs. (2.24) also yield the Gell-Mann-Okubo mass relation [42] \.m 2 K = 3m 2 + m 2 
which has been observed empirically. In addition one can extract the ratio of current quark 
masses 

— - = r— ^ w — . (2.25) 

2m° 2m 2 K - ml 25 K J 

Taking into account isospin breaking m° ^ m° and electromagnetic corrections the rela- 
tions (2.24) can be refined. Nowadays the values of all current masses can be estimated quite 
reliably. However, only the ratios of current masses are uniquely defined because these are 
scale invariant. Citing an absolute value of a current mass hence necessitates reference to a 
scale. At 1 GeV commonly accepted values for the current masses are [43] 



lGeV) p* 5MeV 



m° d (lGeV) w 9MeV 

m°(lGeV) w 160MeV. (2.26) 

One should, however, note that these data are to some extend model dependent because they 
are extracted form low-energy meson phenomenology [44] . 
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2.3 



Chiral anomaly and the QCD vacuum 



The properties of the QCD vacuum are still poorly understood. From hadron phenomenol- 
ogy one deduces that non- vanishing condensates have to exist. Examples are the gluon and 
the quark condensates, (G a G a,lv ) and (qq). The relation of the latter to the spontaneous 
breaking of chiral symmetry has been discussed in the preceding section. In this section we 
will look at it from another point of view: The quark condensate is related to the mean density 
of eigenvalues of the quark Dirac operator. To see this we consider the quark propagator for 
a fixed gauge field A^ xed 

(S F (x, y)) A ^ d = (Tq(x)q(y)) A ^ = £ "^""M (2.27) 

where u n (x) and X n are eigenfunctions and eigenvalues of the Euclidean Dirac operator 

pu n (x) = X n u n {x). (2.28) 

Except for the zero modes the eigenfunctions occur in pairs of opposite chirality with corre- 
sponding eigenvalues iX n . Therefore one obtains from eq. (2.27) 

{ d A x{q{x)q{x)) Al ^ = - 2 ^- , \2 ( 2 - 29 ) 

A„>0 



V Jv y " v V ^ n (m°) 2 + X 2 



where the zero-mode contribution has been neglected. To arrive at the quark condensate one 
has to average over all gauge field configurations and then take the limit V — > oo. In this 
limit the spectrum becomes dense and gives a non-vanishing condensate if the mean number 
of eigenvalues in the interval dX is proportional to the volume: 

^ = - 2 -° r dx jj^h (2j °> 

where p(X) is the mean spectral density. Taking now the limit m° — > one arrives at [45] 

(qq) = -wp(0). (2.31) 

i.e. the quark condensate is related to the level density at zero virtuality, A = 0. 

In deriving this result the infinite volume limit is very important. In a finite volume there 
is no spontaneous breaking of chiral symmetry, and therefore the quark condensate would 
vanish. The spectrum is discrete for finite V and the sum in eq. (2.29) does not develop an 
infrared singularity. If the chiral limit had been taken at finite V the chiral symmetry would 
have been restored. Obviously, the two limits m° — > and V — > oo are not interchangeable. 

In order to analyze the situation at finite volume one has to consider the dimensionless 
quantity x := —Vm°(qq). (Eq. (2.31) implies that the level spacing for small A is AA 
ff/V(qq).) Hence the denominator (m ) 2 + A 2 varies only slowly for neighboring levels if 
x ^> 1. In this case it is a good approximation to replace the sum in eq. (2.29) by an integral, 
i.e. eq. (2.30) stays true as long as the quark current mass is larger than 1/V(qq). This suggests 
that spontaneous chiral symmetry breaking is related to the appearance of small eigenvalues 
of ^)such that X n oc l/V. 

Recently the spectrum of this operator has been investigated [46] by studying the role 
played by the winding number v of the vacuum gauge field configuration. A topologically 
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non-trivial gluon field configuration necessarily gives rise to quark zero modes. For small 
(or vanishing) quark masses these zero modes tend to suppress the fermion determinant like 
(m°)^ N f (Nf is the number of light flavors). However, as emphasized in ref. [46] this sup- 
pression of the winding number v configuration is always accompanied by an enhancement 
proportional to V^ Nf . Hence, in the physical situation x 3> 1 there is no suppression at all. 

A further very interesting result obtained in ref. [46] is a set of sum rules for the spectrum 
of the Dirac operator p. For a gluon field configuration with winding number v there are 
\v\ zero modes of p. The few lowest non-zero eigenvalues are proportional to 1/V(qq) as 
anticipated from the discussion above. Their distribution is sensitive to the winding number, 
and the levels are pushed up if \v\ increases. The sum rules of ref. [46] relate the inverse 
moments of the eigenvalue distribution to the quark condensate. E.g. the lowest one is 



K 4 W\ + N f 



These sum rules reflect the fact that for a finite volume the eigenvalues with —1/V(qq) ^ \ n <C 
Aqcd are the one related to spontaneous chiral symmetry breaking and the occurrence of the 
quark condensate. 

On the basis of these results one can conclude that the winding number is irrelevant as 
long as a; > 1 which is most likely the case in the real world. However, the winding number 
density fluctuations, i.e. the topological susceptibility 



X = V = / d 4 x(GG(x)GG(0)), (2.33) 



V (32tt 2 

has measurable consequences for hadron physics: it generates the rj' mass. Phrased otherwise, 
despite the fact that the winding number is irrelevant the axial U(l) symmetry is still broken 
in an anomalous fashion, and this is reflected by the QCD vacuum. The strength of this 
breaking is determined by the topological susceptibility \- Also here Leutwyler and Smilga 
[46] found an astonishing result, 

— Tm°(qq) 

x = M o^~\~ ( 2 - 34 ) 

iVjT — m u (qq) 

where r is the would-be topological susceptibility in the absence of quark fields. If one 
considers large Nc and a fixed quark mass, r has a finite limit whereas (qq) oc Nc- Then one 
obtains x = r, i.e. in the large Nc world the topological susceptibility would be given only by 
gluons. However, for Nc = 3 and a small current quark mass (m°(qq) <^ r) the topological 
susceptibility is dominated by the quark condensate: \ = —m°(qq) /Nj (for unequal quark 
masses m°/Nf is replaced by the reduced mass). The mean square winding number is large, 
(z/ 2 ) = —Vm°(qq)/Nf } and for an infinite volume all winding numbers are equally likely. 

The physical implications of these considerations become more transparent when one con- 
siders the anomalous Ward identity 

Nf ~ 

dS = 2im°j 6 - ^- 2 G%G a ^ (2-35) 

where 

Js = (97^759) and j 5 = (q-f 5 q) (2.36) 



c For the simpler case of an abelian anomaly the derivation of the anomalous Ward identity is presented in 
Appendix A. 
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denote the axial singlet and pseudo-scalar currents of the quarks, respectively. Considering the 
correlator of the axial singlet current this anomalous Ward identity may be used to derive an 
expression for the rj' mass [47] 



Since for Nc = 3 the topological susceptibility x i s governed by the quark condensate one 
might have anticipated this to be the case for the rj' mass, too, thereby contradicting the result 
of ref. [47]. Note, however, that the rj' mass is still dominated by the winding number density 
fluctuations of the purely gluonic theory. 

Summarizing this chapter we would like to emphasize the following points: First, consid- 
ering the limit of a large number of colors Nc indicates that baryons may be described as 
solitons (Witten's conjecture). Unfortunately, the explicit expression of the effective meson 
theory (2.4) is unknown. Second, from the approximate chiral symmetry of QCD and its 
spontaneous breaking we know that low energy hadron physics is dominated by the would-be 
Goldstone bosons of chiral symmetry, the pions. Thus, at least, these meson fields have to 
appear in an ansatz for the effective meson theory. Third, due to the chiral anomaly in the 
flavor singlet channel we know that the chiral symmetry breaking is an inherent property of 
the small eigenvalues of the quark Dirac operator in the QCD vacuum. 

These ideas will at least partially be applied in the proceeding chapters. We will present a 
model which is meant as a simplified version of the low-energy quark dynamics of QCD, the 
NJL model [17]. It displays dynamical breaking of chiral symmetry and yields quite a successful 
description of meson physics. Having noticed this common feature with the underlying theory 
we will again take advantage of the large Nc arguments and attempt a description of baryons 
as solitons in the NJL model. 
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3 The Nambu— Jona— Lasinio model 



In this chapter we introduce the N JL model as an effective chirally invariant theory ol quark 
flavor dynamics. Originally, it was proposed to describe the pion as a massless bound state of 
the nucleon and the anti-nucleon [17]. Nambu and Jona-Lasinio studied a local four-fermion 
interaction in the scalar-isoscalar and pseudoscalar-isovector channel. By construction the 
original NJL model possesses a global U(l) X SU(2)^X SU(2)# symmetry. Nowadays, it is 
common to call all chirally invariant models with local four-fermion (or even six-fermion) 
interactions an NJL model. 

3.1 Description of the model 

To be specific we will consider the model described by the Lagrangian 



£njl = q(i$ - rh°)q + 2d J2 [ i^^f + (?y*75?) 

8 = * 



Nj-1 



~ 2G 2 E [iqjl.q) 2 + (q^lsl.qfj ■ (3-1) 
Here q denotes the quark spinors and m° the current quark mass matrix. The matrices 



A 8 /2 are the generators of the flavor group (A = y 2/NfE), Nf being the number of flavors 
under consideration. Note that the coupling constants G\ and G2 have dimension [energy] -2 . 
These coupling constants may take different values, G\ 7^ G2, without spoiling the global 
chiral Ui(iV/) X \]r(Nj) symmetry, i.e. the two sums in eq (3.1) are independently chirally 
invariant. 

The special form of the Lagrangian (3.1) can be motivated from the one-gluon-exchange. 
In the local limit of the one-gluon-exchange the quark interaction is given by the current- 
current interaction 

3*M (3-2) 
where is the color octet flavor singlet current of the quarks, 

j^q^q, a=l,...,iV c 2 -l = 8, (3.3) 

the A"/2 being the generators of the color group SU(3) in the fundamental representation. 
Fierzing the interaction (3.2) into the color singlet channel leads to the interaction of eq. (3.1) 
with the additional relation G2 = G\/2 [20]. Since we are considering the NJL model as an 
effective theory we will relax this condition. As mentioned in section 2.1 the interaction (3.2) 
also contains quark-quark interactions leading to diquarks. These are, however, suppressed 
for large Nc [35]. 

The invariance of the Lagrangian (3.1) under chiral rotations 

q^U v q, q^qU v , U v = exp (i9 v ) , V = 9 l y — 

X 

q -> U A q, q -> qU A , U A =exp(i9 A ), 9 A = 9 A —-f 5 (3.4) 
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in the limit m° — > is easily verified. Expressing the quark fields in terms of right- and 
left-handed quark fields defined by (c/. section 2.2) 

qn = Pnq, qL = Piq (3.5) 

the Lagrangian (3.1) decouples in a sum of two Lagrangians containing only right-(left-) 
handed fields, respectively. Only the mass term 

£ mass = -qm°q = -(q R m°q L + q L rh°q R ) (3.6) 

couples right- to left-handed fields and thus breaks the chiral symmetry. 

The local four-fermion interaction of (3.1) with a dimensionful coupling constant is obvi- 
ously not renormalizable. Therefore it is only completely defined when supplemented with a 
regularization scheme in order to cut off momentum integrals and thus avoiding ultraviolet 
divergencies. In QCD these divergencies would have been absorbed in the renormalization 
procedure. On an operational level one may interpret the occurrence of the cut-off as a very 
crude way of mimicing the asymptotic freedom of QCD. Various schemes have been discussed 
in the literature": 0(3) and 0(4) invariant sharp cut-offs, Pauli-Villars regularization and so 
on. Here we will exclusively use the proper time regularization scheme proposed by Schwinger 
[48]. This procedure has the advantage of being gauge invariant if (external) gauge fields are 
coupled to the model. Especially, for the calculation of electromagnetic and weak form factors 
this scheme is superior to sharp cut-offs because it allows one to appropriately manipulate 
the momentum space integrals. A further benefit of this regularization scheme is the fact that 
it may be defined at the level of the action rather than being defined via its application to 
Feynman integrals. This automatically guarantees that different quantities are regularized in 
a consistent manner. We will apply this regularization to the bosonized action, see below. 



3.2 Bosonization 

The aim is to rewrite the quark (fermion) theory (3.1) into an effective meson (boson) 
theory. It is convenient to use a compact notation 

X 

h a = — ®v a , i = o,...,]^ 2 - i, r a e {1^75^,^75} (3.7) 



an 



d 



)aP= UG 1 6 a ^ for T a G {l,z 75 } 

\4G 2 S^ for T a G {7^75} ' 



(3.8) 



The auxiliary field $ = $ a A a , which is introduced via the identity, 

exp f qA a qQ afJ qA^ = f D$ exp f ^(Q" 1 )^ - i f <S> a qA a q^J , (3.9) 

contains therefore (in the case of three flavors) nonets of scalar, pseudoscalar, vector and 
axial-vector meson fields. Using eq. (3.9) the generating functional 



Znjl = J T>qVqexp(—i J d 4 xC 



NJL 



"For a recent review see ref.[18] 



18 



may be written as [31] 

Z NJL = J X>$exp J QQ- 1 ^ Z F [$], 

Z F [$] = J VqVqexp (-i j q(i$ - m° - $A)g) . (3.10) 
For the interpretation of Z F [$] we note that it is equivalent to 

Z F [$] = lim (0|e^ T |0) (3.11) 

1 -^oo 

where T denotes a large Euclidean time interval and h = J d 3 xq^hq is the second quantized 
form of an one-particle Dirac Hamiltonian. In Minkowski space this Hamiltonian is given by 

h = a •p + /3(m° + $) , p = -d, (3.12) 

i 

as can be verified from the relation i(j) — rh° — $ = f3(id t — h). 

For subsequent considerations it is convenient to remove the current quark mass from the 
Dirac operator by shifting $ — > $ — m°. This yields the functional 

z NJL = /d$ e -i/^($-^)g-i($-^) r VqVq - e -ifd**M-*)^ (313) 



Furthermore we decompose the generic meson field $ into irreducible Lorentz tensors 

$ = S + i^P-if-i^. (3.14) 

S is a scalar, P a pseudoscalar, V a vector and A an axial-vector field. All these fields are 
flavor matrices, i.e. S = S t (X t /2) etc. 

Using eq. (3.14) the chirally invariant interaction term of the NJL model is written as 

J($ - m )^ 1 ^ - m°) = -^tr((S - m ) 2 + P 2 ) + -Ltr^V + A^). (3.15) 

Z L\j\ Z(j2 

For subsequent considerations it is also convenient to introduce the angular decomposition of 
the scalar and pseudoscalar meson fields by defining a complex field M 

M = S + iP = dx£ R , (3.16) 

which in turn defines the Hermitian field E and unitary fields £l and £r. The decomposition 
(3.16) is not unique; £l and £r are rather related by a "gauge" condition. 

In eq. (3.13) the quark field appears bilinearly in the exponent and can therefore be inte- 
grated out. Using that Det = exp Tr log one obtains 

Z NJL = JvQe"™, 

A[$] = -^y > J 4 x($-m )g- 1 ($-m°) + Trlog(z^-$). (3.17) 

Note that the action A.[<&] is a non-linear, even non-polynomial function of the meson field $. 
Even more, the term Trlog(z^ — $) is non-local. The quantum theory defined by eq. (3.17) 
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is, however, equivalent to the underlying NJL model defined by the Lagrangian (3.1). On the 
other hand, the generating functional (3.17) has the advantage that it may be treated in a 
semiclassical approximation because the vacuum expectation value of the bosonic field $ can 
be different from zero, and in general will be, whereas the vacuum expectation value of the 
fermionic quark field is necessarily vanishing in the absence of external fermion sources. 
Symmetry currents are constructed by adding external gauge fields 

a (v) = a H«A«'/2 and 

= a^ l \ l /2 to the action for the vector (v) and axial-vector (a) symmetries, respectively 



q{i(j) -<$>)q^ q [i$ - i^' - - $J q. (3.18) 

The symmetry currents are then identified as the terms coupling linearly to the external gauge 
fields 



(3.19) 

a O,a) M = 

The gauge fields may be eliminated from the fermion part of the action by transforming the 
(axial) vector fields accordingly 

V^V^ + aW and A,^A, + a^. (3.20) 

This allows one to straightforwardly compute the derivatives in eq (3.19) yielding the current 
field identities [37] 

tf ) = 2^ and rt ) = W A »- (3 - 21) 



3.3 Dynamical breaking of chiral symmetry 

The vacuum expectation value (VEV) of the auxiliary field $ is found from the stationary 
point of the action (3.17): 

= -Q afJ {qA p q). (3.22) 

Especially, the quark condensate (qq) is related to the VEV of the scalar field via the cou- 
pling constant G\. Alternatively the effective action A[<&] (3.17) yields the Dyson-Schwinger 
equation 

® a {x) = -Q aP tv{G^{x, x)A p ). (3.23) 

The solution of this equation determines the VEV of the meson field $. Hereby G$ is the 
quark propagator in the background of the $-field. It is defined by 

Gl 1 {x,y) = {i$-<$>)6{x-y). (3.24) 

Before being able to derive explicit expressions for (3.23) in the NJL model a regularization 
scheme has to be imposed. Using the decomposition (3.14) the effective action (3.17) may be 
cast into the form 

A F = Trlog(z|)) = Trlog(z^ + zf+z^y 5 - (P R M + P L M^)), (3.25) 
Am = j d\ (-^-tr(ikftjkf - m°(M + Aft) + (m ) 2 ) + ^MV^ + A,A^ . 
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As this action is equivalent to the non-renormalizable NJL model it is only completely de- 
fined if a regularization scheme is provided. As already stated in subsection 3.1 we will use 
Schwinger's proper time regularization [48] which introduces an (9(4)-invariant cut-off A after 
continuation to Euclidean space. For this regularization procedure it is necessary to consider 
the real and imaginary part of Af separately 

A F = A r + Af 

Ar = ^Trlog(lilfe), 

Ai = ^Trlogaiir 1 ^). (3-26) 

The real part Ar diverges for large momenta p whereas the imaginary part Ai does not contain 
ultraviolet divergencies, i.e. it is finite without regularization. 6 Therefore one has the option 
of keeping Ai unregularized, or to regularize it in a way consistent with the regularization of 
Ar. Note that this defines two different models. 

For the real part of the action the proper time regularization consists in replacing the 
logarithm by a parameter integral 

1 Z" 00 ds f + \ 

Ar^--\ — Trexp {-sp E pE , (3.27) 

2 Jl/A 2 s \ J 

which for A — > oo reproduces the logarithm up to an irrelevant constant. Since the operator 

$ePe is Hermitian and positive definite this integral is well defined. 

For the issues discussed in this section it is sufficient to only inspect Ar. Varying the 
regularized effective action with respect to the scalar and pseudoscalar fields yields the Dyson- 
Schwinger or gap equations 



m t = m° l -2G 1 (qq) t 
No 
4tt 2 



(qq)i = -m^r(-l,m 8 2 /A 2 ) (3.28) 



which is the regularized version of eq. (3.23) specialized to the case of the scalar meson field. 
The dynamically generated constituent quark masses m 8 -, i = u } d } (and also the quark con- 
densates (qq)i ) are equal, m u = and (uu) = (dd) } respectively, if and only if the quark 
current masses are equal, m° = m°. Throughout this review we will restrict ourselves to this 
isospin symmetric limit, and in the following we will use the notation m := m u = and 
m° := m° = m°. Of course, for the strange quark we will consider larger current masses 
according to (2.25) leading to larger constituent masses. For these quantities we will keep the 
flavor index and denote them m° s and m s , respectively. 



'For a more detailed discussion of a correct definition of Ai and its properties see sect. 6.1. 
c The functions 

T(u, x) - 

are known as incomplete T-functions. Especially, 

T(0,x) = -logx + j + O(x) for x 0+ 
where j = 0.57721 is Euler's constant. 
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Solution of the gap equation 

A=630MeV 




0.0 10. 20.0 30.0 40.0 

G, A 2 

Figure 3.1: The solution of the gap equation (3.28) for vanishing current mass m° = (solid 
line) and m° = 18MeV (dashed line). For the chosen value of the cut-off, A = 630MeV, 
a constituent quark mass m = 400MeV reproduces the phenomenological value of the pion 
decay constant, /„. = 93MeV. 

As can be seen from figure 3.1 in the chiral limit the quark condensate and therefore also 
the quark constituent mass is zero when the coupling constant G\ stays below a critical value. 
Above this critical value the trivial solution coexists with a non-trivial one. Examining the 
effective potential for a constant scalar held in the chiral limit d 

V(E) = ^-E 2 + ^ (Vr(0, - (E 2 - A 2 )AV S2 / A j (3.29) 

one concludes that the non-trivial solution is energetically favored. 

3.4 Chiral rotation and hidden local symmetry 

The freedom in the choice of in the parametrization (3.16) of the scalar and pseu- 
doscalar held reflects the local hidden symmetry SU(3)^. Under SU(3)l X SU(3)# X 811(3)/! 
the fields and E transform as 

d The first part follows trivially from A m whereas the second part can be obtained from Ar by evaluating 
the functional trace. 
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£r(x) — > h(x)( R (x 

E(x) -> /i(a;)E(a;)/it(a;) (3.30) 

where 

X G SU(3) L ,i? G SU(3) fi and h G SU(3) /l . (3.31) 

The additional degrees ol Ireedom contained in ^l jR can be gauged away as in the usual Higgs 
mechanism [16]. Later on we will adopt the unitary gauge 

tR = d = Z (3-32) 

where £ is an element ol the coset-space SU(3)l X SU(3)^/SU(3)y . Gauge conditions like 
(3.32) lurthermore fix h(x) in terms ol L(x), R(x) and £(x). 

For subsequent considerations it is convenient to perform a chiral rotation Irom the current 
quarks qL,R to the constituent quarks (c/. eq (2.6)) 

Xl,r = ihMLfi- (3.33) 
Under a chiral rotation ol the original quark fields 

q L ->■ Lq L , q R ->■ Rq R (3.34) 
the constituent quark fields transform according to the hidden gauge symmetry 

Xl,r h(x)xL,R- (3.35) 

Alter the chiral rotation 

qipq = Xipx (3.36) 

the chirally rotated Dirac-operator 

ip=TipT^ with T = P R ( R + P L ( L (3.37) 

becomes 

ip=i(p+f+ A 75 )-E. (3.38) 
It contains the chirally rotated vector and axial-vector fields 

%TK = Cl,r(V, TA, + d,)d, R . (3.39) 

These fields transform under a lelt-right-transformation according to the hidden symmetry 
group 



V,TA,^ h(x)(V, T A, + d,)h\(x). (3.40) 

From these relations it follows that the vector field transforms as a gauge field ol the hidden 
local symmetry 811(3)/! defining a covariant derivative 

D^ = d^ + V^. (3.41) 
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On the other hand the chirally rotated axial-vector field transforms homogeneously under the 
hidden local symmetry 



A^(x) -> h(x)A^(x)h\(x). (3.42) 
In terms of the chirally rotated fields the Lagrangian becomes [31, 32] 



E 2 -mo(d% + d^L 



1 tT F \{V,-v,) 2 + {A,-a,) 2 }. (3.43) 



4G 2 ^ 

Here we have defined the vector and axial-vector fields 

= \(tndA + tLdA) (3-44) 

a, = liUdA-tLdA) (3-45) 

which are induced by the chiral rotation. The Lagrangian (3.43) is the basis for the investi- 
gation described in section 6.2.5. 

The chiral rotation of the quark fields does not leave the fermionic integration measure 
invariant but introduces a non-local Jacobian [49] 

J DqDq = J J D x Dx, (3.46) 

which contains the anomaly. Its physical implications will be discussed in section 4.4. 
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4 Effective meson theory 



In the previous chapter the NJL model has been converted into an effective meson theory. 
In the present chapter we will examine its implications for meson physics. 



4.1 



Gradient expansion 



The quanta of the small amplitude fluctuations of the composite meson fields $ (see section 
3.2) around their vacuum expectation values represent the mesons. For the extraction of these 
(free) mesons the effective action has to be expanded up to second order in the fields $. This 
will be done in section 4.3 and leads to the Bethe-Salpeter equation. At the moment we will 
concentrate on a simplified description of the mesons at low energies. In this section we will 
follow ref. [31]. 

The effective meson action obtained in the previous section by bosonizing the NJL model 
is a highly non-local object due to the presence of the quark determinant, Tr log(z^ — $). Since 

$] 7^ this term contains infinitely high derivatives of the meson field. Obviously, for the 
determination of the low-energy (e.g. static) properties of the mesons we do not need to keep 
all the higher order derivatives but a gradient expansion of the quark determinant seems to be 
appropriate. In leading order gradient expansion of the real part of the effective meson action 
one finds 



A 



R 



d 4 xl - C (m/A) + C 2 (m/A) [^(VM^M) - ^t(F^F v ^) - ^t(F^F a ^) 



1 



1 



1 



(4.1) 



where 



V'M = i) M + [V\ M] - 75 {A 1/ , M} 



(4.2) 



denotes the covariant derivative of the scalar - pseudoscalar field M. Furthermore, 



F, 



F, 



A 



d,V v 



d^ + r^KHK,^] 
-^A A1 + [A A1 ,K] + [y A1 ,A i ,] 



(4.3) 



are the vector and axial-vector parts of the field strength tensor. The quantity C^m/ 'A) 
depends only on the ratio ra/A and diverges logarithmically as A — > oo. Eq. (4.1) defines the 
gauged linear cr-model. To cast this model into its standard form one has to perform field 
renormalizations and redefinitions in order to remove the tt — A\ and a — V mixings from the 
Lagrangian. One then finds for the vector and axial-vector masses 



4G 2 



and m\ = m v + 6m^ 



(4.4) 



where 



N c ( m 2 
24^ I ' 



-1/2 



(4.5) 
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is the universal vector meson coupling constant. Furthermore the scalar meson mass is given 
by m 2 a = 4m 2 + m\ (see also eq. (4.46)). The pseudoscalar meson mass is given by 

ml = 77V (4-6) 

The resulting effective meson Lagrangian satisfies the well-known and phenomenologically 
very successful current algebra results". This comes with no surprise since our starting quark 
theory has strict chiral symmetry for m° = 0. From eqs (4.5) and (4.6) one furthermore finds 
the relation 



m 



2 _ 2 f 2 

V 



a= (4.7) 



which for a = 2 represents the KSRF-relation [50]. Furthermore, for this "magic" value of a 
Weinberg's relation [51] niA = y/2mv also holds true. All these low energy theorems are quite 
reasonably reproduced by the experimental data. 

So far we have discarded the imaginary part of the effective meson action (Ai in eq (3.26)). 
This part allows for a similar gradient expansion where again subsequent terms are suppressed 
by powers of the ratio ra/A. As already mentioned in chapter 3 the imaginary part of the 
effective meson action stays finite when the cut-off is removed. For A — > oo only the leading 
order term of the gradient expansion survives and this term is then given by the (gauged) 
Wess-Zumino-Witten action [52] (see section 4.4) 

Ai = Awzw + ■■■ (4.8) 

Therefore in leading order gradient expansion we find for the effective meson action on the 
chiral circle (E = (E) = m) 

A = A n i a + Awzw (4.9) 

where A n \ a denotes the gauged nonlinear cr-model displayed in eq. (4.1). This model La- 
grangian has stable chiral solitons provided that the (axial-) vector mesons are retained 6 . 



4.2 Relation to the Skyrme model 

For the soliton description of baryons the chiral field is crucial. Let us therefore concentrate 
now on the chiral field while setting all the other meson fields to their vacuum value. Aside 
from the pseudoscalar mass term the effective meson action of the chiral field is then given 
by the quark loop Trlog(z^ — m(U) 15 ) where U = ClCr 1S the chiral field, see eq. (3.16). In 
leading order its derivative expansion is given by the nonlinear cr-model Lagrangian 

C nla = -lftrL,L^ L, = U^8,U (4.10) 

in agreement with eq. (4.1). As is well-known this model does not allow for stable solitons 
as can be seen by applying Derrick's theorem [53]. To obtain stable solitons the non-linear 
cr-model has to be supplemented by stabilizing higher order derivative terms like the Skyrme 

"For a compilation of current algebra results see ref.[37] and references therein. 
'For a review on mesonic solitons with vector mesons see ref. [10]. 
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term [8], which has four derivatives. One would expect that these stabilizing terms emerge in 
higher order of the gradient expansion of the effective meson action. In next to leading order 
of the derivative expansion one finds in the limit A — > oo [31] 



The first term is in fact the desired Skyrme term which shows up with a coefficient e = 2tt. 
This value favorably compares with the phenomenological value e = 5.45 obtained by Adkins 
et al. [9] by fitting the Skyrmion mass (after semiclassical quantization, cf. section 7.1) to 
the nucleon mass. Unfortunately the last two terms destabilize the soliton. The last term 
leads to a negative definite energy for any chiral field configuration and forces the soliton to 
collapse. The second term which provides a tachyonic pole in the pion propagator is even 
more dangerous because it leads to a vacuum instability [54]. Even if one goes to the next 
(sixth) order gradient expansion one does not find stable solitons [55]. This certainly signals 
the break-down of the naive gradient expansion in the soliton sector of the effective action of 
the chiral field. These rather frustrating results had for some time given rise to the fear that 
the chiral soliton of the effective action might not at all posses soliton solutions. 

Some hints for the existence of stable chiral solitons are provided by the heat kernel expan- 
sion of the effective chiral action (3.25). The heat kernel expansion represents a semi-classical 
type of asymptotic expansion which corresponds to a resummation of the gradient expansion 
[56]. In next to leading order this expansion gives, besides the forth order terms given by 
equation (4. II), also a term with six derivatives [31] 



which in fact overcomes the forth order terms and stabilizes the chiral soliton [57]. 

Let us also mention that the next to leading order terms of the derivative expansion consid- 
erably improve the prediction for the 7T7r-scattering length towards their experimental values 
[58]. In fact the coefficients of these terms which are predicted by the gradient expansion of 
the quark loop are in rather good agreement [31] with the phenomenological values determined 
in chiral perturbation theory [59]. Thus the destabilizing terms are obviously relevant for the 
low energy meson physics. 

Nevertheless, we will show below that the Skyrme model does arise as low energy ap- 
proximation to the non-local effective meson theory provided one includes the vector and 
axialvector mesons in an appropriate manner. The importance of the vector mesons for the 
stabilization of the chiral soliton against the collapse should come with no surprise since the 
scalar and axialvector mesons are much heavier than the pseudoscalar mesons and therefore 
should control the short distance behavior of the soliton at low energies. 

For the study of the soliton sector of the effective meson theory it is obviously advantageous 
to regard the chirally rotated vector and axialvector fields introduced in section 3.4 as the 
physical fields. The corresponding treatment in the soliton sector will be discussed at length 
in subsection 6.2.5. Note that after the chiral rotation described in section 3.4 all higher than 
leading order gradient terms of the chiral field U are the fully absorbed into the rotated vector 
and axialvector fields. The elimination of these gradients is obviously crucial in The soliton 
sector since the gradient expansion seems not to converge for the chiral field, at least not in 
next-to-leading order and next-to-next-to-leading order. On the other hand for the fields of 

c Note that Atkins et al. considered f w as a free parameter. Their fit to baryon masses provided f w = 64MeV. 
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the heavier vector and axialvector mesons a gradient expansion might still be appropriate at 
low energies, which, however, has to be justified a posteriori. 

In terms of the rotated fields the leading order of the derivative expansion the vacuum part 
of the action is given by the following Lagrangian [32] 



wzw 



(4.13) 



where 



■ 1(7 



2g 



tr 



F v ) 2 + (Pa 2 



my 
9v 



tr 



V 



+ A 
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a — 1 

, 2 



f>A 2 



a 



+ -flraltv 



U + 



(4.14) 



Note that we have assumed the chiral circle condition E = rh. Cwzw is the Wess-Zumino- 
Witten [52, 15] term expressed now in terms of the rotated fields V^, A^. The vector fields 
Vfj, and have been defined in eqs (3.44,3.45). Furthermore the field strength tensors 
{Fv)fj,ut (Fa) hv are defined in equation (4.3) with the unrotated fields replaced by the rotated 
fields Vf^Af,. 

Since the (axial-) vector meson masses my,m^ are large compared to the pion mass it is 
tempting to integrate out the (axial-)vector fields V^, A^ in the above Lagrangian in the static 

limit. Then one might consider the kinetic energy ([Fy) 2 ^, (Fa) 2 ^ and higher order terms as 
well as the Wess-Zumino term as perturbations. The resulting equation of motion yield [32] 



K 



A„ 



A 1 1 

a — 1 



(4.15) 
(4.16) 



For these field configurations the field strength tensors reduce to 



■ I) [aVI 
([d^ + uV] + [a' 1 ,d v + v v ]) = 



(4.17) 
(4.18) 



and the Lagrangian (4.14) precisely becomes the Skyrme Lagrangian 
C = -\f>L,L» + g^tr^, L„] + \m 2 J 2 MU + 

with 



(4.19) 



9v 



\2a - II 



(4.20) 



Thus in the static limit of the vector and axial-vector mesons the total effective meson La- 
grangian in fact becomes the celebrated Skyrme model [8]. 
Since the lowest lying axial-vector mesons ( tti a ~ 

1260 MeV) is considerably heavier than 



the lowest lying vector mesons (trip ~ 
mesons at low energies. Then one obtains 



770MeV) one might discard the axial-vector 



(4.21) 
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£e = ~^lB,B» (4.22) 



Using the experimental value of g v ~ 6.0 [40] for the universal vector coupling constant, the 
Skyrme parameter is not far from its phenomenological value g v = 5.45 obtained by Adkins 
et al. [9] by fitting the nucleon and A-masses. 

In the derivation of the Skyrme model we have used the static limit for the vector and 
axial-vector mesons considering both the kinetic energy as well as the Wess-Zumino-Witten 
term as perturbations. When the static equations of motion (4.15,4.16) are used the Wess- 
Zumino-Witten term vanishes. Alternatively, one could, however, also include the Wess- 
Zumino-Witten term into the definition of the static equation of motion and leaving only 
the kinetic energy as a perturbation. This adds a current-current interaction to the Skyrme 
Lagrangian 

1 

— i 
2 

from the coupling of the cu-meson to the topological current in the Wess-Zumino-Witten 
term. In the context of the NJL model the coupling strength is obtained to be [60] 

2 _ 6vr 2 iVc 

66 " m2r(0,m 2 /A 2 )' (4 ' 23 ^ 

The current-current interaction (4.22) modifies the short distance behavior of the soliton 
resulting in a reasonable description of the short range behavior of the central potential in the 
nucleon-nucleon interaction [61]. 

Alternatively, if the (axial) vector mesons are kept as dynamical fields there will be no 
reason to neglect the Wess-Zumino-Witten action (contained in the full Lagrangian) in com- 
parison to the normal parity terms. In fact, the Wess-Zumino-Witten term is then absolutely 
necessary to ensure the existence of stable solitons. On the other hand, if one integrates out 
the (axial-) vector mesons in the static approximation in the way it was done above, the 
Wess-Zumino-Witten term vanishes for the corresponding static field configurations of two 
flavors. In this case stability of the soliton is provided by the induced Skyrme term. 

4.3 Bethe— Salpeter equations 

In this section we will discuss the determination of the meson masses as functions of the 
NJL model parameters. However, this time we will go beyond the gradient expansion, i.e. we 
will take into account the gradients of the fields in an exact manner. This can be done by 
extracting the Bethe-Salpeter equations for the mesons starting with the regularized action, 
see eqs. (3.25) and (3.27). In order not to disguise the used method by technical difficulties 
we will display it for the case of pions in the isospin limit only. For the other cases we will 
only cite the results. For details we refer the reader to Appendix A of ref. [62] d . 

The physical meson excitations are given by small amplitude fluctuations around the trans- 
lational invariant vacuum field configuration obtained by solving (3.28). Expanding the effec- 
tive action in Euclidean space 

A. — A-m ~\~ A-F 
A m = [d 4 x( ', t t i \ / t \ / ///„i.\/ • .W 1 ! • /I/', 

A F = -- -Trexp l-sp E 
2 Ji/a 2 s V 

/| A = // J • >U/VW • /'/.W 1 • /V\/ f .\/ • l'i.\IM j (4.24) 



furthermore, we will restrict ourselves to proper time regularization. 
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up to second order in the pion fluctuations of the pseudoscalar field Q(x) = 2tt(x) allows us 
to extract the inverse propagator for the pseudoscalar mesons. 6 Since we are not interested in 
the fluctuations of the scalar field we substitute it by the vacuum expectation value (VEV) 
(E) = mJT. The complex field M (3.16) is then given by 

M = mU } U = me l&(x) = me l2v(x) . (4.25) 

Obviously, the VEV of the pseudoscalar field Q(x) is zero. The chiral field U is expanded for 
small-amplitude fluctuations of the pseudoscalar meson field tt = 7t 8 (t 8 /2) as 

U = e 2 ^ = 1 + 2ztt - 2tt 2 + . . . . (4.26) 

The factor 2 is hereby introduced for later convenience. 

As the stationary point of the action occurs for tt = there is no linear term. First, we 
will consider the mass term of the mesonic action A m 

A m = ^- [ d 4 x ((m - m ) 2 + 4m°m tt 2 ) + 0(tt 3 ). (4.27) 

Introducing the Fourier transform of the fluctuating field, ff(q), the bilinear term of (4.27) is 
obtained to be: 

d 4 q 1 „• , , ; , , m°m 



r 9 uEo'I-^tt. (4-28) 

In order to expand the fermion determinant Af we rewrite the operator 

plpE = A + A 1 + A 2 + . . . (4.29) 

where Ak is of k-th order in the field tt. Using eqs. (4.24), (4.25) and (4.26) we obtain 

A = d 2 + m 2 
A x = 2 75 m(^7r) 

A 2 = 0. (4.30) 

The last relation is valid in the isospin limit only. Obviously, only derivatives of tt can occur 
as a consequence of the chiral invariance of Det^^fe. Using now 



1 ds ( \ 
A F = -- — Trexp (splpE 

2 Ji/A 2 s V 



2 Jl/A 2 s 

= _I I™ ^ Tre -^o + I r ds f 1 d (TTe-< Ao A 2 e- s ^ Ao 

2 Jl/A 2 S 2 Jl/A 2 Jo 

— / dss / d( d^Tre- 8 ^ A^- 8 ^-^ A ie - sCAo + 0{tt 3 ) (4.31) 
2 Jl/A 2 Jo Jo 

we can systematically expand this part of the action. The flavor (or isospin) trace only gives 
an overall factor 2, the color trace a factor N r and the Dirac trace a factor 4. The functional 



2 At this order the imaginary part the action does not contribute. 
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trace will be done using momentum eigenstates Tr(q). As A 2 = in the isospin limit we have 
to calculate the following parameter integral 



s((+v)(k 2 +m 2 ) e -s(l-(- V )((k+ q y +m 2 ) 
+ e -s«+v)(k 2 +m 2 ) e - 

dae- sa{k2+m2 h- s{1 - a ^ k+q ^ +m2 \ (4.32) 



-s((+ V )(k 2 +m 2 ) -s(l-(-n)((k+g) 2 +m 2 ) 



Since the momentum k is only appearing in this exponential we may shift k — > k — (1 — a)q 
in the second term without changing the value of the integral yielding 



^ ^ da e -.((*+(l-") 9 ) 2 + "(l-") 9 2 +- 2 ) = r da e -s(k 2 + a { l-a) q 2 + m 2 )_ (433) 



Using 



4iV c / -^-e^ 2 = ^ (4.34) 



2tt 4 4tt 



as well as the definition of the incomplete T-function, the term bilinear in pseudoscalar fields 
arising from Af reads 

d 4 o 3 1 • • N r 1 

^7 E ^'("^M?) ("« V'^ I daT (°, + "C 1 - «)9 2 ]/ A2 )- (4-35) 
This allows us now to extract the inverse pion propagator, 

D;\q*) = -^-Il( q % (4.36) 
where the polarization operator H(q 2 ) is given by 

f( q 2 ) = m 2 ^- f 1 dx r(0,[m 2 + x(l -x)q 2 ]/A 2 ). (4.37) 

47T 2 JO V 1 

The Bethe-Salpeter equation which determines the physical meson masses m v is equivalent 
to the condition that the meson propagator has a pole: 

D; 1 (q 2 = -m 2 T ) = 0. (4.38) 

Note that f 2 (q 2 = —m 2 ) then is the corresponding meson decay constant. We want to empha- 
size here that the Bethe-Salpeter equation (4.36) - (4.38) is the one in ladder approximation 
and that no other approximations as e.g. gradient (subsection 4.1) or heat kernel expansions 
have been made. This also implies that the decay constants are evaluated on the corresponding 
meson mass shell. The pion decay constant is then given by 

fl = m 2 ^- f 1 (forfo, [m 2 -x(l - x)ml]/A 2 ). (4.39) 
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Table 4.1: Mass parameters fixed in the meson sector of the NJL model. The kaon decay 
constant fx is predicted. 



m (MeV) 


m s (MeV) 


m° s /m° 


MMeV) 


350 


577 


23.5 


104.4 


400 


613 


22.8 


100.3 


450 


650 


22.5 


97.4 


500 


687 


22.3 


95.5 



One clearly sees that in the chiral limit (m^ = 0) the expression 

/ 2 = m 2 ^r(0,(m/A) 2 ) (4.40) 

calculated by means of a derivative expansion becomes exact. Using eqs. (4.36) and (4.37) it 
is trivial to show that in the chiral limit m° = the Bethe-Salpeter equation (4.38) is solved 
by setting q 2 = m 2 = 0, i.e. the pions are Goldstone bosons. The mass shell condition (4.38) 
can be expressed as 

mlfl = = 2m°(uu) — ^— - « 2m°(uu) (4.41) 

G\ m — m u 

where we used the gap equation (3.28) to express the coupling constant G\ in terms of the 
quark condensates. This approximate relation is the Gell-Mann-Renner-Oakes relation [41], 
see eq. (2.24). 

The generalization of the above calculation to the case of unequal quark masses, and 
especially the case of three flavors, can be found in Appendix A of ref. [62]. For the inverse 
propagator for the pseudoscalar mesons one obtains 

D7U* ) = - — 7^ " n* (? 2 )Wiy. (4.42) 



4d 

The polarization operator Hij(q 2 ) is given by 

IItf(? 2 ) = q 2 fUq 2 ) + ~ m^ffa 2 ) - ^(m? - m 2 )(^l _ Ml) (4 . 43) 

J ZL J \ ryi ■ ryi ■ / 



wherein 

fW) = 7( m « + m i) 2 T^ I' dx r ( ' K 1 " x ) m ^ + xm2 i + X ( X ~ X h 2 ]/ A2 ) ( 4 - 44 ) 
4 471"^ Jo v 7 

while the quark condensates are defined in eqs. (3.28). Obviously, these expressions coincide 
with the one given above if the flavor symmetric limit is chosen. For the kaon one simply 
uses i = s and j = u(d) } e.g. the kaon decay constant can be obtained from eq. (4.44) by this 
choice of indices and evaluation at q 2 = —m 2 K which in turn is obtained as the root of (4.42). 

Numerical results are displayed in table 4.1. As input serves m v = 135MeV, rriK = 495MeV 
and f v = 93MeV. The up constituent quark mass is chosen as independent variable. One 
observes that the NJL model in the proper time regularization scheme underestimates the 
experimental value fx = 113MeV. 
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For the other mesons (scalars, vectors and axialvectors) we will only give the propagators 
for the case of two flavors in the isospin symmetric limit. For the inverse propagator of the 
scalar mesons a and a one obtains 

D-'iq 2 ) = - (q 2 + 4m 2 )./V). (4-45) 

Note the similarity to the pion propagator (4.36). As a consequence the scalar meson mass 
can be written as 

f2/_ m 2\ 

ml = 4m 2 + ml J —^ ^V + mJ (4.46) 

P{-ml) 

where for the last relation we have assumed f 2 (q 2 ) (4.37) a slowly varying function. However, 
there is a problem. From (4.46) one sees that the scalar meson mass lies above the quark- 
antiquark threshold, i.e. m 2 > 4m 2 . In this case / 2 (— m 2 ), or equivalently -D" 1 , acquires an 
imaginary part-^ 



Imp; 1 ) = (q 2 + Am 2 ) — m 2 Jl + 4m 2 /? 2 for 4m 2 < -q 2 < 4(m 2 + A 2 ). (4.47) 

47T V 



Therefore the extraction of meson masses via the Bethe-Salpeter equation are obscured by 
the unphysical quark-antiquark threshold which is present in the (non-confining) NJL model. 

For the calculation of the vector meson propagator it is convenient to split the expressions 
in the expansion of the fermion determinant in longitudinal and transverse ones. A naive 
extraction of the longitudinal part of the inverse propagator yields 

m 2 r(-l, m 2 /A 2 ) - £ dx {{m 2 + x(l - x)q 2 )T(-l, Y) + 2(x - T(0, Y)j (4.48) 

where 

Y = (m 2 + x(l - x)q 2 )/A 2 . 

This expression identically vanishes as can be shown by a formal expansion in powers of q 2 
and evaluating the Feynman parameter integrals for each power of q 2 separately. Note that it 
is not sufficient to argue that the proper time regularization is gauge invariant (and therefore 
the longitudinal part has to vanish). The small amplitude expansion for vector mesons which 
keeps only terms up to second order in the vector (gauge boson) field but all derivatives does 
explicitly violate (local) gauge invariance. 

As the expression (4.48) vanishes only the mesonic action A m contributes to the longitu- 
dinal part of the inverse vector meson propagator. The latter is given by 

m-r r(q) = (V - <^ M U _ t^L±_ u 49) 



where 



g 2 (q 2 ) 47r 2 Jo 



N r 1 

dxx(l - x)T(0,(m 2 + x(l - x)q 2 )/A 2 ). (4.50) 



•^For — q 2 > 4(m 2 + A 2 ) additional (even more unphysical) imaginary parts appear. 
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The vector meson mass is determined by the on-shell condition 

m 2 v = ^T9 2 (/ = -m 2 v ) (4.51) 

and g v = g 2 (q 2 = —m v ) is the universal vector meson coupling constant taken at the on-shell 
mass. 

Including the axialvector meson changes the formulas for the pion propagator. This is due 
7r — a\ mixing, i.e. the occurrence of a term 

f 2 {q 2 )iq^{q)A^-q) 

in the small amplitude expansion. This modifies the inverse pion propagator: 

D-Hq 2 )-q 2 ^ — (4 52) 

Especially, the pion decay constant is then given by 

fl = l + 4G 2 P(-ml) (453) 

where the pion mass is determined by the root of eq. (4.52). As a further consequence of tt — a\ 
mixing the axialvector propagator is different from the vector one. The axial-vector mass, 

m\ = g 2 (q 2 = -m 2 A ) (±- + f 2 (q 2 = -m 2 A )^ , (4.54) 

is significantly larger than the vector mass and generally lies above the quark-antiquark thresh- 
old. 

This subsection may be summarized as follows: The propagators for the scalar and pseu- 
doscalar fields may generically be written as 

q l + m l \q l ) 

while those for the vector and axialvector fields read 

D^(q 2 ) = — ( 8 ^-J^-\ (4.56) 

with the functions Z(q 2 ) and m 2 (q 2 ) given in table 4.2. For the scalar and axialvector mesons 
their masses lie above the quark-antiquark threshold making their determination via the 
Bethe-Salpeter equation doubtful. It should be remarked that a procedure, which extrapolates 
the polarization tensors from below to above the quark-antiquark threshold, has been proposed 
in ref.f 
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Meson 




m 2 (q 2 ) 


cr,a 


/V) 


Am 2 + m°m/(G 1 f 2 (q 2 )) 


7T 


f 2 (q 2 )l(l + AG 2 f 2 (q 2 )) 


m°ml{G 1 f 2 {q 2 )) 




1 la 2 (a 2 \ 






i/g 2 (q 2 ) 


g 2 (q 2 )(l/±G 2 + f 2 (q 2 )) 



Table 4.2: The functions Z(q 2 ) and m 2 (q 2 ) which determine the meson propagators, see 
eqs. (4.55) and (4.56). The quantities f 2 (q 2 ) and g 2 (q 2 ) are given in eqs. (4.37) and (4.50), 
respectively. 
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Chiral anomaly 



In this section we will discuss how the chiral anomaly (see section 2.3) is represented in the 
effective meson theory (3.25). Defining the quantum theory via path integrals the anomalous 
symmetry breaking is realized by the non-invariance of the functional integral measure [49]; 
despite the fact that the classical action appearing as 'weight' factor under the path integral 
is invariant under the considered transformation. The integral measure 



VqVq 



(4.57) 



is not invariant under chiral rotations 



q(x) -> x(x) = (U(x)y 5 q(x) = e l&(x) ^q{x) 
q(x) -> x{x) = q(x) (U(x)y 5 
but rather acquires a phase J(0) ^ I 

V X V X = J(&)VqVq . 



e = e a 



(4.58) 



(4.59) 



This implies that for the effective meson theory (3.25) the anomaly is contained in the 
fermion determinant Af- The proof of anomalous chiral symmetry breaking, reflected by the 
relation J(0) 7^ I, is given in Appendix A for the simplified case of a Dirac operator containing 
besides the usual kinetic and mass term only a vector field 



In this case one obtains 



J(0) = exp 



.Nc_ 

; 8tt 2 



m 



m 



d 4 xQ(x)tr F (FF) 



(4.60) 



(4.61) 



where is the field strength tensor corresponding to the vector field and F^ is the 
dual tensor. Note that the Jacobian J(0) is a pure phase factor, |J(0)| = 1, and therefore 



35 



contributes to the imaginary part of the effective action. In the presence of a vector field this 
factor will in general differ from one. 

The expression (4.61) is the Jacobian for infinitesimal Abelian chiral rotations. Addition- 
ally, we want to know the Jacobian J(0) for finite non-abelian chiral transformations. It can 
be calculated by functional integration of the differential anomaly 

J(0) = j V® Sj j§^ = exp(z^lwzw). (4.62) 

The term Awzw represents the logarithm of the integrated anomaly and is the Wess-Zumino- 
Witten (WZW) term [52]. In general it is a very complicated functional. We will therefore 
discuss in more detail only the special case of a flavor singlet vector meson (cu-meson) which 
couples to the baryon number and a SU(2) chiral field, 



V^X) = LO^X^lp 



U(x) = e l&{x \ Q(x) = <d a {x) T —. (4.63) 
The WZW term is then given by 

Avzw = N c j d 4 x lo^B^x) (4.64) 

where 

B^x) = -^-e^\v{L v L K L x ), L v = IjhjJJ (4.65) 

is the winding number current related to the chiral field U(x). Due to its topological nature 
the current B^ is conserved independent of the explicit form of the chiral field Q(x). 

In order to reveal the physical nature of this current we consider the quark generating 
functional treating the vector field as an external source 

exp(*4p) = j 'Dq'Dqexp (^j d 4 xq(i((J) + Y) — rh°)q^j 
= J{e)exp(A[V}) 

= exp ^iN c j d 4 xV^x)B^x) + A[V}^ (4.66) 

where A[V] contains the non-anomalous terms of the effective action. From the first equation 
we obtain 

( ^ ={q{x)* 1 »q{x))=:j%{x) (4.67) 



iSV^x 



MW/ v=o 

i.e. the baryon current of the quarks. The last equation allows us to relate the baryon current 
to the topological current 

f B (x) = B»(x) + . . . (4.68) 

where the dots indicate the contributions from the "normal" terms «4.[V]. This leads to the 
following interpretation: The baryon current of the original fermion (quark) theory is carried 
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by the topological current of the chiral field in the bosonized effective meson theory, at least 
in leading order gradient expansion. As mentioned above this topological current is conserved 
independent of the specific form of the chiral field whereas the baryon current of the underlying 
quark theory is the (conserved) Noether current of flavor singlet phase transformations. So, 
the dynamical property of the quark theory has turned into a purely topological property in 
the effective theory. 

It should be remarked that continuing to Euclidean space is equivalent to using Feynman 
boundary conditions in Minkowski space. This implies that the generating functional corre- 
sponds to the "vacuum to vacuum" transition amplitude. For infinitely large Euclidean times 
jg thus represents the vacuum part of the quark baryon current only. 

This immediately leads to the question how the pseudoscalar meson field Q(x) which is 
bosonic and spinless can describe properties of baryons which are fermions with spin s = 1/2. 
We will see in section 5.2 that a sufficiently strong topological non-trivial chiral field with 
winding number n polarizes the vacuum or Dirac sea of the quarks so strongly that n of 
the valence quark levels are bound tightly enough that they join the negative Dirac sea: 
their energy is negative. As the physical vacuum is defined as the state with lowest energy 
all negative energy levels are occupied in the vacuum. The physical vacuum carries a non- 
vanishing baryon number if quark states are bound in the Dirac sea. So it is not really the 
chiral field itself which carries the baryon number but rather the polarized vacuum. As we 
cannot observe the vacuum but only the polarizing meson field we relate the baryon number 
to the chiral field. In that sense meson fields carry baryonic charge. This fact is the underlying 
feature for the description of the baryons as chiral solitons as will be discussed in more detail 
in the following chapter. 

Let us close this section with a comment. Using the invariance of the classical action and 
the Jacobian (4.61) the anomalous Ward identity 



can straightforwardly be derived (see Appendix A). Note that even for m° = the axial 
singlet current j£ is not conserved if vector fields are present despite the fact that the classical 
Lagrangian is invariant under chiral rotations for m° = 0. The axial anomaly as formulated 
by the anomalous Ward identity (4.69) is known as Adler-Bell-Jackiw anomaly [36]. Usually 
it is required that in fundamental (gauge) theories anomalies should be not present or should 
be canceled by other effects. As we are considering effective low-energy models anomalies can 
be present. They even have measurable consequences as e.g. the decay 7r° — > 27. 

Historically, the anomaly was discovered in perturbation theory. The one-loop diagram 
which is responsible for the decay 7r° — > 27 is shown in figure 4.1. The path integral derivation 
of the anomaly is a non-perturbative one and demonstrates that higher order terms do not 
contribute to the anomaly. If one calculates the Adler-Bell-Jackiw anomaly in a renormaliz- 
able theory (as QCD) in which the cut-off A is taken to infinity at the end of the calculation 
only the triangle diagram figure 4.1 contributes. However, if one works within an effective 
non-renormalizable theory where the ratio ra/A has to be kept finite, not only the triangle 
but all higher order diagrams contribute. Remember that the limit A — > 00, i.e. r — > 0, (see 
Appendix A) was necessary in order to make the contributions from h n} n > 3 vanish. In 
the literature [64] there are claims that the NJL model is not appropriate for the treatment 
of the anomaly and the 7r° decay because the triangle diagram gives only about two thirds 
of the experimentally determined decay amplitude using a finite cutoff. However, together 
with the finite cut-off also higher order terms should be included and presumably yield the 




(4.69) 
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Figure 4.1: The triangle diagram which is responsible for the anomaly. 

missing third of the decay amplitude [65]. Note also that not regularizing the imaginary part 
of the action leads to a suppression of higher order terms for the 7r° decay, i.e. in this case the 
triangle diagram gives already the correct decay amplitude. 
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5 Chiral Solitons 



This chapter is devoted to a brief discussion on the interpretation of solitons as baryons. 
Also an outlook is provided on calculations which are performed at length in chapters 6 and 
7. For these introductory remarks on soliton physics we assume that some effective chirally 
invariant meson theory is given which allows for stable soliton solutions. When studying low 
energy properties it is suggestive to take into account only those mesons with the largest 
Compton wave-length, i.e. with the smallest mass. As explained in the preceding chapters 
these are the pseudoscalar would-be Goldstone bosons tt(x) of spontaneous chiral symmetry 
breaking. Again the non-linear realization will be adopted 

U(x) = exp(iQ(xj) (5.1) 

which defines the chiral field U(x). 
5.1 Topological properties 



For two flavors a static chiral field maps the coordinate space into the group of isospin 

U :R 3 -> 577(2). (5.2) 

Solitons are finite energy field solutions to the Euler-Lagrange equations. Finiteness of the 
energy requires the energy density to vanish for r = \r\ — > oo. This implies that the soliton 
configuration U(x) asymptotically approaches a constant value (independent of the orientation 
r) which, as a consequence of chiral symmetry, can be chosen to be unity: 

U(x) — > 1 for r — > oo. (5-3) 

Hence all points at spatial infinity are identified thereby compactifying 1R 3 to a three dimen- 
sional sphere S 3 . Since the group manifold of SU (2) is also S 3 a static chiral field configuration 
with the boundary condition (5.3) represents a mapping from S 3 to S 3 

U :S 3 -> S 3 . (5.4) 

These mappings are distinguished by the winding number v £ X which is a topologically 
invariant functional of the chiral field 

u[U] = J d 3 xB°(x). (5.5) 

Here B^(x) is the topological current (4.65). 

A non-trivial mapping Uq(t) is given by the so-called hedgehog ansatz (cf. subsection 
6.2.1) 

U (r) = exp (it ■ r 0(r)) . (5.6) 

For this field configuration the isospin vector points into the radial direction. Uniqueness of 
the chiral field at r — > requires the chiral angle to obey the boundary condition 

0(r = O) = -n7r, neZ. (5.7) 
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For the hedgehog configuration (5.6) the spatial components of the topological current vanish 
while 



o / \ 1 s sin 2 0(r) 

B » = 7T- 2 &(r 5.8 

Ztt r z 

Without loss of generality one may choose 0(r — > oo) = to satisfy (5.3) yielding 

v[U] = n. (5.9) 

Let us emphasize at this point that the winding number v = n is a purely topological property 
of the chiral field, which a priori is not related to any physically relevant quantity. 

5.2 Emergence of the soliton 

In general the effective meson theory results from integrating out the quark and gluon 
degrees of freedom contained in QCD. Then the baryons, which are originally built from Nc 
quarks, emerge as solitons of meson degrees of freedom. In order to explain how the baryonic 
character is encoded in the topological properties of the soliton we examine the spectrum of 
constituent quarks in the background of the chiral field. The corresponding Dirac Hamiltonian 
reads 

h = a ■ p + (3m(Uy 5 , (5.10) 

where m denotes the mass of the constituent quarks. The hedgehog field configuration (5.6) 
violates the spin and isospin symmetries' 1 , i.e. 



but preserves the grand spin 



[h,j] = -[h,t]^0, (5.11) 



G = j + t (5.12) 



symmetry, i.e. [h,G] = 0. Hence the eigenfunctions of h carry good grand spin and parity. 
In order to examine the spectrum of (5.10) we parametrize the chiral angle by [66] 

^ . . f 1 — 7t- for r < a , 

e(r) = -"U • ( 1 

This special form is motivated by the fact that for small r the chiral angle is linear whereas 
in the chiral limit it is proportional to 1/r 2 for large r. The chiral angle and its derivative are 
continuous at the matching point r = a, which parametrizes the spatial extension of the chiral 
angle. Figure 5.1 shows the eigenvalues e v of the Dirac-Hamiltonian (5.10) for a chiral field 
with winding number n = 1 as a function of the strength of the chiral field measured by a ■ m. 
As the strength of the chiral field increases, the lowest valence quark state in the G v = + 
channel becomes strongly bound and eventually joins the Dirac sea. Figure 5.2 shows the 
spectrum of constituent quarks in the background field of a hedgehog with winding number 
n = 2. In this case a second valence quark state G v = 0~ becomes bound in the Dirac sea for 
sufficiently large a • m. In general, strong chiral fields with winding number n bind exactly 



"The quantities j and t denote the single quark total angular momentum and isospin operators, respectively. 
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Figure 5.1: The eigenvalues of the Dirac Hamiltonian (5.10) in the background of a chiral 
held (5.13). Displayed are the lowest eigenvalues in the + (full lines) and 0~ (dashed lines) 
channels. 




Figure 5.2: Same as figure 5.1 for n = 2. 
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n valence quark orbits in the negative Dirac sea. Obviously the baryon charge of the valence 
quarks gets encoded in the topological structure of the soliton when describing baryons in the 
framework of effective meson theories. 

By definition, all positive energy states are empty in the vacuum configuration while all 
negative energy states are occupied. Hence for sufficiently strong chiral fields the vacuum 
possesses a baryonic charge which is identical to the winding number. As the polarizing 
meson field is observable rather than the vacuum, the baryonic properties of the charged 
vacuum are attributed to this meson cloud. This justifies the statement that baryons emerge 
as solitons in the effective meson theory. Strictly speaking, however, the baryonic charge is 
carried by the polarized Dirac vacuum. To make this scenario more explicit it is helpful to 
consider the baryon number carried by the vacuum, which is defined as the asymmetry of the 
Dirac spectrum (see eq (6.12)) 

IT ac = -Vsgn( ei ,). (5.14) 

For a vanishing chiral field the Dirac spectrum is symmetric and therefore B vac = 0. However, 
as the strength of chiral field increases the energy eigenvalue of the valence quark orbits 
eventually reverse their signs leading to B vac = n. 

These studies provide the following reasoning for the solitonic picture of baryons: The 
valence quarks of the baryons with baryon number B generate a strongly localized chiral field 
of winding number n = B. In the background of this field the valence quarks are bound in 
such a way that the baryonic charge is carried by the polarized vacuum, i.e. B vac = n. Once 
the valence quarks have joined the Dirac sea, the polarizing chiral meson cloud is interpreted 
as the baryon. 

5.3 Semiclassical quantization 

As already mentioned after eq (5.10) neither spin nor isospin represent good quantum 
numbers of the soliton. Here we will briefly describe a treatment which generates states 
carrying quantum numbers of physical baryons. For details we refer to section 7.1. 

Time-independent rotations in coordinate- and/or isospace do not alter the energy of 
the soliton. In order to perform the semiclassical quantization one requires time dependent 
solutions to the equations of motion. These are approximated by allowing these rotations to 
adiabatically vary in time. Due to the grand spin symmetry the rotations in coordinate- and 
iso-space are equivalent. It therefore suffices to write 

U(r,t) = R(t)U (r)ltf(t), (5.15) 
where R(t) denotes an SU(2) isospin matrix. Defining the angular velocity Q by 

rf(t)-R(t)= % -n-T, (5.16) 
the energy functional E[U] of the rotating soliton configuration can be expanded in powers of 

n h 

E[U = RU Rt] = E[U ] + ^n t Q ll [U ]n 3 + ■■■. (5.17) 

'Since isospin breaking is neglected the rotation matrix R does not appear explicitly. 
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The term linear in J? is missing due to time reflection symmetry which holds for the isospin 
symmetric two flavor model. The moment of inertia tensor 



[U ° l ~ dnm 3 



= S tJ a 2 (5.18) 

f2=o 



is diagonal because no direction in isospace is distinguished. The canonical quantization of 
the coordinates R corresponds to the replacement 

a 2 f2^J, (5.19) 

with J being the collective spin operator. The quantized energies of the baryons are finally 
obtained to be 

E[U] = E[U ] + (5 .20) 

where J( J + 1) is the eigenvalue of J 2 . 

Due to the grand spin symmetry of the hedgehog the spin of the quantized soliton coincides 
with the collective isospin 7 up to a rotation. The quantized hedgehog soliton therefore yields 
a tower (rotational band) of states with I = J. Generalizing this treatment to three flavor 
models shows that for an odd (even) number of colors the spin is constrained to half-integer 
(integer) values [52, 26, 27]. For N c = 3 the nucleon (7 = J = |) and the A(7 = J = |) 
resonance thus represent ground and first excited states, respectively. 

Both, the classical soliton mass and the moment of inertia are of order Nc- Therefore the 
classical soliton mass 7?[77 ] and the rotational energy 7(7 + l)/2a 2 are of order Nc and l/iVc, 
respectively. However, a systematic expansion of the baryon energy in powers of 1/Nc 

E[U] = E ( ~^ + E^ +£W + ••• (5.21) 

contains in addition to (5.20) a term of order (l/iVc) , E^°\ This term, which presumably 
dominates over the rotational energy, is generated from small amplitude mesonic fluctuations 
off the soliton and is associated with the quantum corrections to E [Uo\. Naively one would 
expect 

^ (0) = ^(E^-E^ 0) )> ( 5 - 22 ) 

\ i i J 

where u t and uf ] refer to the eigen-frequencies of the meson fluctuations in presence and 
absence of the soliton, respectively. However, the expression (5.22) is UV divergent and 
thus subject to regularization. Nevertheless, (5.22) is illuminating since it indicates that the 
dominant contribution to E^ stems from the zero-modes (uii = 0) because there are no 
counterparts in the absence of the soliton. Moreover, the contributions of the zero-modes is 
negative and hence causes a substantial reduction of the classical mass. This is a desired feature 
since the classical mass commonly overestimates the baryon masses once the parameters are 
fitted to mesonic data (see chapter 4). Recently, it has been shown for the Skyrme model 
that a renormalized version of (5.22) indeed predicts baryon masses which well agree with 
experimental data [24, 25]. The analogous calculation [67] for the NJL model is presented in 
section 7.4. 
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6 Static solitons of the Nambu— Jona-Lasinio model 



In the present section we will discuss the emergence as well as properties of static solitonic 
meson configurations within the NJL model. The most important ingredient of the effective 
meson action of the bosonized NJL model is the fermion determinant (3.25). In practice we 
have to evaluate this determinant in the presence of non-perturbative meson field configura- 
tions. This determinant is given in terms of the eigenvalues of the Dirac operator p. However, 
since det(z/3) = I one may equally well consider 



i/3p= id t - h 



(6.1) 



which introduces the one-particle Hamilton operator h. For the time being we consider h 
to be static and Hermitian. The discussion of more general cases will be postponed to later 
subsections. For static fields the eigenvalues of (6.1) separate into the eigenvalues of id t 
and h. The fermion fields assume anti-periodic boundary conditions on the time interval T. 
Therefore the eigenvalues of id t are given by the Matsubara frequencies fi n = (2n + 1)tv/T with 
n = 0, ±1, ±2, . . .. Denoting furthermore the eigenvalues of h by e v the fermion determinant 
is obtained as the product" [68, 70] 



Det(# = nn 



2n + 1 
T 



-IT 



v n>0 



(2n + 1W 



(6.2) 



since for static configurations id t and h may be diagonalized simultaneously. The constant 



c = UU 



v n>0 



2n + 1 
T 



-IT 



(6.3) 



does not depend on the dynamical properties of the system and may hence be absorbed into 
the integration measure. The product over n is readily carried out 



Bet(p) = C~[[cos 



C exp 



Y[ (1 + exp [-iT\e p 



(6.4) 



With the introduction of occupation numbers r\ v = 0, 1 the product over v may finally be 
expressed as 



]J (1 + exp [— zT|e^|]) = ^ exp 



(6.5) 



where the sum goes over all possible combinations oirj v = 0,1. Then the fermion determinant 
acquires the form 



Det (p) = C exp [iAo] ^ exp \iA v 



(6.6) 



"See also chapter 9 of ref. [69]. 
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which provides a natural decomposition into vacuum 

A = 1 -TJ2K\ (6.7) 

and valence (anti-) quark 

Ap» } = -Tj2v^\ = -TEp» } (6.8) 

contributions to the fermion determinant. 

In order to equip the occupation numbers r\ v with a physical meaning let us consider the 
baryon number current which is defined as the average 

r(x) = (q(x)Yq(x)) = ^y lo S Det (P- if) L w=0 - (6-9) 
Treating as a perturbation in the eigenvalue problem (h + ftY) 1 ^ = e ^Vv reveals that 

Se " v , {x)=0 = ^( x )^M x ) (6-10) 



with Vv(^) and e v being the eigenstates and -values of h. It is then easy to see that according 
to (6.6) the baryon number current is additive in vacuum and valence parts [70] 

j"( x ) = jS(x)+ft(x) 

joi x ) = -^X^sM^Vv^b^Vv^) 

Jv( x ) = J2vv s g n ( e v)i ) A x h >J '4>A x )- (6-ii) 

Taking into account that the eigenfunctions %\) v are properly normalized one obtains the baryon 
number 

B = E (V " \) sgn(e„). (6.12) 

Thus, considering a special set {rj^} of occupation numbers confines the system to a sector 
with definite baryon number. 

Up to now we have ignored the fact that the vacuum part of the fermion determinant 
(6.7) is divergent and hence needs regularization. We will address this problem in the next 
subsection when extracting the energy functional from Aq. Furthermore, a more rigorous 
treatment of symmetry currents for the regularized theory will be presented in chapter 7. 



6.1 The energy functional 



We have already observed that the contribution to the energy functional due to the explicit 
occupation of the valence quark orbits is given by (cf. eq (6.8)) 

E v = ^2lA e A (6.13) 
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with e v being the eigenvalues of a one-particle Hermitian Dirac Hamiltonian h. This quantity 
may be extracted from the Dirac operator via eq (6.1). 

The vacuum (or ground state) contribution E to the energy is extracted from the func- 
tional integral by continuing to Euclidean times r = ix and observing that for large Euclidean 
time intervals, T — > oo, the ground state provides the dominant contribution, i.e. 

lim Det(z/3|fe) oc exp (-E T) . (6.14) 

T^oo 

Here pE denotes the Euclidean Dirac operator which is obtained from p (6.1) by analytic 
continuation. It is important to note that in Euclidean space r has to be considered a real 
quantity. 

An Hermitian one-particle Dirac Hamiltonian has been the starting-point of the preceding 
considerations. We will waive this assumption from now on and allow h to contain anti- 
Hermitian parts as well. Nevertheless, the Euclidean Dirac operator for static meson fields is 
decomposed into a temporal part and a static Euclidean Hamiltonian 

if3pE = -d T - h. (6.15) 

Commonly the non-Hermiticity of h stems from continuing the time components of (axial) vec- 
tor fields to Euclidean space: Vo — > —iV 4 and A — > —iA 4 . The Euclidean Dirac Hamiltonian 
then reads 

h = a ■ p + V 4 + 75 A 4 + ia. ■ V + i l5 a ■ A + f3(P R M + P L Aft). (6.16) 

Furthermore our notation for the (axial) vector meson implies to take and anti- 
Hermitian. Especially the anti-Hermiticity of V 4 and A 4 causes the eigenvalues of (6.16) 
to be complex e v = + ze^. Hence the eigenvalues \ n ^ v of the operator d T + h are given by 

K, v = -i^n + t v = -ifL n + + iel. (6-17) 

Let us (for the moment) ignore regularization and perform manipulations analogous to those 
described above for the special case of h being Hermitian. The fermion determinant in Eu- 
clidean space can again be shown to separate into vacuum and valence contributions 

Det(|fe) =C exp(-TE ) ^ exp (-T$>*<U (6.18) 

with the vacuum energy E = —(Nc/2) e v . The quantities e v are defined in terms of the 
real and imaginary parts of the eigenvalues of h 

e„ = |^| + ;sgn(e?)e£. (6.19) 

The contribution of the valence orbits to the (Euclidean) energy functional can be read off 
(6.18) as 

E V = E* + xE'y = NcYl VA^I + iNc £ ^sgn(e?)e£ (6.20) 

in the case that h contains anti-Hermitian parts. We have also made explicit the dependence 
on Nc- 
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At this point a few remarks on the analytic properties of the eigenvalues e v are in order. 
Considering the generalization of the Euclidean Dirac Hamiltonian (6.16) defined by substi- 
tuting Vo — > zV± and A — > zA 4 the eigenvalues e v actually become functions of the complex 
variable z [71, 72]. It has been observed that for physically motivated field configurations these 
eigenvalues are analytic in z. This has been achieved by numerically verifying that the Laurent 
series for t v {z) indeed reduce to Taylor series [71]. The relevant Cauchy integrals are computed 
by parametrizing the path ((f) = 8e lLp — z . Here z refers to the center of the Laurent expan- 
sion while a variation of 8 provides a means to extract the radius of convergence. Then the 
eigenvalues of the generalized Dirac Hamiltonian are computed along the path (0 < <p < 2tt) 
and substituted into the Cauchy integrals. The eigenvalues are then found to exhibit an ana- 
lytic structure for a radius of convergence of the order unity or even larger. This analyticity 
is not a priori clear 6 because the eigenvalues represent roots of the characteristic polynomial, 
which carries a large degree (infinitely large in the continuum case). For the relevant field con- 
figuration it has furthermore been shown that e v (z)* = e v (z*). Obviously all quantities which 
depend on e^(z,z*) = {t v {z) + e v (z*)) /2 and el(z,z*) = —i(e v (z) — e v (z*)) /2 separately are 
no longer analytic functions of z. Such quantities are e.g. e v . Next we consider the special 
parametrization z = e* v , which for < <p < tt/2 describes the path connecting Euclidean 
and Minkowski spaces. In case the sign of e^(z, z*) is reversed along this path the analytic 
structure of E is obviously destroyed. In ref. [71] it has, however, been demonstrated that the 
total energy for a subsystem with unit baryon number, B = 1, can be written as 



The valence quark level (val) refers to the state with the smallest module Thus a change 
of sgn (e^i) does actually not destroy the analytic properties. As the other states (v ^ val) 
vary only mildly along the path z = e* v , analyticity is formally maintained for the unregular- 
ized energy functional when constraining oneself to configurations with unit baryon number. 
Unfortunately, a different type of non-analyticity exists. Level crossings along the path con- 
necting Euclidean and Minkowski spaces will appear if the time components V4 and/or A 4 
are strong enough. This makes a definition of a Minkowski energy functional by analytic 
continuation impossible. Thus, in order to attach physical significance to a given field config- 
uration one has to check that no such level crossing occurs. When these level crossings are 
avoided the unregularized energy functional is analytic and the continuation forth and back 
from Euclidean to Minkowski spaces can be performed. In this context it is then evident that 
models like the chiral quark model with the cu-meson included [74] indeed exhibit analytical 
structures. 

Next we have to face the problem of regularization which has been ignored for the static en- 
ergy functional up to here. Eq (6.18) demonstrates that for T — > 00 only the vacuum (ground) 
state contributes to the functional integral since only the term with all 77^ = survives in the 
sum (6.18). However, the expression for E is divergent and thus needs regularization. As in 
the study of the meson sector (c/. chapter 4) we will employ the proper-time regularization 
scheme [48] which substitutes the logarithm by a parameter integral (3.27). Unfortunately, 
this procedure is only applicable when the argument of the logarithm is positive. We therefore 
decompose the contribution of the fermion determinant to the mesonic action into real (Ar) 

'See e.g. appendix C of ref. [73]. 
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and imaginary (Ai) parts 



A f = Ar + A i = ^Tr log (p ] e Pe) + ^Tr log 



(6.22) 



In terms of the eigenvalues \ n ^ v (6.17) this decomposition is expressed as 
1 



and 



Ai 



(6.23) 



Note that at this point the rules for manipulating the logarithm have been used in the sense 
that Tr log(^fe) and Tr log(^) have been computed independently. This obviously represents 

an approximation because the traces of h and are associated with different Hilbert spaces. 
In ref. [75] it has been demonstrated that such a treatment may induce an error which is of 
quadratic or higher order in the time component of the (axial) vector fields. 

The proper-time prescription can obviously be applied to this form of Ar yielding 



A 



R 



E lo s ((« 



4? + (a 2 



E 



ds 

1/A 2 S 



exp < —s ( (0 



J\2 



+ (a 2 )}- 



(6.24) 
(6.25) 



According to the above discussion the expression (6.24) can only be considered as an approx- 
imation to the action in the presence of V4 and/or A 4 rather than being exact. 

For large Euclidean time intervals (T — > 00) the sum over n in (6.25) may now be replaced 
by an integral 



A 



T 



R 



E 



dz 



ds 



2lT Jl/A 2 s 



exp 



{s (* 2 + (e?) 2 )} 



(6.26) 



where we have furthermore shifted the integration variable z — — > z. As mentioned above, 
the limit T — > 00 allows one to extract the vacuum contribution to the real part of the energy 
functional from Ar — > —TEq which coincides with the corresponding expression derived from 
a Hermitian Hamiltonian [70]: 



E, 



R 



4^ 



N, 



2 ^ 



A 




exp 



(6.27) 



where 




erfc 



r R 



A 



(6.28) 



are the "vacuum occupation numbers" in the proper time regularization scheme, which for 
A — > 00 reduce to M v = 1. 



c Since the integral (6.26) converges absolutely, the sum over v and the integral over z may be exchanged. 
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For the imaginary part (6.23) we obtain (again to be considered with some caution as eq 
(6.24)) 

^ = \{H £ log(A„„) -EE log(A:,j) =\llt lo g|^ (6-29) 

\ u n= — oo u n= — oo / u n= — oo li£, n tj, 

where we have reversed the sign in the first sum over the integer variable n. Next we express 
Ai in terms of a parameter integral 

L v n=-oo J ~ 1 t v lAt u 

In analogy to (6.27) we may carry out the temporal trace in the limit T — > oo: 

i ^ Z" 1 „ _ 1'°° dz e 1 



Ai = —Y dXTV —- r^z = . (6.31) 

2 ^J-i J -co 2vr i(z - Xel) - K ' 

Care has to be taken when performing the ^-integration because only its principle value (V) 
is properly defined. Next the shift in the integration variable z — Ae^ — > z is performed. This, 
of course, also effects the boundaries 



Ai = —Y dXT lim / — — (632) 

2 z ?J-i M^oc,J-M-\el2Triz-€^ v ; 

Taking into account that A covers a symmetric range of integration, which allows one to reverse 
the sign in the A integral without changing the associated boundaries, the shift appearing in 
the ^-integral boundaries can be shown not to contribute as Ai — > oo. This yields 

A! = —Y I 1 d\T lim l M (6.33) 
2 V-'-i M ^°° J-m 2tt iz - K ' 

Now the integral over the parameter A may be done. Due to the principle value prescription for 
the z integration the terms of odd powers in z cancel. This results in a convergent expression 

It should be mentioned that the divergence, which disappears in the principle value prescrip- 
tion, has the drastic consequence that the energy functional is not invariant when shifting the 
to field by a constant u c . In that case the imaginary part ot the Euclidean action is augmented 
by TNc<-o c B where B denotes the baryon number of the original configuration. This result is 
very gratifying, as will be seen below. 

Although Ai is finite in the principle value formulation, the proper time regularization 
may be imposed by expressing the integrand as a parameter integral [76] 

^ Q> ds exp { _s ^ + (ef)2 ^ (6 - 35) 

which does obviously not diverge as A — > oo. Completing the evaluation of Ai in analogy 
to eqs. (6.25-6.27) we find for the contribution of the Dirac sea to the imaginary part of the 

Euclidean energy E* from Ai ^ - iTEl 

I —N c v-^ j , f 1, Ai not regularized , , 
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The upper case in eq (6.36) corresponds to the limit A — > oo. Obviously only the real part of 
the one-particle energy eigenvalue is relevant for the regularization of Ai. Eq (6.36) represents 
a regularization for Ai that only involves quantities which are strictly positive definite. 
The Euclidean "energy" functional 

E R + 1E 1 = E R + E R + i (E T + E v ) (6.37) 

exhibits the interesting feature that, as a quark state becomes part of the vacuum, both real 
(E R ) and imaginary (E 1 ) are continuous under the condition that the baryon number (6.12) 
remains unchanged. For a state \v) to become part of the vacuum means that e R changes 
sign. Assume that for the configuration under consideration one special orbit, say |val), has 
??vai = 1 as long as e R al > 0. In order to stay in the same baryon number sector r] va \ has 
to vanish as e R al reverses its sign (cf. eq (6.12)). Thus we demand f] va \ = (l + sgn(e^ ll )^ /2. 
Noting that erfc(0) = 1 one easily verifies that for e va i the terms which depend on the 
sign of e R al cancel in the sum (6.37). From this discussion it is also obvious that the valence 
quark orbit occupation numbers r\ v are dynamical quantities which functionally depend on the 
soliton configuration. 

The determination of the Euclidean energy functional is now completed and we have to 
obtain a Minkowski energy functional from this. In the discussion proceeding eq (6.20) we have 
already remarked that the one-particle energy eigenvalues t v {z) can be considered analytic in 
the complex variable z. Unfortunately, such a statement was shown [71] not to hold for the 
energy functional 

E(z, z*) = E R (z, z*) + iE\z, z*). (6.38) 

Here the dependence on z and its complex conjugate z* is due to the implicit dependencies 
e R,I (z,z*). The appearance of z* is caused by the fact that the regularization treats real and 
imaginary parts of the action in the complex plane separately. Remember that (formally) the 
unregularized action has been observed to be analytic in z. The functional (6.38) has been 
investigated with respect to its analytical properties [71]. It has been found that the Laurent 
series centered at the Euclidean point (z = i) has vanishing radius of convergence and the 
coefficients of the singular terms are non-vanishing. Stated more drastically: The analytic 
continuation of (6.37) does not exist. Hence the Minkowski energy functional for static, non- 
perturbative field configurations involving time components of (axial) vector fields cannot be 
obtained by means of analytical continuation. 

It has then been argued [71] that the generalized energy functional (6.38) can well be 
approximated by E R (i, —i) + zE 1 ^, —i) at the Euclidean point (z = i). This approximation 
gains further support by the fact that the regularization as discussed above is not without 
ambiguities at the quadratic order in the time components of the (axial) vector fields (see 
section 6.2.4). Adopting this approximation the analytic continuation to the Minkowski point 
(z = 1) is trivial, yielding 

E = E R + Eq + E R + Ey + E m . (6.39) 

The mesonic contribution E m is straightforwardly obtained by substituting the soliton config- 
uration into the purely mesonic part of the action A m - One might wonder whether there is a 
sign ambiguity for Eq v in the definition (6.39). This is not the case for the energy functional 
of a self-consistent soliton configuration because the eigenvalues of h are analytical and hence 
such a sign ambiguity may be absorbed in the definition of the time components V4 and A4. 
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Let us next mention two consistency conditions which have to be satisfied by the Minkowski 
energy functional. As already noted after eq (6.34) shifting the u field by a constant amount 
(u} c ) leads to a change of the (unregularized and thus analytical) energy functional of iNcBu c . 
We demand this to hold for the regularized energy functional in Minkowski space as well. As 
a matter of fact this condition is a consequence of global gauge invariance [i.e. a constant 
Vo behaves like a chemical potential). Furthermore, the current field identity (3.21) for the 
baryon current imposes a normalization on the profile function of the isoscalar u meson. This 
normalization is obtained by integrating the stationary condition for this field over the whole 
coordinate space. As this stationary condition is obtained from extremizing the Minkowski 
space energy functional we have available a second consistency condition. It is important to 
note that the energy functional (6.39) satisfies these two consistency conditions; at least when 
the imaginary part remains unregularized. Thus (6.39) provides a well suited object for the 
investigation of soliton solutions [76, 73, 75]. 

At this point it should be mentioned that two other approaches to include time components 
of vector mesons are discussed in the literature. In ref. [77] these components are not continued 
according to the discussion preceding (6.16). In that case the energy functional in the whole 
z-plane depends on the eigenvalue t v {z) only and is thus analytical. However, due to the 
regularization the global gauge invariance as well as the current field identity for the baryon 
number current are lost in Minkowski space. This approach has therefore to be abandoned 
because of physical reasons. In the second treatment [78, 79] the time components have indeed 
been continued. These authors have expressed the Euclidean energy functional (6.37) in terms 
of t v {z) and its complex conjugate e*(z) = e v (z*). Then the energy functional depends on z and 
z*\ E = E(z } z*). It should be evident from the preceding discussions that the dependence on 
e*(z) and thus z* purely originates from regularization and treating real and imaginary parts 
of the action differently. As in Euclidean space z = i = —z* the authors of ref. [78, 79] have 
then set E = E(z } —z) and claimed that the continuation of E to z = 1 yielded the Minkowski 
energy functional. It should be noted that this approach preserves the global gauge invariance. 
However, as already indicated above, the continuation of E(z } z*) does not exist and hence 
the treatment of ref. [78, 79] has to be considered as a physically motivated definition but not 
as a mathematically correct derivation of the Minkowski energy functional. 

Very recently it has been demonstrated that counting powers of the time components of 
vector fields in the eigenvalues of the Dirac Hamiltonian does not correctly reproduce the 
corresponding expansion for Ar [71]. The reason is that the manipulations leading to the sum 

(6.24) are in fact ill-defined as already indicated because the traces of h and are computed 
in different Hilbert spaces. This topic will further be illuminated in subsection 6.2.4. 



6.2 Self— consistent solutions 



We denote meson field configurations which extremize the Minkowski energy functional 
(6.39) self-consistent. Let <p represent the whole set of meson fields involved in the soliton 
configuration. Then the stationary condition may be expressed as 

6E 

— = 0. 6.40 

The mesonic part of the energy functional, E m , is at least quadratic in <p while the Dirac 
operator provides a Yukawa type coupling of <p to the quark fields. Thus (6.40) relates the 
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meson fields to the quark fields. The latter diagonalize the static Hamiltonian h and the as- 
sociated eigenvalues again enter the energy functional. It is this interplay between eigenstates 
and -values of h which leads to the notion of self-consistency. Commonly the self-consistent 
soliton solution is numerically obtained with the use of an iterative procedure [83, 80]. A test 
profile is employed to diagonalize h. The resulting eigenvalues and -vectors are then substi- 
tuted into the stationary condition (6.40) yielding an updated profile function. This updated 
profile serves again as input for h. This procedure is repeated until convergence is gained. 



6.2.1 The pseudoscalar hedgehog 



The simplest meson field configuration which allows for soliton solutions involves the chiral 
field only. This configuration has already been employed in chapter 5 to discuss general 
properties of the soliton. Then all vector meson fields are set to zero and for the pseudoscalar 
field the celebrated hedgehog ansatz is assumed 

M = m exp (it ■ r 0(r)) . (6.41) 

Obviously the scalar fields are constrained to the chiral circle E = (E) = m. The radial 
function 0(r) is commonly referred to as the chiral angle. Substituting this ansatz into the 
Dirac Hamiltonian (6.16) yields 

h = ol ■ p + f3m (cos0(r) + ijsT ■ r sin0(r)) . (6.42) 

As h is Hermitian the eigenvalues e v obtained from 

h^ v = t v y v (6.43) 

are real. Technically, the discretization of the eigenvalues e v is achieved by restricting the 
coordinate space M 3 to a spherical cavity of radius D and demanding certain boundary condi- 
tions at r = D. Eventually the continuum limit D — > oo has to be considered. We relegate the 
discussion of the special form of the boundary conditions to appendix B where the coordinate 
representations of the states \v) are given. It is nevertheless worthwhile to discuss the struc- 
ture of these states. Due to the special form of the hedgehog (6.41) the Dirac Hamiltonian 
commutes with the grand spin operator 

G = j + ^ = l+° + ^ (6.44) 

where j labels the total spin and I the orbital angular momentum, r/2 and cr/2 denote the 
isospin and spin operators, respectively. The eigenstates of h are then as well eigenstates of G. 
The latter are constructed by first coupling spin and orbital angular momentum to the total 
spin which is subsequently coupled with the isospin to the grand spin [81]. The resulting states 
are denoted by \ljGM) with M being the projection of G. These states obey the selection 
rules 

G + l/2, /={5! +1 
j = { „ • (6-45) 



G-l/2, 1={ G G _ X 
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As the one-particle eigenenergies e v are real, the Euclidean energy does not possess an 
imaginary part. The Minkowski energy functional arising from eq (6.27) is then given by 




(6.46) 



with the mesonic part 



E m = ^mlf^ / drr 2 (1 - cos©(r)) . 



(6.47) 



Also the energy functional associated with the trivial meson configuration © = is subtracted 
in (6.46) and (6.47). The stationary condition 6E[Q]/ 6Q(r) is obtained with the help of the 
chain rule. Whence we require the functional derivative of e v 



St 



6Q( 



m 



r)f3 (— sin0(r) + ijsT ■ r cos0(r)) ^ v (r) 



(6.48) 



Then the stationary condition becomes the equation of motion for [83] 

4tt m 2 Jl ' 



cos0(r) tr J dVl ps(r, r)i^f 5 T ■ r = sin0(r) < tr J dVl ps(r, r) 



Nn m 



(6.49) 



where the traces are over flavor and Dirac indices only. According to the sum (6.46) the scalar 
quark density matrix ps(x } y) = (q(x)q(y)) is decomposed into valence quark and Dirac sea 
parts: 



ps{x,y) 



p v s(*,v) + P°s(x,y) 



p v s( x ,y) 



ps{ x ^y) 



V 

-1 



* J/ (y)sgn(e J/ 



(6.50) 



As indicated above the stationary condition relates the soliton profile function 0(r) to the 
eigenfunctions ^ v of the Dirac Hamiltonian h. 

We are now at the point to discuss the numerical solutions in the unit baryon number sector 
(B = 1). For the eigenvectors (B.4-B.6) one can easily verify that tr / dft ps{f^ r)i^ 5 T-r = 
0. This transfers to the boundary condition for the chiral angle: 0(0) = lir. On the other 
hand, for r — > oo (6.49) has a solution with = 0. In the (B = 1) sector it turns out that 
the soliton profiles with the boundary conditions 0(0) = —it and ©(oo) = (modulu 2tt) 
extremize the static energy functional (6.46). It should be mentioned that for configurations 
which do not obey these boundary conditions the energy remains finite although it is not a 
minimum. This is in contrast to topological soliton models [82] where the energy diverges for 
infinitesimal deviations from these boundary conditions. The NJL soliton is not a topological 
one! 
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Figure 6.1: The radial dependence of the self-consistent profile function 0(r) for various 
constituent quark masses m. The pion mass is m v = 135 MeV. 

Table 6.1: The soliton energy E tot as well as its various contributions according to the sum 
(6.46) as functions of the constituent quark mass m. The pion mass is taken to be m v = 
135MeV. All numbers are in MeV. 



m 


350 


400 


500 


600 


700 


800 


Etot 


1236 


1239 


1221 


1193 


1161 


1130 


E v 


745 


633 


460 


293 


121 


-55 


E 


459 


571 


728 


869 


1012 


1103 




31 


34 


33 


31 


28 


26 



Self-consistent solutions were numerically obtained for m > 325MeV [83]. In Figure 6.1 
the radial dependence of the self-consistent profile is displayed for various values of the con- 
stituent quark mass m. The profile function obviously exhibits only a very mild dependence 
on m. The soliton energy -Etot, the energy functional corresponding to the self-consistent 
chiral angle, shows the same characteristics as can be seen from table 6.1. Obviously the 
soliton mass is larger than the three quark threshold as long as m<420MeV. In table 6.1 also 
the various contributions to the soliton energy are displayed. For small constituent quark 
masses the explicit occupation of the valence quark orbit (?y va i = 1) provides the dominant 
contribution to the energy. As m increases the situation is reversed. For m > 750MeV Ey 
even becomes negative. According to our previous discussions this implies ?y va i = 0. Stated 
otherwise: The valence quark becomes part of the polarized vacuum which then carries the 
baryon number. This situation actually corresponds to the Skyrmion picture of the baryon. 
Furthermore, the contribution E m of the purely mesonic part of the action to the energy is 
shown in table (6.1). This quantity can actually be interpreted as the pion-nucleon E term 
at zero momentum transfer E 7r jv(<? 2 = 0). The results for E m seem to be somewhat lower 
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than the data extracted from a critical examination of the existing data on ttN scattering 
[84]: E jr jv((j' 2 = 0) 45MeV. It should, however, be noted that T l7rN (q 2 ) is claimed [84] to 
have a strong momentum dependence which makes the extraction of the E^jvl*? 2 = 0) from 
ttN scattering data rather involved. 



6.2.2 Beyond the chiral circle 



So far the scalar mesons have been constrained to their vacuum expectation value, the 
constituent quark mass m. The reason for choosing this ansatz is not only that of simplification 
but rather the appearance of a special instability of the soliton configuration when the scalar 
field is allowed to be space dependent; at least in the proper time regularization scheme 
[77, 85]. In this case one can prove the existence of a meson field configuration, which has 
vanishing energy. When varying the spatial extension of this configuration below a critical size 
the system passes from a unit baryon number to a zero baryon number configuration. This 
variation can be parametrized as (M = a + 7r • r) 

a(r) = <M1 + Wf(^) cos 0(-^)], 7r(r) = r<f> Wf(^) sin ©(-^) (6.51) 

with W • (R(f>o) a = const. In ref. [85] a Wood-Saxon shape was chosen for f(r/R) and a linear 
profile for 0(r/i?). The constant <f> ~ m characterizes the vacuum expectation value of the 
scalar field. For very narrow configurations, i.e. R — > 0, the mesonic part of the NJL soliton 
energy, which is maximally quadratic in the scalar field <f> = \J a 2 + 7r 2 (to which we will 
refer to as the chiral radius), can be shown to vanish for a < 3/2. Thus only the fermion 
determinant contributes to the energy in this case. For such a value of a the spectrum of 
the Dirac Hamiltonian has the interesting feature that the valence quark level gets transferred 
from the lower boundary of the positive Dirac spectrum to the upper boundary of negative 
Dirac spectrum. All other levels in the G v = + channel follow as a consequence of "avoided 
crossings" [85]. Obviously the baryon number is carried by the Dirac sea and hence the soliton 
energy is solely given by the vacuum part E . This is plotted as a function of R- <f> in figure 6.2 
for a = 4/3. Evidently the total energy is identical to zero in the limit R — > 0. However, since 
these localized meson field configurations carry no baryon number for small i?, this instability 
is unphysical indicating merely that for R — > the system changes from the B = 1 to the 
B = sector. This is also indicated by the fact that these configurations assume zero energy, 
which is the energy of the ground state in the B = sector. 

One way to circumwent these problems is e.g. indicated the possibility to mock up the 
trace anomaly of QCD in effective meson theories. In such an approach the scalar dilaton 
field x ~ (Gfii/G^) is incorporated as an order parameter to absorb the mass dimension of 
the parameters [87]. Especially one may introduce a coupling to the chiral radius <f> such that 
a fourth order term appears in the meson part of the energy [88, 86] 

E m = 4tt j drr 2 a 2 (^-(f> 4 (r) — ?7Xo^ 2 ( r ) — const^ . (6.52) 

Here Xo denotes the vacuum expectation value of the dilaton field. The two new parameters 
a and r\ are related via the (extended) gap equation. Hence one may equally well consider 
rjXo as the only free parameter (besides the constituent quark mass m). The meson energy 
with the dilation included (6.52) only vanishes for a < 3/4 when the parametrization (6.51) is 
adopted [86]. Glancing again at figure 6.2 one observes that in this case the energy does not 
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Figure 6.2: The energy of the fermion determinant, Ey + E 0} for the variational meson held 
conhguration (6.51). The parameter a is dehned after eq (6.51). (Figure taken from ref.[86].) 




Figure 6.3: The self-consistent soliton conhguration with the scalar held included. "Non- 
linear" refers to the case with the scalar helds hxed at their vacuum expectation values. 
(Figure taken from ref. [86].) 



56 



vanish for R — > but rather becomes very large. Indeed, the solution to the self-consistent 
problem exists and is displayed in figure 6.3. From this figure one also recognizes that the 
collapse appears when the fourth order interaction is switched off. As a consequence of the 
gap equation the limit a — > corresponds to ^Xo/^o ~~ > Figure 6.3 also demonstrates 

that for increasing coupling to the dilaton field, the chiral angle, gets more concentrated at 
r = 0. This reflects the bag formation caused by the dilaton field [89]. 

The above considerations made use of an extended definition of the energy functional in 
order to avoid the collapse of the meson profiles. As the net result of the limit R — > in the 
configuration (6.51) for a = 3/4 is to shift all eigenvalues e v in the G v = + channel by one 
level, the vanishing energy is just a matter of regularization which ignores the asymmetries of 
the spectrum of the Dirac Hamiltonian at large energies. However, the baryon number (6.12) 
is sensitive to these large energies because no regularization is involved. Taking the point of 
view that the imaginary part of the Euclidean action should also be regularized, automatically 
leads to a regularized baryon number 



In the same way as the energy vanishes, goes to zero for the above considered limit R — > 0. 
In ref. [90] therefore the idea has been pursued to constrain B^ to unity in order to avoid the 
collapse. Then the energy functional (6.46) is supplemented by a Lagrange multiplier for B^ 



In this case no quartic term in <f> is needed in the mesonic part of the action. The last term 
in (6.54) is introduced for convenience; eventually the limit ^ — > should be assumed. At 
first place it should be noted that for the chiral soliton with the scalar fields constrained to 
the vacuum expectation values (see subsection 6.2.1) B^ is not exactly unity; but rather only 
0.97. Even by including the constraint (6.54) this number cannot be increased significantly as 
long as <f> is not allowed to be space dependent. Taking <f> = <f>(r) to be a radial function indeed 
increases B^. Numerically, however, a solution with B^ = 1 has not been found for finite^ D. 
For finite D always a non-vanishing lower limit for ( exists below which no solution has been 
found. This limit can be translated into a value for B^ which slightly deviates from unity 
but approaches it as D increases. The associated self-consistent soliton profiles are displayed 
in figure 6.4. The most important result is, of course, that stable solutions do exist when 
the regularized baryon number is fixed. Furthermore, one observes that the large distance 
(r > Ifm) behavior of the profiles is almost uneffected by the constraint (6.54). However, in 
the vicinity of the origin drastic changes occur. For B^ — > 1 the chiral radius, <f> vanishes at 
the origin, i.e. the chiral symmetry is restored. Simultaneously, the slope of the chiral angle, 
increases and eventually goes to infinity. In ref. [90] it has also been shown that constraining 
the regularized baryon number leads to a valence quark dominated picture of the soliton. 
Especially, the corresponding energy eigenvalue approaches the one obtained in the baryon 
number zero sector. Nevertheless, the valence quark wave-function was found to be strongly 
localized. 

A further approach to prevent the meson profiles from collapsing is to include vector mesons 
[79, 75]; especially because the cu-meson provides a sizable repulsion which is supposed to 
stabilize the system. This will also be discussed in the proceeding subsection. 

d D denotes the radius of the spherical box for the numerical calculations, cf. the discussion after eq (6.43). 




(6.53) 




(6.54) 
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Figure 6.4: The self-consistent soliton configuration with the scalar held included and the 
regularized baryon number (6.53) being constrained. Different values for correspond 
to various radii of the spherical box used for the numerical calculations. "Non-linear" denotes 
the case without constraint and the scalar held at its vacuum expectation values. (Figure 
taken from ref. [90].) 

6.2.3 Inclusion of (axial) vector mesons 

In recent years the NJL soliton has experienced several types of extensions. In the present 
section we will describe the effects of various vector and axial-vector meson fields on the soliton 
solution and put some emphasis on the historical development. Proceeding in this manner 
we are also enabled to illuminate the effects of different (axial) vector meson fields separately. 
A convenient parametrization of the anti-Hermitian fields and is given by separating 
isoscalar and isovector parts 

T T 

V tl = iu tl + ip tl - - , A tl = if lfl + ia ltl - -. (6.55) 

The Euclidean Dirac Hamiltonian for the most general grand spin zero held configuration 
consistent with the parity properties of the meson fields is given in eq (6.63). Additionally the 
complete set equations of motion resulting from extremizing the Minkowski energy functional 
(6.39) is displayed in appendix C. The Hamiltonian corresponding to the subsystems which 
will subsequently be discussed below can easily be obtained from (6.63) by setting those fields, 
which are not involved, to their vacuum expectation values. The detailed numerical results 
will be presented at the end of this section. This allows us to directly compare the influence 
of the various fields. 

In ref. [91] the inclusion of the /j-meson via the Wu-Yang ansatz 

Pia = tikJ-kG{r) (6.56) 

has been studied. Except for the pseudoscalar field, for which the hedgehog shape (6.41) 
was adopted, all other fields were set to their vacuum expectation values. Although this 
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configuration does not exhibit global chiral invariance it allows one to investigate the attractive 
character of the /j-meson. Indeed, the size of the soliton, as measured by the baryonic root 
mean squared radius, was found to decrease. Furthermore the value for m, above which soliton 
solutions do exist, decreased to 270MeV. It should be stressed that the incorporation of the 
/j-meson neither alters any of the parameters involved nor introduces additional ones. This, 
however, is not the case when also the ai-meson is included. As the tt — mixing changes 
the relation between the pion decay constant, /„., and the cut-off, A, the value of A for a 
given constituent quark mass, m, is significantly increased [31], see also eq (4.53). For details 
the reader is referred to chapter 4. Needless to remark that the inclusion of the a\ meson is 
mandatory to preserve chiral symmetry. The ansatz for the a\ meson which is consistent with 
the grand spin zero assumption for the static fields and exhibits the pseudovector character 
involves two radial functions 



Self-consistent solutions were found for m > 300MeV [92, 93]. From the equation of motion 
for H(r) (cf. eq (C.9)) one easily verifies that -H"(0) is non-vanishing. Thus the axial-vector 
meson has a direct influence on the valence quark wave-function which is localized at the 
origin. As a result, the energy of the valence quark state is significantly lowered and actually 
becomes negative. From the discussion in section 6.1 it is then obvious that the total baryon 
charge is carried by the polarized Dirac sea and no orbit is explicitly occupied, i.e. all r\ v = 0, 
in the unit baryon number sector. The baryon charge is solely due to j£ in eq (6.11). This 
result strongly supports the Skyrmion picture of baryons as can be seen from the following 
considerations. The gradient expansion of j£ yields in leading order the topological current 
(4.65) [94] which, in the Skyrme model, is identified with the baryon current 6 . A priori the 
point of view that a distorted Dirac sea is responsible for a non-vanishing baryon charge is 
integrated in the Skyrme model. This point of opinion is commonly referred to as Witten's 
conjecture. Obviously, it is supported by the NJL model [92]. Furthermore, the Skyrmion 
phenomenology has established that short range effects are described by either explicit quark 
degrees of freedom or (axial) vector mesons [95, 96]; the "or" being exclusive. In the NJL 
model, however, this matter of opinion represents a direct result! 

Up to now the cu-meson has been left out, mostly for technical reasons. As discussed 
in section 6.1 it introduces a non-Hermitian Dirac Hamiltonian (c/. eq (6.63)) and, more 
troublesome, in a non-perturbative approach a unique extraction of a Minkowski energy func- 
tional is not available. However, the incorporation of this meson is extremely desirable from 
the point of vector meson dominance (the cu-meson has significant influence on the isoscalar 
radius) and the current field identities (3.21) (the isoscalar part of the electromagnetic current 
is proportional to the u field). Additionally the repulsive character of the cu-meson is expected 
to provide stability even when the scalar field is allowed to vary in space. We will therefore 
report on investigations which are based on the physically motivated definition (6.39) for the 
Minkowski energy functional. 

First attempts in this direction have concerned the tt — u system with all other fields put 
to their vacuum expectation values [78, 76]. Only a radial function for the time-component 
of the isoscalar vector meson is allowed by the grand spin symmetry 



e The topological current (4.65) can also be obtained as the Noether current associated with the singlet 
vector symmetry once the Skyrme model is extended to flavor SU(3) and the Wess-Zumino term is included. 




(6.57) 
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As the Hamiltonian is non-Hermitian we have to distinguish between left and right eigenstates 

h\m u ) = e u \m u ) ^l v \h = t v ^! v \ i.e.ht\y v ) = tl\y v ) (6.59) 

with the normalization condition (v^l*^) = 6^. In general it is possible to separate the 
non-Hermitian part by writing 

h = h @ + iuj (6.60) 

with h@ being Hermitian. Assuming the basis (B.4,B.5) it is easy to see that the matrix 
elements of u are real and symmetric. This further implies \^ v ) = It is then useful 

to decompose the eigenstates into real and imaginary parts \® v ) = |*^) + H ere the 

superscript refers to real (R) and imaginary (I) parts of the expansion coefficients V v k in eq 
(B.6). Taking account of the fact that 1^^) = — allows one to extract real and 

imaginary parts of the one-particle energy eigenvalues 

£ = (*^e|*f)-(^|^e|^)-(^M*f)-(*fM^), 

< = (*f|o;|^)-(^|a;|^) + (^^e|^) + (^^e|^). (6.61) 

These relations are suitable to evaluate the functional derivatives of e^' 1 ) with respect to the 
meson fields. E.g. one finds 

-0^ = r 2 /^ (<r|*?)<*» - <r|^)<*») . (6.62) 

Expressions like this enter the equations of motion for the soliton profile functions (c/. ap- 
pendix C). As can be inferred from the equation of motion for the cu-meson (C.6) this profile 
function is directly related to the baryon density. As a matter of fact this causes a non- 
vanishing to field. In particular, the integral / d 3 r u = gy/Aniy is fixed when Ai is not 
regularized and satisfies the normalization condition on the u field imposed by the current 
field identities (3.21). 

Numerically, solutions have been found for m > 350MeV, although the determination of a 
lower bound has not been the central issue in ref. [76]. More importantly it has to be remarked 
that these solutions were obtained for the physical u meson mass, m w = 770MeV. This is in 
contrast to the treatment discussed in ref. [78]-^ where stable solutions appeared only when m u 
was chosen about four times as large. In subsection 6.2.4 we give a possible explanation for 
the non-existence of solutions in that treatment. 

Next it should be noted that the value of the constituent quark mass m at which the 
valence quark energy changes its sign is at 545MeV and thus considerably lower than in the 
purely pionic system. This result is easy to understand since the repulsive character of the 
to meson yields a larger extension of the soliton profile which causes the valence quark to 
be more strongly bound. This repulsive effect is also observed when evaluating the baryonic 
radius although its extraction is somewhat hampered by finite size effects. These effects do 
not show up in the determination of the energy, however, higher moments of the u profile 
function may be obscured. It has also been ascertained that regularization of the imaginary 
part of the action does not lead to qualitative changes of the above presented results. 

In the next step all vector mesons were included [73]. This embraces the ansdtze (6.56,6.57) 
and (6.58), only the isoscalar-scalar field being kept at its vacuum expectation value. The 

•^See also section 6.1 for the discussion of this approach. 
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above discussed structure of the Euclidean Dirac Hamiltonian is not altered, however, it is 
more complex because a larger number of fields is involved. In this case the interesting 
question arises whether or not Witten's conjectures remains valid. The previous explorations 
seem to indicate that effects on the valence quark energy eigenvalue of the isovector and 
isoscalar (axial) vector mesons add up coherently. Numerically self-consistent solutions to 
this extended problem have been found in the interval 300MeV< m <400MeV although the 
authors of ref. [73] do not exclude the existence of solutions in a larger range. Indeed a further 
decrease of the valence quark energy eigenvalue, e^j has been observed, e.g. for m = 350MeV 
e^j changes from -133MeV to -154MeV when supplementing the it — p — a\ system by the to 
meson. This, however, is only the case when the imaginary part of the action is subject to 
regularization; when this regularization is discared e^j increases even slightly as an effect of 
including u. On the whole, however, the influence of regularizing the imaginary part has been 
found to be small. 

Very recently the energy functional (6.39) has been studied for the case that also the 
isoscalar-scalar field is space dependent [75]. This introduces the additional radial function 
<f>(r), cf. subsection 6.2.2. As this represents the most general case we use this opportunity to 
display the associated Euclidean Dirac Hamiltonian 

h = a ■ p + iio(r) + m<f>(r)/3 (cos0(r) + i^T ■ r sin0(r)) 

+ ^((x x r) ■ rG{r) + ^(<r • f)(r • r)F(r) + ^(<r • r)H{r). (6.63) 

As long as the scalar field <f> has not been included as a dynamical degree of the freedom 
the question of whether or not regularizing the imaginary part only played a subleading role. 
Taking, however, <f> to be space dependent leads to a completely different situation. As al- 
ready seen in subsection 6.2.2 a field configuration which is sharply peaked at r = yields 
a vanishing soliton energy. This remains true even when vector mesons are included as long 
as no mechanism prevents these vector mesons from being zero. There is no such mechanism 
for the static p and a\ fields but the static u field is directly correlated to the baryon number 
density. Unfortunately, for the sharply peaked field configuration (6.51) also the regularized 
baryon number vanishes and thus also u = is allowed in this case. 

When the imaginary part is not regularized stable solutions have been obtained for m>350 
MeV [75]. Then the u meson provides sufficient repulsion to prevent the scalar field from 
collapsing. The lower bound for m is somewhat larger (410MeV) when p and a\ fields are 
ignored. Actually the deviation of the scalar field from its vacuum expectation values is almost 
negligible. The alternative method of stabilization (6.54) has not been studied for the case 
that ( 1) vector mesons are included. 

In table 6.2 we summarize the numerical results for the various cases discussed above. 
There the specific contributions to the energy are compared for the constituent quark mass 
m = 350MeV. 

Figure 6.5 is well suited to discuss the effects of various mesons fields on the chiral angle. 
One observes that the p meson is repulsive at small distances and attractive for r > 0.3fm. A 
much stronger attraction is provided by the axial-vector meson. Also the repulsive character 
of the cu-meson is reproduced by the NJL model. In the inner region of the soliton <1.5fm 
the scalar field provides a small repulsion while at larger distances (not shown in figure 6.5) 
the scalar field causes the chiral angle to swell slightly. Although this attraction appears to be 
quite weak it indicates that the NJL soliton with scalar mesons contains at least some inter- 
mediate range attraction called for by the central potential of the nucleon-nucleon interaction. 
Unfortunately, no quantitative statement on this subject can be made presently in the NJL 
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Table 6.2: The soliton energy for various treatments of the NJL soliton. The meson fields listed 
in the first line represent those meson profiles which are allowed to deviate from their vacuum 
expectation values. All numbers are evaluated for a constituent quark mass m = 350MeV and 
m v = 135MeV. 
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7T, id 


7r, p, a x 


7T, LO, p, CL\ 


<f>, 7r,cu, p, a x 


E (MeV) 


1236 


966 


1404 


1011 


1144 


1125 
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459 
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El (MeV) 








-14 





186 


177 


E m (MeV) 


31 


173 





395 


400 


-1323 


e val/ m 


0.71 


0.20 


0.56 


-0.37 


-0.27 


-0.28 


4J m 








0.27 





0.16 


0.15 



model. 

6.2.4 Quadratic expansion for time components of vector fields 

The preceding discussions concerning the analytic continuation of the action in the pres- 
ence of the to meson were based on the analysis of the analytic properties of the eigenvalues 
of the Euclidean Dirac Hamiltonian h (6.63). It has been assumed that the expansion of 
these eigenvalues in powers of the cu-field can straightforwardly be transferred to the action 
functional. However, it has recently been demonstrated that this is actually not the case 
[71]. In order to prove the identity one had to assume that h and can be diagonalized 
simultaneously. Starting off at 

PeA = -d 2 T + h% + 2iuj 4 d T + i [h @ ,u 4 ] + u\ (6.64) 
and imposing the proper-time prescription at the operator level 

1 f°° ds ( I 
A r ^--Tt — exp\-d 2 T + hl + 2iuj 4 d T + i[h @ ,uJ4]+^l} (6.65) 

2 Jl/A 2 s L J 

the energy functional (6.25) can only be obtained with the above mentioned assumption of 

simultaneous diagonalizability because the term linear in d T originates from (h — d T . One 

may obtain eq (6.25) from (6.65) by substituting the coefficient of d T with 2ze^. The before- 
mentioned ambiguities, which arise from the inadequate application of the rules for manipu- 
lating the logarithm to derive (6.24), are avoided when starting from (6.65) without further 
assumptions. It is also obvious that because of these assumptions different expansion schemes 
will lead to different results. It is the purpose of the present section to point out the differences 
between the approach which relies on expanding the eigenvalues e v of h on the one hand and 
the expansion of the operator (6.65) on the other; both expansions are understood in terms 
of id 4. 

Employing techniques which have been worked out in the context of the semi-classical 
quantization of the soliton (section 7.1) [70] and been extended to the treatment of small 
amplitude fluctuations off the soliton (see section 7.3) [97] one obtains a Minkowski energy 
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r/ f m 



Figure 6.5: The profile function 0(r) for the self-consistent chiral angle in various approaches 
to the NJL model. Those fields which are allowed to be space dependent are indicated. The 
parameters used are m = 350MeV and m v = f35MeV. 



functional up to second order in the u held [71] 



£f mk [O,io ] = £ [©] + iVc?7vai(val|coo|val) 



c 



^sgn(e°)erfc 



(v\uj \v) (6.66) 



-N C f] va \ 



| (v\loo |val 



Nr 



^val 



2 47T m 2 , fi 

Z LUJ7T 



drr 2 ul. 



Here -Eo[©] refers to the energy functional associated with the Hermitian part of the Hamil- 
tonian h@. Furthermore e° and \v) denote the eigenvalues and -vectors of h@. Also the mass 
terms associated with the p and a\ mesons are included in _E [0]. The imaginary part of the 
action has been assumed in regularized form. If one wishes to abandon this regularization the 
complementary error function in (6.66) has to be replaced by unity. Finally the regularization 
function 



/(e^e*; A) 



1 sgn^erfc ( e f)~ sgn^erfc ( ^ ) A e"^) 2 - e -( £ "/ A ) 



7T 



nas a smooth limit [70] 



lim / (e^, e v ; A) = 0. 



(6.67) 



(6.68) 



It is actually because of this limit that the consistency conditions for energy functional, which 
are discussed in section 6.1, are satisfied by the functional (6.66). As was to be expected this 
regulator function is identical to the one for the moment of inertia [70]. In both cases the real 
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Energy 

m. = 770 MeV 
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Meson Profiles 

m = 500 MeV, A=0.73 
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Figure 6.6: Left: The total energy as a function of the scaling parameter A defined in eq 
(6.69) for various values of the constituent quark mass m. Right: The profile functions which 
minimize the energy functional (6.66) under variation of the scaling parameter A. The quark 
baryonic density b ~ q is artificially scaled such that the spatial integrals over uo/m and 
b coincide. Here the constituent quark mass is assumed to be 500MeV. (Figure taken from 
ref.[7f].) 

part of the action is expanded in terms of static and anti-Hermitian operators (in the case of 
the moment of inertia these are the isospin generators). 

In ref. [71] several parametrical descriptions have been adopted for the profiles other than 
uo(r) and subsequently the latter has been determined from the stationary condition corre- 
sponding to (6.66). As an example we quote the scaling ansatz 

A (r) = s . c .(Ar) (6.69) 

for the model consisting only of the pion and the to meson. In eq (6.69) s . c .( r ) refers to 
the self-consistent soliton solution to the problem without the to meson. Subsequently the 
stationary condition associated with the functional (6.66) was solved in order to obtain the 
profile function uo\(r). This procedure of first parametrizing the chiral angle and subsequently 
solving for the to profile exactly, is also justified by the fact that commonly the stationary 
condition for uo(r) represents a constraint rather than only an Euler-Lagrange equation of 
motion. Then the profile functions as well as the energy functional depend on the parameter 
A. This dependence is displayed in figure 6.6. This figure indicates that the functional (6.66) 
indeed possesses a unique local minimum for constituent quark masses m>400MeV. This 
figure also contains the profile function which minimize (6.66) when m = 500MeV is adopted. 
Obviously the to profile gets support at larger distances than does the quark baryon current 
b ~ q. This confirms the repulsive character of the uo meson. The minimal energy values are 
given as functions of m in table 6.3. Also shown are the contributions due to the linear and 
quadratic terms in uo. Note the virial factor, -2, between these two pieces. The presence of the 
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Table 6.3: Soliton energy 

^Mink (g_ 66 ) when q and 

to are the only space dependent fields 
for various constituent quark masses. Also shown are the contributions from the terms linear 
(lin.) and quadratic (quad.) in u. All numbers are in MeV. (Results are taken from ref. [71].) 



m 


^ Mink 


lin. 


quad. 


400 


1585 


650 


-325 


450 


1685 


734 


-367 


500 


1769 


856 


-428 



to meson then increases the soliton energy by about 300 — 400MeV. This contribution grows 
with the constituent quark mass. It should be noted that the quadratic term splits into two 
pieces, one from the purely mesonic part of the action (— 296MeV for m = 500MeV) and one 
stemming from the expansion of Ar ( — 129MeV for m = 500MeV). The smallness of the last 
quantity compared to the total energy gives evidence for the fact that this expansion of Ar 
in terms of cu 4 converges quickly. Furthermore, it has been shown that including the u meson 
via the expansion (6.66) indeed provides a sizable repulsion, a feature commonly attributed 
to this meson field. This is expressed by the fact that the value, which A assumes at the 
minimum of the energy, is significantly smaller than one. 

The approach having just been described has to be compared with an expansion of the 
eigenvalues of the whole Dirac Hamiltonian h@ + zcu 4 (6.60) in Euclidean space. Applying 
standard perturbation techniques results in 

r v-^ H^4^)(H^4k) , / 7 1 1 \ (a 7n\ 

2s ^Z~o ••• and t v = i{v\u) A \v) . . . . (6.70) 

Substitution of this expansion into the expression for real and imaginary parts of the Euclidean 
energy (6.27 and 6.36) and subsequent continuation to Minkowski space yields at second order 
an energy functional similar to (6.66), however, with / (e^, t v \ A) replaced by 

{1 sgn(e I/ )erfc(|5f|)-sgn(e M )erfc(|^ i |) 
2 ^ (6.71) 
0, e u = 

We will refer to the corresponding energy functional by E 2 . In ref. [71] it has been demonstrated 
that for a well motivated field configuration E 2 deviates only by about 3% from (6.37) which 
results from diagonalizing (6.60). 

However, the regulator functions / and / differ significantly. First of all, they only agree 
in the limit A — > 00. As the cut-off A is quite small this difference is sizable numerically. Next 
it has to be noted that / is discontinuous as — > e v . This has the unpleasant consequence 
that the u contribution to the energy is not positive definite for configurations which satisfy 
the stationary condition. Furthermore, the second order contribution associated with / is 
— 553MeV; to be compared with the above mentioned — 129MeV for /. Most importantly, 
however, E 2 does not possess a local minimum. This explains why for empirical parameters 
the authors of ref. [78] did not find a self-consistent solution in the model containing only the 
pion and the u degrees of freedom. 

These investigations demonstrate that counting powers of cu 4 in the Euclidean energy func- 
tional (6.37) may not be the correct approach to derive a Minkowski energy functional in 
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the presence of isoscalar field u. Rather the functional (6.66) appears to be the appropriate 
starting point. This statement is also supported by previous calculations in the context of 
the semi-classical quantization of the chiral soliton as well as other nucleon observables (see 
chapter 7). These computations all impose the proper-time regularization at the action level 
in Euclidean space (6.65). 



6.2.5 Local chiral rotation 



In sections 3.4 and 4.2 the effectiveness of the chiral transformation ^ = with T = 
PlCl + PrCr has become evident when exploring the meson sector of the NJL model. In 
particular, it provides the link between the massive Yang-Mills and hidden gauge approaches 
to vector mesons. It is therefore interesting to illuminate the role of this transformation in 

the soliton sector as well. For this study the unitary gauge Cr = Cl = C wm be adopted. This 
implies T(0) = cos(0/2) + i^T ■ r sin(0/2). The Dirac Hamiltonian for the chirally rotated 
quark fields ^ then becomes [32, 98] 

h R = T(Q)hT^(Q) = a-p + Pm-^((r-r)(T-r)(Q'[r 




-^-(<r • r) sin 0(r) - ■ (r x r) sin 2 ( ) . (6.72) 

Again one observes that the chiral field has been eliminated from the mass term at the expense 
of induced vector mesons. Before discussing the physical implications of this transformation 
in the soliton sector a few remarks on its technical feasibility are in order. First of all, one 
observes that the coordinate singularities do not disappear at r = as long as 0(r = 0) = lir. 
These singularities appear since T(0) is topologically distinct from the unit transformation. 
For T(20) these coordinate singularities are absent, however, the spectrum obtained from 

diagonalizing T(20)/Wt(20) in a finite basis is not identical to that of h. In the G v = + 
channel the lowest state is missing while an additional one shows up at the upper boundary 
in momentum space. This shift of eigenstates is repeated for each n in T(2n0), with n > 2. 
Hence, one observes another reflection of the topological character of the transformation T. 
Secondly, one should remember that for the diagonalization (6.43) certain boundary conditions 
on the eigenstates have been imposed (c/. the discussion after eq (B.6)). It should be obvious 
that these boundary conditions are affected by the chiral rotation T which is not uniquely 
defined at r = 0. With the corresponding redefinition of the basis spinors the coordinate 
singularities are canceled. The numerical diagonalization 

h R y v = (6.73) 

then indeed yields e u = e u and = Tt(0)vl/ J/ with e u and ^ being the solutions to the 
eigenvalue problem (6.43). 

This formalism has in turn been employed to examine the situation when the chiral rotation 
is performed on the ( 1) vector fields [31] 

%-K = t L (d, + V,-A,)d. (6.74) 
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Table 6.4: Contributions to the energy for self-consistent solution in various treatments of 
the axial-vector meson in the NJL model. Those meson fields which are allowed to be space 
dependent are indicated. The constituent quark mass m =400MeV is common. All numbers 
are in MeV. (Results are taken from ref. [98].) 





TT — p 


7r — p — a\ 


7T - V 


£val 


313 


-222 


-351 


J^det 


711 


543 


240 


E m 


149 


393 


175 


Etot 


861 


937 


415 



and the transformed axial-vector field is neglected, i.e. = 0. This ansatz does not violate 
the local chiral symmetry in contrast to the approach = 0, cf. eq (3.42). This configuration 
is of special interest since this approximation has frequently been applied to extended Skyrme 
models 3 . The important question arises which (or whether at all any) features of the unrotated 
axial-vector field (A^) are maintained. For these studies the effects of the u meson have been 
ignored. After carrying over the chiral rotation onto the quark spinors only vector meson 
degrees of freedom are contained in the Dirac Hamiltonian. For these the Wu-Yang ansatz has 
been assumed yielding [98] 

G(r) 

h R = a ■ p + /3m + — —a ■ (r x r). (6.75) 
2r 

The coupling of the vector and pseudoscalar fields is then completely contained in the mesonic 
part of the energy functional [98] 

E m = ^- J dr |(G(r) + 1 - cos0(r)) 2 + K 2 (O'(r)) 2 + sin 2 0(r)| . (6.76) 

Obviously 0(r = 0) = —it implies that G{r = 0) = —2. One thus has to deal with a Dirac 
Hamiltonian which contains a topologically non-trivial vector meson field. Thus (6.75) is 
singular and cannot be treated using the standard basis [81] but rather by employing techniques 
which are analogous to those developed to diagonalize the rotated Hamiltonian (6.72). 

It has been demonstrated that self-consistent solutions which minimize the total energy 
exist [98]. Here the total energy is defined as the sum E to t = E m + E det with E det being the 
contribution of the fermion determinant (6.27) in terms of the eigenvalues of (6.75). Also the 
valence quark contribution has to be added if necessary to accommodate unit baryon number. 
In table 6.4 the results of these calculations are compared with two different treatments of the 
axial-vector degrees of freedom in the unrotated formulation. Regarding the strong binding 
of the valence quark state it is clear that this important effect, which is commonly asserted to 

^ 0, is retained in the model with A^ = 0. Furthermore, the valence quark wave-functions 
exhibits features which are associated with a localized antiquark dominating lower 

component. Unfortunately the total energy comes out quite low making an application of 
this model for the description of baryons doubtful. In ref. [98] is has also been demonstrated 

3 For a compilation of relevant articles see ref. [12]. 
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that the smallness of the total energy is linked to the missing repulsion in this model. This is 
also reflected by the meson profile functions which possess only a very small spatial extension. 
Again this is an effect originating from the presence of an axial-vector field in the unrotated 
frame, ^ 0. Furthermore the effects asserted to the cu-meson were modeled by the inclusion 
of a 6 th order term (4.22) usually used in the Skyrme model to simulate the cu-meson. Then 
a sizable repulsion was obtained and the energy of the valence quark orbit decreased even 
more to approximately the negative constituent mass [98]. Thus this treatment provides a 
strong support of Skyrme type models which rely on the assumption that the valence quark 
has joined the negative Dirac sea. 
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7 Baryons 



In the preceding sections we have examined the soliton solutions in the NJL model. Except 
for the baryon number these solutions do not carry the quantum numbers of physical baryons. 
In particular, these solutions are neither eigenstates of the angular momentum J nor of isospin 
i\ As these solitons have vanishing grand spin (G = 0) they represent linear combinations of 
states with \J\ = In order to describe physical baryons the soliton configuration has to 
be projected onto states with good angular momentum and isospin quantum numbers. As al- 
ready mentioned in section 5.3 this is commonly achieved by introducing collective coordinates 
which describe the orientation of the soliton in coordinate- and isospace. Subsequently these 
coordinates are canonically quantized yielding the angular momentum and isospin operators. 
As the grand spin symmetry of the hedgehog (6.41) causes the equivalence of rotations in 
coordinate- and isospace only one set of collective coordinates is needed. For convenience, one 
chooses the isospin orientation. This procedure has first been applied by Adkins, Nappi and 
Witten to the SU(2) Skyrmion [9]. In the proceeding section we will explain the analogous 
treatment for the NJL chiral soliton of pseudoscalar fields [70]. The generalization to three 
flavors is more involved not only due to the fact that the moment of inertia tensor (5.18) is no 
longer proportional to the unit matrix but also because of symmetry breaking being present. 
Two complementary approaches ([62] and [99]) will be presented in subsection 7.5. 

7.1 Quantization of the chiral soliton 

The time-dependent collective coordinates describing the isospin orientation of the soliton 
are parametrized with the help of an 2 X 2 unitary matrix R(t) 



where M (r) denotes the static hedgehog configuration (6.41). The scalar-pseudoscalar part 
of A m does not contain any time derivatives. Thus it is independent of collective coordinates 
as long as symmetry breaking is ignored. Then the dependence of the action on R(t) and its 
time derivatives completely originates from the fermion determinant 



The functional trace is most conveniently evaluated by transforming to the flavor rotating 
system q — >■ q' = Rq [70] 








(7.3) 



which introduces the angular velocities 




d 



(7.4) 
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In eq (7.3) h denotes the static Dirac Hamiltonian denned in eq (6.42). In the adiabatic ap- 
proximation, i.e. when the time derivatives of J? are neglected, one might define an intrinsic 
Hamiltonian h' = h-\-^T ■ f2 which could formally be treated analogously to the preceding cal- 
culations as long as regularization is ignored. It is then obvious that the fermion determinant 
can be decomposed into vacuum and valence quark parts. Unfortunately, the eigenvalues of h' 
are not known and thus a perturbation expansion in J? has to be carried out. For the valence 
quark part this is a straightforward application of standard perturbation theory resulting in 

A v = -TE V + I E Qlf^ttj + O (tt 4 ) . (7.5) 

In a symmetric two flavor model only even powers of S7 are allowed. 0™ 1 denotes the valence 
quark contribution to the moment of inertia [70] 

ovai_^v- n _ x (H t »»H t jI^) N c \ - (val|r 8 |/i)(/i|r J |val) 

The latter equality arises because only the valence quark orbit is occupied. The fermion 
determinant gives also rise to a vacuum contribution to the moment of inertia 0™ c [70]. In 
order to extract this part from the fermion determinant one again has to continue to Euclidean 
space and consider the limit T — > oo. In this context it is important to note that J? represents 
the time component of an induced vector field. Thus J? has to be continued as 

Q^riQ E (7.7) 

with Qe being Hermitian. Due to isospin symmetry only even powers of Qe appear in an 
expansion. Thus we only need to consider the real part of the fermion determinant. It is 
useful to define 

K(s, z) = exp (-s(-d T ~ h')(-d T - h')^ \ g ^ = exp (sA(z)) (7.8) 

where we have indicated that the temporal part of the trace is substituted by a spectral integral 
over z for T — > oo. This allows us to express the real part of the fermion determinant 

1 f°° ds f°° dz .„ . 

2 Jl/A 2 S J-oo 2tt 

in the proper-time regularization. In (7.9) the trace includes spatial and internal degrees of 
freedom. Furthermore, 

A(z) = (-d T - h')(-d T - h')^ . = ( z + l -T-n E ) 2 + l[r-f2 E ,h] + h 2 . (7.io) 

In order to extract the moment of inertia we first expand K(s,z) up to second order in Qe 
and define the coefficient of the quadratic term as 



d 2 K(s,z 

(s,z) 



dtt E dtt E 



(7.11) 

J7e=0 
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This derivative may be expressed with the help ol a Feynman parameter integral 



tr (K 2 b (s } z] 



s~ I dx tr 
/o 



lZT a + 



2 ' 
+tr 



) K°(s(l -x),z) (izn + f ,^ ) K°(sx, z) 

(7.12) 



, . K"(s,z 
2 ' 2 ' v ' 



where K°(s, z) = exp (— s(z 2 + h 2 )) denotes the zeroth-order heat kernel. The spectral integral 
in (7.9) is of Gaussian type and may be carried out straightforwardly. The remaining trace is 
performed using the eigenstates of h (6.43). This then allows one to also carry out the Feynman 
parameter integral. The vacuum part of the moment of inertia may finally be extracted from 
A R = -TE - ^Qlt c tt a E tt b E + O (f2 4 E ) in the limit T -> oo [70] 



N c 



(7.13) 



with the cut-off function (which actually is the same as in (6.67)) given by 

A e"WA) a - e -(-M) 2 sgn( ei ,)erfc ( ^ ) - sgn( eA1 )erfc ( c f) 
fe(e ll ,e 1/ ;A) = — — — . (7.14) 



el - el 



2^ - e h 



Due to isospin invariance no direction in isospace is distinct and the total moment of inertia 
is isotropic 



a 



(7.15) 



a 2 is a functional of the chiral angle since the eigenstates \v) and -values e v functionally 
depend on 0. Numerical results for a 2 may e.g. be found in refs.[100, 101] a . Up to quadratic 
order in J? the collective Lagrangian L(f2) in Minkowski space is finally given by 



L(f2) 



-E[Q] + \a 2 [ 



}f2 2 



(7.16) 



One easily verifies that the infinitesimal change of the meson fields M(r,t) under spatial 
rotations may be written as 



[M(r,t),r x d] 



dM(r,t) 



Thus the total spin is given by the Noether charges 



dL 



\ dM c)Q • I c)Q 



a 2 fl. 



(7.17) 



(7.18) 



The hedgehog structure of the soliton relates isospin / and spin J via the adjoint representation 
Da = £tr (t.RtjR^) of R 



Ii — ^ Dij Jj . 



(7.19) 



a Our numerical results agree with ref.[100] but disagree with ref.[101]. 
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In this respect the spin of the soliton may be regarded as the isospin in the isorotating frame 
and vice versa. As in the Skyrme model the resulting collective Hamiltonian is quantized as 
a rigid top 

H = J- — -L = E + — J 2 = E + —I 2 (7.20) 
oil 2a 2 2a 2 

which yields a tower of baryons with identical spin and isospin like e.g. the nucleon or the 
A resonance. We will postpone the discussion of numerical results for the baryon masses 
Mb = E + J ^a^ ^° sec ti° n (7.5.1) where the application of this approach to three flavors 
is discussed. However, the quantization (7.18) is equally important for the study of static 
nucleon properties which will be discussed in the next section. Then frequent use will be 
made of one more relation between collective coordinates and operators [9] 



UJj = -\Dn (7.21) 



which, however, is valid only when sandwiched between nucleon states because then J 2 
I 2 = 3/4. 



7.2 Static nucleon properties 



In this section we will discuss the results for nucleon observables such as magnetic moments 
and charge radii (subsection 7.2.1) and the axial charge of the nucleon, §a (subsection 7.2.2). 
As these quantities correspond to certain moments of symmetry currents the calculations are 
performed in two steps. Firstly, these currents are constructed from the NJL model action 
(3.25). Secondly, matrix elements of these currents with respect to nucleon states are evaluated 
using the apparatus developed in section (7.1). 

The results of such calculations have already been reported in another review article [102]. 
Therefore we only briefly report and comment on these results. The reader may consult that 
reference for more details. 

The symmetry currents Jp(s) are extracted by introducing external gauge fields a T v (x) and 
identifying their linear coupling to the meson fields, see also eqs (3.18) and (3.19). Y denotes 
the symmetries under consideration. These are: (1) the vector symmetry in electromagnetic 
direction Q = ^r 3 + | and (2) the axial symmetry. Formally we write 



6a 



x 



(7.22) 

aY,=0 



where <p denotes the set of meson fields involved. In the present case this only refers to the 
chiral field (7.1). The gauge fields a T v (x) only appear in the fermion determinant 

X( x ) = ^T7^ Trl °g (p-iM(r,t) + aF p Y 



Sal(x) 
6 

JaUx) ° V~ l 2 



a =0 



Tr log (id t --T-n-h- iR^a T D YR 



(7.23) 

aY,=0 
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where the transformation to the flavor rotating frame has been performed (c/. section (7.1)). 
As demonstrated for the baryon number current (6.11) the currents are additive in valence 
and vacuum parts. The valence quark part ivair^) i s obtained by a perturbative expansion 

of the single particle eigenvalues of h + |r • J? + R^ fta T p ^ p R up to linear order in both J? and 
aF . This ends up to 



??valiVc 



6a 1 



+ E 



(val|i?t/3a^ 7 p i?|val) 

{val\T a \fi)(fi\R^a r j p R\val) 



-tt n 



/i^val 



£val 



(7.24) 



aY,=0 



where "val" labels the valence quark orbit. The evaluation of the vacuum part of the currents 
ivacr^) Proceeds in a fashion similar to the computation of the moment of inertia [70] (see 
section 7.1), i.e. we again have to continue to Euclidean space and consider the limit T — > oo. 
However, in the present case also the imaginary part of the Euclidean action contributes. 
Besides the angular velocities J? also the time component of the external gauge field has to 



ia 



4 1 



with af Hermitian. Furthermore, in order to apply the proper time 



be continued as cXq 

regularization scheme, one has to distinguish between real and imaginary parts of the fermion 
determinant. For the real part the regularization is defined in eq (6.25) while for the imaginary 
part a procedure analogous to (6.35) is employed. The part of the fermion determinant which 
does not contain the angular velocity receives contributions proportional to a r (k = 1,..,3) 
from the real part and proportional to from the imaginary part. Taking into account that 
t ■ Qe (7.7) also behaves like a time component of a vector field one finds that the parts 
which are linear in the Euclidean angular velocity have ft E af stemming from the real part 
and fl E ap from the imaginary part. Finally, the total expression for the vacuum part of the 
current is continued back to Minkowski space 



c 



6 



2 8aF v (x) 



Y,^{e,)eric(^)(p]R^fja T p fR\ij) 



Y,fe(^e^A)^n a ,(ij]T a \^^\R^a r p YR\^) 



(7.25) 



aY,=0 



The regulator function /© is defined in eq (7.14). Its appearance in (7.25) is not accidental but 
rather guarantees the proper normalization of the nucleon charge when the quantization rule 
(7.18) is employed. The expressions (7.24,7.25) were obtained in ref. [101] for the evaluation 
of magnetic moments and the axial charge of the nucleon. 



7.2.1 Electromagnetic properties of the nucleon 



In order to extract the electromagnetic current we put a T v (x) = a®' 111 ' = Qa^x) where Q 
is defined above eq (7.22). Obviously Q contains isovector ~ r 3 /2 and isoscalar ~ 1/6 parts. 
Accordingly, the electromagnetic (x) current may be decomposed into isovector jy(x) and 
isoscalar js( x ) contributions. 

Let us first discuss the magnetic moment operator. As in all static soliton models these 
are obtained as the spatial integral of the space components of (x) via fi = (1/2) / d 3 rr X 
Jem ( x )- ^ i s t nen obvious that fi involves the matrix elements of r X f3~f = r X a between 
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Table 7.1: The magnetic moments in nucleon magnetons as functions of the constituent quark 
mass m compared to the experimental data. (Taken from ref. [101]). 



m(MeV) 


350 


400 


500 


600 


700 


800 


expt. 


fJ-s 


0.68 


0.61 


0.54 


0.51 


0.48 


0.46 


0.88 


fJ-v 


2.87 


2.64 


2.29 


1.99 


1.75 


1.56 


4.70 


H 


1.78 


1.63 


1.42 


1.25 


1.12 


1.01 


2.79 




-1.10 


-1.02 


-1.38 


-0.74 


-0.68 


-0.55 


-1.91 



the eigenstates ^ of h (6.43). It is most convenient to consider fi z . For the isovector part of 

the magnetic moment one has t 3 R = D^Ti. For the isoscalar part the explicit dependence 
on R drops out. The quantization rules (7.21) and (7.18) (for isovector and -scalar parts, 
respectively) are used to replace the explicit dependence on the collective coordinates by 
collective operators which act on nucleon states with spin-projection J 3 = 1/2. It is then 
obvious that the magnetic moments operator can be written as fj, = + l 3 [iy . Both n$ 
and fiy may be decomposed into valence and vacuum parts according to (7.24,7.25) ps,v = 
pLgy + p^sv- Explicit evaluation shows that only the first terms on the RHS of eqs (7.24,7.25) 
contribute to fiy 



lTy l = -N C T] vld ^^(val\T 3 (r x a) 3 |val) 

/4 aC = Nc ~f~ IZ(H T 3 (r X a) 3 |//)sgn(e /i )erfc 



A 



(7.26) 



where Afjy = 939MeV denotes the experimental nucleon mass 6 . On the other hand the second 
terms on the RHS of eqs (7.24,7.25) contribute to ps 

vai M M N ^ (vaI|r 3 |//)(//| (r X q) |val) 
= -^c^vaiT,— r 2^ z 

^ ac = -Ncf^ J2M^^ A WT3\v)(v\(rx<*) 3 \p). (7-27) 

in/ 

In all cases use has been made of the fact that the choice of the ^-component for the magnetic 
moment operator projects out the ^-component in matrix elements like (//|t 4 -|z/) after summing 
over the grand spin projection. 

In table 7.1 the numerical results as obtained in ref. [101] for the isovector and -scalar 
moments are displayed for various constituent quark masses m, the only free parameter in 
the baryon sector of the model. Both ps and \iy are seen to decrease with increasing m. 
Furthermore the corresponding values for the proton and nucleon magnetic moments p, p , a = 
(l/2)(/is ± fiy) are compared to the experimental data. Although the isoscalar part of the 
magnetic moment is reasonably well reproduced, the isovector part comes out too small. It 
should be added that in the calculations of ref. [103] a somewhat larger fiy is obtained since 
these authors do not regularize terms which stem from the imaginary part of the action. We 
will comment on a possible solution to the problem of too small a fiy in section 7.2.3. 

In order to obtain the charge radii one needs the second moment of the electromagnetic 
charge density (r 2 ) = / d 3 rr 2 j® m (x). Again a decomposition into isoscalar and -vector parts 



b The magnetic moments are measured in nucleon magnetons 
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Table 7.2: The mean squared charge radii (in fm 2 ) as functions of the constituent quark mass 
m compared to the experimental data. (Taken from ref. [103]). In parenthesis the results for 
a regularized imaginary part are given. 



m(MeV) 


370 


420 


450 


expt. 


(r 2 )s 


0.63 (0.64) 


0.52 (0.50) 


0.48 (0.45) 


0.62 


(r 2 )v 


1.07 


0.89 


0.84 


0.86 


(r 2 ) P 


0.85 (0.86) 


0.71 (0.70) 


0.66 (0.65) 


0.74 


(r 2 )n 


-0.22(-0.22) 


-0.18(-0.19) 


-0.18(-0.19) 


-0.12 



can be carried out 



(r 2 ) = ^(r 2 )s + h{r 2 )v 



(7.28) 



with 



(r 2 )s 



(r 2 )v 



N ° ?7vai(val|r 2 |val) - ^ ^(/i|r 2 | / u)sgn(e At )erfc 



3 

Nc_ 
6a 2 



(7.29) 



??val 

/i^val 



(val|r 2 r|/i) • (//|r|val) 



£val 



1 



+ 2 J2M e ^ e ^ A)(//|r 2 r|i/) • {v\r\fi) 



(7.30) 



These expressions correspond to the case that the imaginary part of the Euclidean action is 
regularized. Noting that (r 2 )y originates from those terms in (7.24,7.25) which are linear in J? 
it should be obvious that only (r 2 )s receives contributions from the imaginary part. Choosing 
not to regularize this part then corresponds to replacing the complementary error function 
in (7.29) by unity. In that case numerical results are also available [103] which are displayed 
in table 7.2. Obviously the question of regularizing the imaginary part only plays a minor 
role. For the range of m displayed in table 7.2 the valence contributions strongly dominate. 
Especially the vacuum contribution to (r 2 )s is only about 1-2%, however it increases with m 
for (r 2 ) s and (r 2 ) v [103]. 

Nucleon Compton scattering provides access to electromagnetic polarizabilities which thus 
may be taken as a measure for the square of dipole transitions between baryon states. Tech- 
nically, the polarizabilities may be extracted from the response of the soliton to external 
electromagnetic fields. For this purpose the action for the rotating soliton (7.1) is expanded 
up to quadratic order in the source = A^Q, i.e. one higher order than for the current 
(7.23). It is sufficient to only consider homogeneous electromagnetic fields in this context: 

= (Ez, ^tijkXjBk)^. When expanding the action the electric (magnetic) polarizabilities 
a (/3) are then read off from the coefficients of E 2 (B 2 ). In the NJL model so far only the 
isoscalar electric polarizability of the nucleon, a/=o = \{ol p + a n ) has been investigated [104]. 
This quantity is sensitive to the large distance behavior of the chiral angle 0(r). Further- 
more, the leading term in a chiral expansion is known [105] to be proportional to g\ (qa 
denotes the axial charge of the nucleon.). The authors of ref. [104] also took into account the 
subleading (in 1/Nc) contributions to a/=o, which unfortunately are subject to ordering ambi- 
guities in the quantization procedure and might be artificial. In any event, these contributions 
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Table 7.3: The axial charge of the nucleoli function of the constituent quark mass m. 



m(MeV) 


350 400 500 600 700 800 


expt. 


9a 


0.80 0.76 0.70 0.66 0.62 0.59 


1.26 



have the desired effect of rendering the experimental value for §a (see section 7.2.3). For 
the constituent quark mass m = 420MeV the isoscalar electric polarizability was found to be 
a i = o = 19 x 10~ 4 fm 3 to which the subleading terms contribute about 3 X 10~ 4 fm 3 . Although 
this result overestimates the experimental value (9.6 ± 1.8 ± 2.2)fm 3 by about a factor of 2 
it is still significantly smaller than in other soliton models which predict §a correctly as e.g. 
the cr-model [105]. It should be added that the calculation of ref. [104] contains the further 
simplification that the 1/Nc subleading terms have been approximated by the leading order 
expression of a gradient expansion. The numerical result a/=o = 19 X 10~ 4 fm 3 may further 
be lowered when one takes account of the quantum character of the pion fields. In that case 
the seagull contribution to a/=o has been claimed to be absent as a consequence of current 
conservation [106]. In the calculation of ref. [104] the seagull terms contributes about 60%, 
though. 



7.2.2 Axial charge of the nucleon 



In order to extract the axial properties of the nucleon we put a^(x) = a J/ (x) a 75y. Eqns 
(7.24,7.25) then provide the axial current j^' a . Due to isospin invariance the axial charge §a 
of the nucleon may be obtained as the matrix element of 2J5' 3 between proton states with spin 
projections +|. Considering the angular velocities S7 as commuting c-numbers one observes 
that there are no contributions to g& from the terms linear in J? (see, however, subsection 
7.2.3). Then one only needs to evaluate the matrix element (p f |Z) 33 |p f) = —1/3 as a 
consequence of (7.21). Further use of /3^ 3 ^ 5 = cr 3 yields [108] 



9A = -^TT \ ?7val(val|cr 3 T 3 |val) - \- ^( / u|cr 3 T 3 | / u)sgn(e At )erfc 




A 



(7.31) 



As gA is obtained to be the matrix element of spatial components of a current and does 
not contain any J? dependence it is obvious that gA completely stems from the real part of 
the fermion determinant and thus necessarily undergoes regularization. It should be noted 
that in the chiral limit (m^ = 0) the RHS of eq (7.31) acquires an additional factor (3/2) 
corresponding to the symmetric zero momentum transfer limit [9]. 

One may as well employ PCAC (2.17) in order to relate gA to the profile function 0(r) 
[9, 109] 



4tt „ „ Q <90 



8tt 



9a = —f: lim i? 3 — _ d - —timl Jrr 3 sin0. (7.32) 



t=r 9 



As a matter of fact these two formulas for evaluating gA provide an excellent check on the 
accuracy of the numerical results which are displayed in table 7.3. The results extracted from 
eqs. (7.31) and (7.32) differ by 1-2% only. Obviously the experimental value is underestimated 
by about 40% for all values of the constituent quark mass. 

The calculation of nucleon axial charges has been extended to those corresponding including 
strange quark currents. Then there are two more axial charges of the nucleon which recently 



76 



received special attention in the context of the "proton spin puzzle" [107]. The first one, g A 
is obtained by assuming a T v (x) = af / (x)^ 5 \g/2 in the definition of the currents (7.22) with A a 
being the Gell-Mann matrices. The second, g A is related to the axial singlet current defined 
via al(x) = a°(x)7 5 /3. 

In order to compute g A and g A the model has to be extended to flavor SU(3). This is 
achieved most straightforwardly by choosing the collective coordinates R(t) in eq (7.72) to 
be 3 X 3 unitary matrices. This is often referred to as the collective approach to include 
strangeness in chiral soliton models. It will be discussed extensively in section 7.5.1. Here we 
will not go into detail but rather remark two essential differences to the two flavor case. First, 
one has to take into account that the SU(3) nucleon wave-function differs from the SU(2) case 
especially in the flavor symmetric case. One consequence is that (p f |Z) 33 |p f) = —7/30 [110] 
i.e. a significant reduction from (7.21). Thus the axial charge for nucleon f3 decay is altered 
when considering an SU(3) collective quantization even without changing the chiral soliton. 
We will refer to this axial charge by g A . Secondly one has to include the flavor symmetry 
breaking arising from the different quark masses. This gives rise to further changes of the 
nucleon wave-function as will be described below . Furthermore the extension to three flavors 
leads to additional terms for the currents not contained in (a three flavor generalization of) 
(7.24,7.25) [112]. These effects have been taken into account in a perturbative scheme with 
the strange current quark mass as the perturbation parameter while non-strange and strange 
constituent quark masses were identified [113]. This scheme differs somewhat from the more 
elaborate one discussed below, however, as the baryon spectrum comes out similar (when 
physical values for the parameters are assumed) there is enough reason to take the results on 
the axial charges seriously. From g A , g A and g A one may finally extract the contributions of 
the individual flavors to the axial charges^ Am, A d and As. Assuming a constituent quark 
mass m = 420MeV the authors of ref. [113] obtain 

Am = 0.64 Ad =-0.24 As = -0.02. (7.33) 

These should be compared with the "experimental" data given by Ellis and Karliner [114] 

Am = 0.81 Ad =-0.44 As = -0.12. (7.34) 

The absolute values for Am and Ad turn out to be somewhat too small as a result of the 
too small g A . When comparing (7.33) with (7.34) one should be somewhat careful since the 
extraction of the latter from semi-leptonic hyperon decays involves the assumption of flavor 
symmetric baryon wave-functions [115]. Within such a framework usually the absolute value 
of As comes out too large. This assumption has been dropped in the collective approach. 
From studies in Skyrme type models it is well known that nevertheless the branching ratios 
for the semileptonic hyperon decays are well described in the collective approach together with 
a small |As| [116]. As most chiral soliton models [96, 44] the NJL model reasonably explains 
the smallness of the quark contribution to the proton spin Am + Ad + As 0.38. 

7.2.3 Remarks on 1/Nc corrections 

In all preceding considerations we have considered the angular velocities J? as commuting 
c-numbers. For the computation of nucleon observables J? is replaced by the spin operator 

c For infinitely large strange quark masses one has g A = g A , i.e. the two flavor limit [111]. 
d E.g. g\ = Am - Ad. 
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J according to (7.18). However, in non-leading contributions ordering ambiguities exist and 
different results may be obtained because of non-vanishing commutators 



or 



[Ji, D a b 



£ e lbm D am . 



(7.35) 



m = l 



This substitution leads to non-vanishing matrix elements between nucleon states 



e lmn (N\D tl tt m D rn \N) 



a 



(7.36) 



The prefactor on the right hand side is ambiguous, in particular, it is vanishing in the ordering 
assumed in section 7.1. Eqn (7.36) then yields additional contributions to \iy and §a which we 
denote by fiy^ and g^\ respectively. These corrections are subleading in an 1/Nc expansion 
because the moment of inertia a 2 is of order Nc- Since the starting point for the computation 
of the baryon observables is completely classical the ordering between the angular velocities 
and the rotations matrices in the prescription (7.36) is ambiguous. Adopting a "natural" 
ordering [117], which is suggested by the perturbation expansion of the quark wave function 
in the flavor rotating frame 



R(t) 



(h\t ■ S7\v) 



R(t) 



(fi\r ■ J/a 2 \v) 



(7.37) 



the numerical results for fiy^ and g A ^ are about 30% of the leading order [118]. Thus it 
is suggestive but not obvious that the series has already converged at this order. Similar 
calculations have been performed [113] for the three flavor axial couplings discussed in the 
preceding subsection. For these quantities the quantization prescription (7.35) leads to 1/Nc 
corrections of the same order. 

At the time when these corrections were first observed they were highly welcome since they 
helped to solve (among others) the problem of the too small g^ (see table 7.3). Nevertheless 
the applicability of (7.35) from the very beginning remains doubtful, especially since a special 
ordering of the collective coordinates and operators has to be adopted. In order to obtain 
some restrictions on possible orderings the compatibility of the associated 1/Nc corrections 
to gA with symmetries of the underlying model have been examined. 

As the 1/Nc corrections appear to be a special feature of the semiclassical quantization 
procedure the self-consistent soliton profile remains uneffected. Thus g^ evaluated with eqn 
(7.32) assumes the leading order value. The correction g A ^ may then be interpreted as to 
violate PCAC by 30% which, of course, represents an undesired breaking of one of the funda- 
mental symmetries of the model. It has been shown that this problem can be solved by simply 
adopting PCAC as the equation of motion instead of the stationary condition (6.40) [109]. 
This then leads to an extended soliton profile function causing the moment of inertia and the 
radii to increase. Unfortunately, until now this equation of motion has not been derived from 
an action principle. 

Furthermore, requiring the proper behavior of the symmetry currents under particle con- 
jugation may provide additional constraints on orderings of the collective coordinates. This 
has recently been done in the framework of the chiral quark model [119], which at least for 
the valence quark part of the action is similar to the NJL soliton model. It has been shown 
[120] that the axial current transforms under this symmetry properly only when an ordering 
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is chosen such that the correction g A actually is zero. For example the Hermitian ordering 
prescription 

D ab n 3 ^- (D ab Jj + J 3 D ab ) . (7.38) 

leads to symmetry currents which transform properly under particle conjugation while there 
are no corrections to qa and \iy at the subleading order in the 1/Nc expansion. In ref. [112] an 
expression for g^ A has been presented, which involves the absolute value l/(|e M — e^l) rather 

than l/(e M — e v ) as in the chiral quark model. Hence this result for g^ A complies with particle 
conjugation. It should, however, be remarked that the alternative quantization prescription 
(7.38) leads to vanishing 1/Nc corrections in the NJL model as well. 



7.3 Meson fluctuations off" the chiral soliton 



In order to explore baryon properties like e.g. tt — N scattering 6 , electro-magnetic po- 
larizabilities (beyond the isoscalar electric channel) [105] or for the investigation of quantum 
fluctuations it is mandatory to go beyond the zero-mode quantization described in section 
(7.1). For applications like these time-dependent fluctuations off the soliton have to be in- 
cluded and quantized canonically. This formalism is also of relevance for the discussion of 
quantum corrections to the soliton mass [25] (section 7.4) as well as the description of hyper- 
ons in the bound state approach [29] (compare subsection 7.5.2). 

In order to formally introduce the fluctuations we consider the ansatz [97] 

M = ^/(£R/Co. (7.39) 

Here £ = exp [^t ■ r Q(r)j denotes the hedgehog soliton configuration while the space-time 
dependent fluctuations rj a (x) are given by 

Z f (x)=exp(ij2ri a (x)\ a /2y (7.40) 

The main task now is to expand the action up to quadratic order in the fluctuations rj a . In 
order to avoid problems with stability the meson fields have been constrained to the chiral 
circle. Although the parametrization (7.39) deviates from the commonly adopted unitary 

gauge M = ^j- (S)£r, ^j- = £r it has proven to be convenient when matrix elements between 
eigenstates of the static Hamiltonian (6.42) are computed. This can be observed easily by 
considering the Euclidean Dirac operator pE 

iPp E = -d r -CL-p- T/3 (iflftifPR + i) Tt (7.41) 

wherein r = ix is the Euclidean time. The unitary matrix 

T = Z P l + 4pr (7-42) 



3 For a review on 7r — N scattering in soliton models see ref. [11]. 



79 



contains all the information on the chiral soliton. Thus, whenever a matrix element involving 
the fluctuations has to be evaluated it can be simplified by chirally transforming the states 
between which the operator is sandwiched. 

As the first step towards expanding the action we write the Euclidean Dirac operator as 



ifipE = -d T -h = -d T - [h {0) + h {1) + h {2) + ■■■) (7.43) 

wherein the subscript labels the power of the meson fluctuations. The dots indicate powers of 
the fluctuations larger than two. Upon noting that as given by (6.42) is time independent 
one obtains for the argument of the real part of the fermion determinant Ar (6.22) 

plpE = -d 2 T + h\ 0) - [d T , h (1) ] + {h (lh h (0) } - [d T , h (2) ] + {h {2) , h (0) } + h 2 (1) + • • (7.44) 

When the valence quark contribution as well as the terms originating from the mesonic part of 
the action are included it is obvious that the zeroth-order (in the fluctuations) just renders the 
static energy functional while the expression linear in rj a (x) vanishes subject to the equation 
of motion for the static soliton. Therefore only the second order expression is of interest 
for the current discussion. The corresponding contribution from the real part of the fermion 
determinant reads 

~j p oo ^ /* oo /* s 

4? = 2 Tr y i/A2 dsI ^ V)l + h h) - 4 Tr J ds I ds'k Q {3 - s 1 ) 



x ([d T , h {1) ]K (s')[d T , h {1) ] + {h {1) , h {0) }K (s'){h { i), h {0 )}) ■ (7.45) 



The introduction of the zeroth-order heat kernel K (s) = exp ys(d 2 — ^(o))J nas turned out to 
be useful for this computation. For the imaginary part of the fermion determinant, Ai } the 
cut-off A is introduced by the substitution 

((<9 T - h {0) - h {1) )(d T + h {0) + fy^)) 

— > - ^ /a2 ds exp (s(d T - h (0) - h {1) )(d T + h {0) + (7.46) 

which is legitimate since the argument is negative definite, i.e. it converges for large momenta. 
The resulting imaginary part starts off at second order in the fluctuations 

ds ds'K (s - s')d T h {1) K (s')h {0) h {1) + ■■-. (7.47) 

To perform the temporal part of the functional trace in Euclidean space the meson fluctu- 
ations are Fourier transformed 

f+°° dio 

Va (r } -it) = / —fj a (r, iu)e-™ (7.48) 

J -co ZTT 

This transformation directly transfers to the Hamiltonians: 

/ + OO (fjjj „ 
— h {1) (r,iuj)e~ tuJT and 
-oo Z7T 

/+°° cluj /*+°° duj' ~ 
— / — h {2) (r,zu,zu')e-^^ (7.49) 
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since the chiral transformation (7.42) is time independent. In order to project out the vacuum 
part of the fermion determinant, *4 vac , one considers the limit of infinitely large Euclidean 
times. Noting that [70] 



(t\K (s)\t') = -^=exp(-sh 2 {0) )exp y — j (7.50) 

the temporal part of the trace then amounts to carrying out Gaussian integrals involving the 
Fourier frequency u. The spatial part of the trace as well as the traces over Dirac and flavor 
indices are evaluated using the eigenstates of the static one-particle Hamiltonian (6.43). 
Finally, the frequency u has to be continued back to Minkowski space in order to obtain 
physically relevant expressions. The second order (in meson fluctuations) contribution to the 
vacuum part of the fermion determinant finally results in [97] 

I Jl A 2 x/AlTS ,, J-oo Z7T 



2 Jl/A 2 V47TS 

N c f 00 , l~s~^ f +co dto 



+ 4 ii/A^VI^ l-oo (7-51) 
X I + [u? - (e^ + e u ) 2 ]Ro(s;u, e^, e v ) - icoe^R^s; to, e^, e v ) > . 

The information on the orderings of the operators in eqns. (7.45,7.47) is contained in the 
Feynman parameter integrals 



R,(s; u , £( „ £ „) = /' .VI, exp (-.,[( 1 - + * - s(l - *M) (7.52) 



which represent moments of the quark loop in the presence of the soliton. 

Besides the polarized vacuum configuration also the explicit occupation of the valence quark 
level contributes to the action as long as the associated energy eigenvalue e va i is positive. Since 
no regularization is involved the computation is completely performed in Minkowski space. 
Treating the meson fluctuations as time-dependent perturbations the associated first order 
change 6^ veL \ of the valence quark wave-function ^vai is obtained to be 

6V ysi (r,t) = (id t - fyo)^))" 1 h {1) (r,t)V ysi (r,t). (7.53) 
The corresponding contribution to the second order part of the action reads [97] 

(2) f + °° did ( 

*4vai = -VvdNc / — f (val|^ (2) (r,cc;, -^)|val) 

+ j2 ( val l^(i)( r ^)l/ i )(/ i l^(i)( r ^~ cu )l val ) \ ( 7 54 ) 



Here ?y va i = 0, 1 again denotes the occupation number of the valence quark and anti-quark 
states. 

The complete second order contribution to the action is then given by 

^ (2) = 4 a ) c + 4 a ) 1 + ^ ) - (7-55) 
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The mesonic contribution, is obtained by substituting the ansatz (7.39) into the expres- 

sion for the mesonic part of the action (3.25) 

+ + — )( cos ° +^)E^H^(-^)} (7-56) 

wherein use has been made of the relation G\ = m m/m'^f% (4.41). Obviously A^ contains 
terms of odd powers in u. These correspond to the imaginary part in Euclidean space and 
have the important property of removing the degeneracy between solutions with ±cu [29]. In 
the Skyrme model these terms originate from the Wess-Zumino action [52] which is identical 
to the leading order term of the gradient expansion of the imaginary part [31], cf. section 4.1. 

In general the second order contribution to the action can be written as a functional of the 
fluctuations r\ a 

+ j d\ ^(r)r} a (r,-u)r} b (r,u)}. (7.57) 

The local and bilocal kernels have to be computed as mode sums involving the eigen- 

states and -functions of (6.43). 

Using the Wigner-Eckart theorem it can be shown that due to the symmetry of the soliton 
under grand spin and parity transformations, the fluctuations decouple with respect to their 
grand spin and parity quantum numbers. I.e. the kernels are diagonal in these quantum 

numbers. This property turns out to be helpful when investigating quantum corrections to 
the soliton mass [67] and hyperons in the bound state approach [99] to the NJL soliton. 



7.4 Quantum corrections to the soliton mass 



The action functional (7.57) of the meson fluctuations in the background field of the NJL 
chiral soliton has originally been derived to make possible a description of hyperons in the 
NJL model within the so-called bound state approach (cf. subsection 7.5.2). Very recently it 
has been shown in the two flavor reduction that the action (7.57) also allows one to estimate 
the quantum correction to the soliton mass [67]. Then the bilocal kernel $( 2 ) only depends on 
cu 2 , i.e. the contributions to the action corresponding to odd powers in the frequency vanish. 
This has the consequence that the solutions to the Bethe-Salpeter equation-^ 

J d 3 r' ^{r,r',u)i\ h {r',u) + $$(r)»/ 6 (r,u;) = (7.58) 

appear in pairs ±iu 8 -. Denoting the corresponding wave-functions which solve (7.58) by 
f}^\r^Ui) the fluctuating field may be decomposed as 

r,(r,t) = "7= {a#(^'K^ + 4f7«(r,^) e -^4 . (7.59) 



method for numerically solving this equation is provided in ref.[121]. 
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After canonical quantization (i.e. assuming the commutation relations 
Hamiltonian associated with the fluctuations is that of an harmonic oscillator 

f 



ajcii + 



6; A the 



(7.60) 



This result (7.60) is quite non-trivial since the action (7.57) involves all orders of time deriva- 
tives rather than terminating at quadratic order as e.g. in the Skyrme model. This new 
feature leads to an involved energy functional (in coordinate space) as well as a complicated 
orthonormalization condition for the meson fluctuations. E.g. the normalization condition for 
a solution fj(r,Ui) to the Bethe-Salpeter equation reads 



^3 [J3 >~{i)( ,d*f(T,T',U>) 

a r a r rj K a >(r,uji)- 



du 2 



Vb 



(i)t 



1. 



(7-61) 



Unfortunately, the orthogonality condition for solutions with different frequencies cannot be 
presented in such a closed form. 

It has also formally been shown that the vacuum contribution to the energy (^C;^/2) is 
obtained to be of the form (7.60) whenever the background field is static and the eigenvalues 
of the Bethe-Salpeter equation appear in pairs ±iu 8 -. 

The form (7.60) implies the existence of an operator acting in the space spanned by the 
eigenstates of (7.58) 



H 2 = H 2 + V 



(7.62) 



with eigenvalues ±iu 8 -. This operator can be expressed as the sum of the corresponding one 
in the absence of the soliton, Hq } and a "perturbation", V, which depends on the soliton. In 
ref. [25] it has been shown for the Skyrme model that a finite (renormalized) energy correction, 
AE } is obtained from these operators via 



AE 



I, 



Tr [H - H 



(7.63) 



which actually represents the generalization of the quantum corrections to the kink mass [122] 
to 3 + I dimensions. Although the NJL model is quite different from the Skyrme model the 
trace may similarly be computed by expressing it in terms of the corresponding eigenvalues 



AE 



I 



E 



i 



3 



(0) 



)) 



(7.64) 



3 + 6 



(0) 



(0) 



,(°) 



and denote the eigenfrequencies and -wave-functions to the Bethe- 

Salpeter equation (7.58) in the absence of the soliton. The overlap in eq (7.64) has occurred 
because the operators H and H act in distinct Hilbert spaces. A reasonable definition of 
these overlaps is gained by including the metric 



Q$ab 



du 2 



M ab (r,r\uj t ) 



3 

d 3 x ^2 

c=l 



be 
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into the wave-function 



<f> a (r,u,i) = J d 3 xJ2(VM) Ca (x,r,u,i) U^^i)- (7.66) 

c=l 

The modificatied wave-function cj) turns out to be independent of the parametrization (7.39) 
[123]. Finally the relevant overlap matrix element is given by 

(77(r,^)|^ (0) (^S (0) )) : = / d " r <f>( r >"i) ■ <£ (0) (^ 0) ) (7-67) 

where cf)^ denotes the analogue of q) in the absence of the soliton. 

From eq (7.64) it is intuitively clear that the zero modes (u>i = 0) provide the major 
contribution to AE because no counterpart exists in the absence of the soliton. The zero modes 
arise because the soliton breaks the rotational 3 and translational invariance. In that sense the 
zero modes may be regarded as Goldstone bosons. The associated wave-functions, ff zia (v), 
correspond to a wave-function in the P-wave channel for the rotational zero mode while the 
translational zero mode contains S- and D-wave parts. According to these structures ansdtze 
can be made for the wave functions in the zero mode channels. Then the Bethe-Salpeter 
equation reduces to (coupled) homogeneous integral equations for purely radial functions. As 
this procedure is quite technical we refer the reader to ref. [67] for details on this calculation 
as well as on the explicit construction of the modified wave function cj). As already noted the 
contribution of the zero modes to the energy correction 



^E-f <^,n>)l^ (0 V,S (0) )> 2 (7-68) 

3 



is negative definite. Hence it automatically leads to the desired result of reducing the total 
energy. 

Before presenting the numerical results one remark concerning the non-confining character 
of the NJL model has to be made. Obviously the transformation from Euclidean to Minkowski 
space is only well defined as long as the exponent in the Feynman parameter integrals (7.12) 
vanishes for large s along the path connecting these two spaces. In the absence of the soliton 
this is only the case for u < 2m. Beyond this threshold the regularization functions develop 
imaginary parts which measure the decay of the meson fluctuations into quark-antiquark 
pairs. Consequently the quantum correction to the soliton mass which stems from the zero 
modes is estimated by the truncated sum 

AP, m . = --| £ uf\(f, z . m Xr)\^\r^))\ 2 . (7.69) 

16 (0 , 

In order to judge the quality of this truncation it is useful to define a sum of overlaps 

$= £ \(v,. m Mv {0) (r,uf))\ 2 (7.70) 

Jp<2m 

which should approach unity if the model were insensible to the truncation. 



3 Due to the grand spin symmetry of the hedgehog ansatz coordinate space rotations and isospin transfor- 
mations are equivalent. 
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Table 7.4: The quantum corrections to the soliton mass due to the rotational zero mode [67]. 



m(MeV) 


m v = 
400 500 600 


m v = 135MeV 
400 500 600 


s 

A£ z . m .(MeV) 


0.88 0.94 0.95 
-201 -274 -290 


0.83 0.89 0.91 
-244 -297 -323 



Table 7.5: The quantum corrections to the soliton mass due to the translational zero mode. 
The contributions stemming from the S(l = 0)- and D{1 = 2)-waves are disentangled [67]. 



m(MeV) 


m v = 
400 500 600 


m v = 135MeV 
400 500 600 


s 

A£ /=0 (MeV) 
A£ /=2 (MeV) 
AE z . m .(MeV) 


0.32 0.41 0.52 
-18 -22 -32 
-127 -140 -207 
-145 -162 -239 


0.28 0.37 0.46 
-12 -22 -32 
-82 -128 -187 
-94 -150 -218 



The numerical results obtained in ref. [67] for the contributions of the rotational and trans- 
lational zero modes to the energy correction are displayed in tables 7.4 and 7.5, respectively. 
Also shown are the corresponding sums of overlaps, S. The chiral limit (m^ = 0) as well as 
the physical case (m^ = 135MeV) are considered. Note that the number of states which lie 
below to = 2m decreases as the pion mass increases. In case of the rotational zero mode S 
approaches unity close enough to consider AE —250 ~ 300MeV as a reliable estimate of 
the quantum correction to the energy. The situation, however, is not as good for the trans- 
lational zero mode. The corresponding wave-function is more strongly peaked at r = than 
in case of the rotational zero mode. Therefore the Fourier transform acquires contributions 
from modes carrying higher frequencies which lie above 2m. Hence S hardly reaches 0.5. The 
corresponding AE can then only be considered as a lower bound. As in the Skyrme model it 
turns out that the D-wave contributions strongly dominate the S'-wave part. 

In ref. [67] it has furthermore verified that the scattering modes provide a negligible con- 
tribution of a few MeV. 

Taking the numbers displayed in tables 7.4 and 7.5 seriously (improving the model such 
as to overcome the problem of non-confinement would in any event also alter the predictions 
on the classical mass, moment of inertia, etc.) one obtains a mass formula for baryons 

M = E cl + AE+ J{ i\ l) (7.71) 

with AE estimated by the sum of the rotational and translational zero mode contributions. 
The last term in (7.71) arises from the semi-classical cranking procedure discussed in section 
7.1. In general there are also quantum correction to this rotational term. However, these are 
of the order 0(Nq ) and hence they are omitted. The numerical results for the baryon masses 
are given in table 7.6. Reasonable agreement with the experimental data for the masses of the 
nucleon (939MeV) and the A-resonance (1232MeV) are only obtained for constituent quark 
masses tyi ~ 400MeV. Unfortunately S then is as low as 0.4 for the translational zero mode. 

As a conclusive statement of this section we would like to mention that the mass of the 
nucleon is significantly lower than the soliton mass due to quantum corrections. In general 
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Table 7.6: The predictions for the masses of the nucleoli (iV) and A-resonance. The empirical 
data are 939MeV and 1232MeV, respectively. (Results taken from [67].) 



m(MeV) 


m v = 
400 500 600 


m v = 135MeV 
400 500 600 


^cl(MeV) 

AE(MeV) 

a 2 (l/GeV) 

Mjv(MeV) 

M A (MeV) 


1212 1193 1166 
-302 -436 -525 
6.26 4.73 3.87 
970 836 738 
1210 1153 1126 


1250 1221 1193 
-338 -448 -538 
5.80 4.17 3.43 
976 863 764 
1236 1223 1201 



these corrections reduce the masses of baryons by a few hundred MeV as compared to the 
cranking result (7.20). The effect of the quantum corrections on observables other than the 
masses has not yet been studied in the NJL model. 

7.5 Hyperons 

In this section we will describe the treatment of strange degrees of freedom within the NJL 
model of pseudoscalar fields. The main goal is to find a description of the hyperon spectrum in 
this model. As has already been pointed out, the projection of the soliton onto states with good 
spin and flavor quantum numbers can only be performed approximately. This is already the 
case for the (symmetric) two flavor model. For strange degrees of freedom the situation is even 
worse due to the presence of symmetry breaking. Calculations in Skyrme type models provide 
two different approaches which are frequently viewed as opposite limits of symmetry breaking. 
Both treatments represent approximations to the exact time dependent solution which is yet 
un-known. The first one represents a generalization of the zero mode quantization of SU(2) 
and requires the introduction of collective coordinates describing rotations in the whole SU(3) 
flavor space [26, 27]. We will therefore refer to this treatment as the collective approach. As 
these collective coordinates describe large amplitude fluctuations the restoring force is assumed 
to be small. Stated otherwise, the flavor symmetry breaking is considered to be small and 
an expansion in the parameters measuring symmetry breaking is performed. The collective 
treatment has undergone considerable improvement when Yabu and Ando [28] observed that 
the resulting collective Hamiltonian including symmetry breaking terms can be diagonalized 
exactly by numerical methods. Later on this treatment was seen to represent admixtures of 
states from higher dimensional SU(3) representations to the basic octet and decuplet states 
[111]. In the complementary treatment, which was initiated by Callan and Klebanov [29] 
and to which we will refer as the bound state approach, the starting point is to consider 
symmetry breaking large, although the treatment yields the correct results in the symmetric 
limit as well. Here only small amplitude fluctuations are allowed (c/. section 7.3). These 
fluctuations can be quantized canonically. The bound state approach heavily relies on the 
special feature that a solution to the Bethe-Salpeter equation (7.96) emerges in the zero- 
mode channel when the symmetry breaking is switched on. The eigenfrequency of this mode 
is different from zero but also significantly lower than the kaon mass, i.e. it represents a bound 
state. Additionally the bound state wave-function is well localized. I.e. the bound state is the 
"would-be" Goldstone boson of the strange flavor transformations in the soliton background. 
By construction the occupation number of this bound state is identical to the strangeness of 
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the baryon under consideration. Finally the real zero modes corresponding to spin and isospin 
are treated within the collective approach yielding the hyperfine splitting and thus removing 
the degeneracy of baryons with identical strangeness as e.g. the E and the A [29]. 

Here we will discuss both approaches in the framework of the NJL model and critically 
compare the results. 



7.5.1 Collective rotational approach 



In the collective approach symmetry breaking is considered to be small and consequently 
strange degrees of freedom are introduced as if they were zero modes. On top of this, effects 
due to symmetry breaking are treated by expanding the fermion determinant in terms of the 
difference of the constituent quark masses. It should be stressed that for the mesonic part of 
the action A m (3.25) no expansion is performed. 

According to the above discussion, collective coordinates for rotations are defined in the 
whole flavor space in order to approximate the time dependent solution. As the effects due to 
flavor symmetry breaking have to be taken into account we have to go beyond the ansatz for 
the two flavor case (7.1). We consider [62] 

M(r,t) = R(t)£ (r)ltf(t)(X)R(t)£ (r)ltf(t) R(t) £ SU(3). (7.72) 

Here £ = exp (it ■ r 0(r)/2) denotes the static soliton configuration. Obviously, only the 
pseudoscalar fields rotate in flavor space while the scalar fields are kept at their vacuum 
expectation values. For the ongoing exploration it is helpful to define eight angular velocities 
fi a , (a = 1, .., 8) which measure the time dependence of the flavor rotation' 1 

l -J2Ktt a = Rht)R(t). (7.73) 
a=i 

These are the extensions of the previously introduced angular velocities (7.4) to the SU(3) 
flavor group. Again it is convenient to transform to the flavor rotating system in order to 
evaluate the fermion determinant: q = Rq' . This allows to eliminate the "outer" rotations in 
(7.72) at the expense of an induced rotational part 

1 8 

Kot = -J2 x - n - ( 7 - 74 ) 

(1=1 

In the rotating frame the Dirac operator acquires the form 

if3p' = id t - /i(o) - hot - h S B- ( 7 -75) 

Here represents the static one-particle Hamiltonian defined in eq (6.42). The introduction 
of the chiral transformation T (7.42) allows one to easily display the symmetry breaking part 
in the Dirac operator 



h SB = T/3(^i?t(E)i?-(E)jTt 

-TP (^DsiXi + AkX + (A* - l)\ s ) Tt. (7.76) 



m — m. 



\i=l a=4 



l The dot indicates the derivative with respect to the time coordinate. 
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We have also indicated the SU(2) invariant pieces. Furthermore use has been made of the 
adjoint representation for the flavor rotations D a b = | (^X a RXf,R^ . 

The main task consists of expanding the fermion determinant Af in terms of h rot and hsB- 
In order to regularize Af a continuation to Euclidean space is required. In this context it is 
important to take into account that h rot corresponds to the time component of an induced 
vector field. Hence h rot has to be considered as an anti-Hermitian quantity. The real part of 
Af therefore contributes terms of even power in J? to the effective action. Since the expansion 
is constrained to the second order in h rot + hss we consider 

pip' E = -d 2 + h 2 0) + {h (oh h SB } + h 2 SB + [h (oh h rot ] - h 2 rot . (7.77) 

In general there could also be contributions from the commutator [d T} hsB] } however, these 
will contribute to the moments of inertia at second order in symmetry breaking and have to be 
discarded for consistency. Since the imaginary part of the fermion determinant, Ai } contains 
only terms of odd powers in the angular velocity, it will receive contributions from the mixed 
terms of the form h rot hsB- Noting that the action is only expanded up to second order in the 
angular velocity, Ai may be obtained via 

M = ^Trlog^)-!^ 

= ^Tr | [(d T - /j(o) - 1i S b) (~<9t - h(Q) - hss)] {h rot , fyo) + ^sb}} + • • • .(7.78) 

The time-dependence of hss actually yields a non-vanishing commutator [d T ,hsB] which, 
however, will not contribute to Ai below (9(0^) and may therefore be discarded for the 
further calculation. Although the imaginary part is finite the ongoing evaluation of Ai can be 
made consistent with the proper-time regularization of the real part. This is achieved by the 
replacement (7.46). 

Now the functional trace can be evaluated. The temporal part is performed by introducing 
eigenstates 8 exp(zcu n r) of d T which satisfy anti-periodic boundary conditions in the Euclidean 
time interval T. The spatial part as well as the traces over Dirac and flavor indices are 
evaluated using eigenstates of the static one-particle Hamiltonian (6.43). 

The leading term in the expansion of the fermion determinant is the vacuum contribution 
to the soliton mass and stems from the real part (c/. chapter 6). The imaginary part starts 
off with an expression related to the anti-commutator {h rot} /i(o)} 



^ Tl * //A2 dS 6XP f S ^ T ~ ^ (0) ) ^ T + ^ (0) ^ ^ r °*' ^ (0) ^ 

-r^sign^erfc ( ^ ((i\h rot \(i). (7.79) 



l/A 2 

N Cr 



The sum over the eigenstates \fj) of obviously projects out the grand spin zero part of h rot 
which is related to the eighth component of the angular velocity 



^Tfig^Tsign^erfc 



A 



TB vac tt 8 (7.80) 



(2n + 1)/T are the Matsubara frequencies, cf. chapter 6. 
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wherein B vac denotes the contribution to the baryon number originating from the polarization 
of the Dirac sea. The remainder of the imaginary part is linear in both, the angular velocity 
and the constituent quark mass difference m s — m and we denote it by (see appendix B for 
the definition of A a &) 

^W«. (7.81) 

The feature, which was already observed for the imaginary part, that the first order in the 
expansion of the real part in terms of + hsB only allows for grand spin symmetric expres- 
sions also holds for the real part. Taking into account that Ar is even in the angular velocities 
there remains only one grand spin symmetric term involving D 88 

since the deviation of the D-matrix from unity is considered to be the dynamically relevant 
quantity (7.76). Furthermore the real part of the fermion determinant contains terms which 
are quadratic in either the angular velocity or the symmetry breaking 

^TO^OA, l -TT2 c D 8a D 8b . (7.83) 

Again we refer to ref. [62] for the actual calculation of the quantities 0™ c and T™ c (see also 
section 7.1 for definitions and appendix D for the explicit expressions). Needless to mention 
that the resulting expression for the action has to be continued back to Minkowski space in 
order to extract the collective Lagrangian. 

Up to this point we have considered the limit of large Euclidean times which has projected 
out the vacuum contribution to the action. Additionally the coefficients of the collective 
quantities fi a and D 8a receive contributions from the explicit occupation of the valence quark 
level. These are obtained by consideration of the extented Dirac equation 

h( ) + hot + hss) ^val = Cval^val (7.84) 

in stationary perturbation theory. Since the valence quark part of the action is not regularized 
there is no need to continue forth and back to Euclidean space. Thus the relevant calculations 
are completely performed in Minkowski space. The resulting expressions are displayed in 
appendix D. Finally there remains the contribution from the meson part of action A m - This 
can be evaluated straightforwardly by substituting the ansatz (7.72) into eq (3.25) yielding 

\A m = "^7 m (l " A*) - ^T^DsiDsi - ^ D saD 8a - \rf(l - D 88 D 88 ) (7.85) 

with the coefficients containing the direct information on the current quark masses m° and 
m° s . Again we refer to appendix D for their explicit forms. 

We have now collected almost all ingredients for the collective Lagrangian. We still need 
to discuss meson field components which vanish classically but are induced by the collective 
rotation into strange direction. Parametrizing the corresponding kaon fluctuation by (c/. the 
previous section 7.3) 
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the imaginary part of the Euclidean action provides couplings which are linear in both, the 
kaon field K as well as the angular velocity fi„ (a = 4, ..,7). Expanding the real part of 
the action in terms of K only contains even powers of the kaon field. Hence a non-trivial 
solution for K(r) may exist. A suitable ansatz which carries the appropriate flavor and parity 
quantum numbers reads [124] 

A» = r.rmr) (£:;£). (7.87) 

In general W(r) is a complex radial function. The framework of the collective approach 
obviously requires to expand the action up to quadratic order in fi a . This procedure adds 

jfr[W] E ^ a (7.88) 

Z a=4 

to the action. The equation of motion obtained from the variation Sfij/ 6W{r) = for the 
radial function W contains a source term stemming from the linear coupling to the collective 
rotation discussed aboved. This source term is completely given by the static (hedgehog) 
field. Hence the solution requires W ^ 0. It turns out that only the real part of W is 
excited. In the NJL model up to now the excited fields have not been treated exactly but 
rather in the gradient expansion. Then the imaginary part of the action is approximated 
by the Wess-Zumino term. The real part of the action is approximated by the non-linear 
cr-model allowing, however, the kaon decay constant fx to be different from /„.. For fx the 
prediction of the NJL model is used (c/. table 4.1). Whenever the chiral angle shows up 
the self-consistent solution is substituted. The reader may consult ref. [62] for details on this 
calculation. The important result is that the induced part of the strange moment of inertia 
fix is not negligible. It should also be stressed that there are no analogous excited fields for 
the moment of inertia associated with rotations in iso-space since the imaginary part of the 
action vanishes identically in a two flavor model as long as only pion fields are included- 7 . 

The derivation of the collective three flavor Lagrangian is now terminated and its final 
form may be discussed. Isospin and rotational invariance provide relations between certain 
components of the moment of inertia 

a 2 := O n = 22 = 33 , /3 2 := 44 = 55 = 66 = 77 and 88 = (7.89) 

while all other components vanish. The last equation in (7.89) stems from the vanishing 
commutator [/i( ),A 8 ]. Analogous relations hold for A a b and T a b. 

In the next step the collective Hamiltonian associated with the Lagrangian (D.13) has to 
be constructed. This is achieved by generalizing the treatment of section (7.1) to three flavors. 
This derivation is also described in appendix D; it results in the mass formula for baryon B 
carrying spin J 

E B = E» + I(i - j z )J(J+D - ± + ±*ss. (7.90) 

The moments of inertias a 2 and f3 2 may be found in appendix D as well. The quantity ess 
denotes the eigenvalue of the SU(3) operator (D.18) which contains the information on the 

J An iso-singlet 77 field, however, may yield a non- vanishing but small contribution to the non-strange 
moment of inertia. 
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Table 7.7: The mass differences of the low-lying | and | baryons with respect to the nucleon. 
We compare the predictions of the collective approach to the NJL model with the experimental 
data. The up-quark constituent mass m is chosen such that the A-nucleon mass difference is 
reproduced correctly, see text. The last column refers to the case that the symmetry breaker 
7 is scaled by (fi* pt '/ff[ ed ') 2 All data (from ref. [99]) are in MeV. 





A fixed 


fx fixed 


Expt. 


Xcorr. 

Jx 


A 


105 


109 


177 


175 


E 


148 


151 


254 


248 




236 


243 


379 


396 


A 


293 


293 


293 


291 


E* 


387 


391 


446 


449 




482 


489 


591 


608 





576 


586 


733 


765 



spin and strangeness of baryon B. As discussed in appendix D the eigenstates of the collective 
Hamiltonian leading to (7.90) are constrained to carry half-integer spin, i.e. they are fermions. 

We are now enabled to discuss the numerical results on the hyperon mass spectrum. As 
already explained in section 4.3 the NJL model underestimates the ratio fx I 'fix by about 15% 
for non-strange constituent quark masses of the order m 400MeV. This shortcoming of the 
model is also expected to show up in the baryon sector. In order to demonstrate that the ratio 
fi</fir rather than the absolute value for fx is the ingredient to which the results are sensitive 
we compare calculations for two different sets of parameters. First, /„. = 93MeV is kept at its 
empirical value and m = 407MeV is chosen such as to reproduce the experimental A-nucleon 
mass difference (293MeV). For this set fx = 99.8MeV is predicted, i.e. fx/fir = 1-07. The 
second set of parameters is fixed such as to correctly give fx = 114MeV. Again m = 433MeV is 
determined by demanding the experimental A-nucleon mass difference. Then /„. = 104.9MeV 
is increased considerably, however, the ratio fx/fir = 1-09 remains almost unaltered. The 
corresponding numerical results for the baryon mass difference are displayed in table 7.7. 
Obviously the change in the mass differences for the two sets of parameters is not larger 
than the change for the ratio fx/fir, *- e - 2%. As expected the SU(3) symmetry breaking in 
the baryon mass differences is considerable underestimated. It has been demonstrated that 
this effect is strongly correlated to problem of incorrectly predicting fx- The leading order 
terms in a gradient expansion to the dominating symmetry breaking parameter 7 are given 
by (7.82,7.83) 

- cos©) + ^^cosO (O' 2 + + . . . j 

cos©), (7.91) 

where the latter relation is verified numerically. From Skyrme model calculations it is well 
known [124] that this scaling of 7 with fx represents the most important dependence on the 
kaon decay constant. In order to estimate the effects of this scaling in the framework of the NJL 
soliton the corresponding quantities (7,^,5,8) m the collective Lagrangian (D.13) are scaled 
by (fx Pt '/fx ed ) 2 ~ (114/100) 2 and the hyperon spectrum is re-evaluated [62]. Again m is 
adjusted to reproduce Afjy — Ma- The results for the mass differences, which are also shown 



vac plus m 4 f j3 J / 2 r2 

7grad. exp. = 3 J d T \ \™ K fx 

4 r 
~ ^xfx j d 3 r (1 - 
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in the last column of table 7.7, agree excellently well with the experimental data. Of course, 
it should be stressed that this estimate has the purpose of relating the too low prediction 
for the SU(3) symmetry breaking in the baryon sector to the meson sector rather than being 
considered a restrictive prediction of the NJL model. We may conclude this section by stating 
that the NJL model in the collective approach provides the correct hyperon spectrum as far as 
the quantum numbers and the ordering of the hyperon states within a specified spin multiplet 
are concerned, however, the quantitative SU(3) breaking is underestimated due to the too 
small prediction for fx- 

It should be stressed that the collective approach together with the above described ap- 
proximation is not sensitive to the absolute value of the strange current quark mass but rather 
to the ratio m° s /m° as can be seen from eqs (D.10, D.ll and D.12). This procedure avoids 
uncertainties related to the choice of a special regularization prescription because the ratio is 
known to be insensitive on the regularization prescription while the absolute values for the 
current quark masses are not fc . 

There have been similar calculations by another group using the collective approach to the 
NJL model [125]. These authors choose not to perform the shift S — > S — rh° for the scalar 
fields (c/. section 3.2). Of course, when expanding the fermion determinant to all orders that 
treatment should be identical to the one presented here. However, this is not the case when 
approximations are involved. In ref. [125] furthermore flavor independent constituent quark 
masses (ra s = m) are assumed and the determinant is expanded in terms of m° s . It should be 
mentioned that an expansion of the action in m° s in the sense of chiral perturbation theory 
is claimed not to converge for the "physical" value of m° s [18]. In any event, it is obvious 
that this treatment is sensitive to the absolute value of m° s and thus also to the regularization 
prescription. In order to accommodate the empirical value m° s 150MeV the authors of 
ref. [125] employ a double step proper-time prescription. The approximation (ra s = m) also 
implies the identity fx = fix (4.44). It is therefore obvious that in that treatment the hyperon 
mass splittings cannot be described properly. In order to get reasonable agreement with the 
experimental spectrum m° s had to be increased by about 20%. The relation between m° s and 
the kaon mass rriK is almost linear. The gradient expansion (7.91) reveals that this increase 
of m° s just corresponds to substituting the physical kaon decay constant. Furthermore the 
double step regularization scheme of ref. [125] involves two more undetermined parameters. 
Not surprisingly this makes possible a further increase of the symmetry breaking parameter 
7. It should be remarked that within the proper time regularization a consistent computation 
of the parameters of the collective Hamiltonian yields a larger 7 when the shift S — > S — rh° 
is actually performed'. This indicates that the results on the hyperon mass splittings are 
somewhat sensitive to the regularization prescription rather than the parametrization of the 
rotating fields. 

7.5.2 Bound state approach 

The bound state approach represents an alternative treatment for the description of hy- 
perons in soliton models. It has originally been applied to the Skyrme model [29] and later 
on experienced sizable extensions [126, 127, 128]. In the case of the NJL model the bound 
state approach provides, besides the description of hyperons, an apparent application of the 

& For a compilation of relevant articles see ref. [18]. 

'Our numerical computations in the proper time scheme show that the 0{m\~) contribution to 7 is about 
10% larger when the shift is performed and different constituent quark masses are assumed. 
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formalism for treating mesonic fluctuations off the chiral soliton [97]. The general aspects of 
this formalism have been explained in section 7.3. 

As shown in ref. [99] and mentioned in section 7.3 the mesonic fluctuations decouple with 
respect to their grand spin and parity quantum numbers. Since kaons carry isospin 1/2 the 
grand spin of these fluctuations is half-integer. It is therefore sufficient to only consider kaon 
modes in eqs (7.40, 7.48) and disregard other pseudoscalar fluctuations i.e. 



^ ( K(r,uj) , 

V (r } u) = and £ rj a (r, u)X a = ^ A V n ' ; ) (7.92) 



wherein K(r,u) is a two-component isospinor. The bound state is expected to appear in 
exactly that channel which contains a zero mode in the symmetric (m° = m°) case. An 
infinitesimal vector transformation in strange direction may be associated with 

K (r) = r .rUo^^y (7.93) 

Here Uq denotes an arbitrary 2x2 space-time independent unitary matrix fixing the isospin 
orientation. This form of the "would-be" zero mode" 1 suggests the following ansatz for the 
kaon bound state wave-function 

k(r,u) = r -rO(cu,r) with fi(r, u) = ( ) (7.94) 

wherein 0(r, u) is a two-component isospinor (not to be confused with the angular velocity J?) 
which only depends on the radial coordinate r and the frequency u. The angular dependence is 
completely given by r ■ t and thus (7.94) represents a P-wave kaon. Due to isospin invariance 
the kernels (7.57) corresponding to are unit matrices in iso-space 

A {2) [Sl] = f +CO —I ( drr 2 [ dr'r' 2 & 2 \u; r, r')rf(r, co)fi(r', u) 



+ 



drr 2 $ {1) (r)tt^(r } uj)rL(r } uj)}. (7.95) 



Explicit expressions for the kernels $( 1,2 ) are given in appendix E. 

Varying (7.95) with respect to yields a homogeneous linear integral equation 

r 2 jy dr'r' 2 ^ {2) (uj;r } r')n(r' } uj) + $ (1) (r)0(r, cu)| = (7.96) 

which in fact is the Bethe-Salpeter equation for the kaon bound state in the soliton back- 
ground. It is the analog of the bound state equation in the Callan-Klebanov approach [29] 
to the Skyrme model. The fact that the kernels $( 1 )>( 2 ) are unit matrices in isospace leads to 
identical Bethe-Salpeter equations for the isospinor components a(r,u) and &(r, u) implying 
that both have the same radial dependence r) w (r). Thus Fourier amplitudes a 4 -(u;) may be 
defined as 

a(r,u) = i] UJ (r)a 1 (ijj) and b(r, u) = rj CJ (r)a 2 (u}). (7.97) 



m Actually, the bound state is nothing but the "would-be" Goldstone boson for flavor rotations into strange 
direction. 
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In the process of quantization the Fourier amplitudes a 4 -(u;) acquire the status of creation- and 
annihilation operators. 

In order to determine the normalization of r) w (r) (which cannot be deduced from the Bethe- 
Salpeter equation (7.96)) it is helpful to compute the strangeness charge S associated with 
the kaon bound state. At the microscopic level S is defined as the expectation value 



s 



j VqVq j d 3 r q^(-S)q exp(z J d 4 xC) (7.98) 



wherein C is the NJL - Lagrangian of eq (3.1) and S = diag(0,0, 1) is the projector onto 
strangeness. Note that in the standard convention a strange quark possesses strangeness — 1. 
The expression (7.98) is expanded up to second order in the fluctuations. This calculation has 
been performed in ref. [99] and may be summarized as 

S = j J dr j Jr'$ 5 (cu;r, r')?y* (r)^(r') ^aj(cu)ai(cu) + a\(u)a 2 {io)^ . (7.99) 

The symmetric bilocal kernel $s(iu; r, r') is given in appendix E. Since the strangeness charge 
has to be quantized to be integer-valued, eqn. (7.99) provides a suitable normalization con- 
dition for the radial function r) w (r) 

dr J dr'$s{w; r, r')?y* (r)r) w (r') =1. (7.100) 

Obviously, the sign of the expectation value S is determined by the dynamics. In nature 
only hyperons with negative strangeness are observed. Thus a major criterion which the 
NJL model has to satisfy is the existence of a bound state with negative strangeness while 
the Bethe-Salpeter equation (7.96) should not possess solutions with \u\ < rriK and positive 
strangeness. 

In figure 7.1 the radial dependence of the bound state wave-function normalized according 
to (7.100) is displayed for several values of the constituent quark mass m while the corre- 
sponding bound state energy lo is listed in table 7.8. 

Different occupations of the bound state remove the degeneracy of states with different 
strangeness. However, the system consisting of the static soliton and the kaon bound state 
has still to be projected onto states with good spin and isospin quantum numbers in order to 
also remove the degeneracy of baryons with identical strangeness, say the E and the A. This 
projection is achieved by the usual cranking procedure for the exact zero modes, which corre- 
spond to isospin rotations. The associated collective coordinates are introduced analogously 
to the procedure explained in section 7.1, however, also the fluctuations have to be taken into 
account. This is accomplished by writing 

M = R(t)ZoZf(E)ZfZ ltf(t) with R(t) G SU(2). (7.101) 

Obviously, a global isospin rotation U g corresponds to the substitution R(t) — > U g R(t). This 
indicates that the total isospin is carried by R(t). Since (E) commutes with the isospin 
generators the ansatz (7.101) is equivalent to 

Co -> R^ZoR^t) and K^R{t)K. (7.102) 

This demonstrates that K has lost its isospin. It should be remarked that although the ansatz 
(7.101) is quite different from the one assumed in the original Skyrme model calculation [29] the 
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Bound State W a v e Function 




Figure 7.1: The radial dependence of the bound state wave-function r) w (r) for various con- 
stituent quark masses m. r) w (r) is normalized according to (7.100). (Taken from ref. 



substitution (7.102) is identical. As for the two flavor version of the model angular velocities 
J? are introduced, see eq (7.4) 

% -t ■ Q = R^(t)R(t). (7.103) 

Although the collective rotations are not the only time-dependent held components when 
fluctuations are present, the identity 



\r„M 



DM 

D »m < 7 - 104 > 



3 



nevertheless holds. The rotation matrix Dij for two flavors is defined in section 7.1 before eq 
(7.19). Adopting the previous argumentation on Noether charges (c/. eqs (7.18, 7.19)) the 
total isospin is given by 



d L(R, O ) 



h = -D^ \^ ' (7.105) 



wherein L is the Lagrange function which is a functional of the chiral angle as well as the 
kaon fluctuations K . It is furthermore instructive to take into account the hedgehog structure 
of the soliton and define the momentum conjugated to the angular velocity J? 

Js = ^ (7.106) 

as the spin carried by the soliton. This leads to the identity I 2 = J 2 & . Still Jq needs to be 
related to the total spin. On the microscopic level the latter is defined as the expectation 
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value 

(J) = J DqDq J d 3 r q^ Jq exp (z^njl) (7.107) 

wherein J is the spin operator for a Dirac spinor and *4njl = / g? 4 ^£njl denotes the action 
associated with the NJL Lagrangian (3.1). Since the spin operator commutes with the iso- 
rotations R(t) the transformation into the rotating frame q = Rq' is straightforward 

(J) = J Dq'Dq' J d 3 r q'^Jq' exp {iA' mij ) • (7-108) 

Here *4n JL represents the NJL action in the rotating frame which also contains the Coriolis 
term 

*4njl = / d 4 x (c NJL - l -q'^r ■ Oq^ . (7.109) 

Substituting the definition of the grand spin G yields 

(J) = j Dq'Dq' j d 3 r 9 't (g - ^ q' exp (tA' N]L ) . (7.110) 

In this expression the soliton contribution to the spin Jq may be identified by differentiating 
^njl with respect to the angular velocity J? 

r 1 r) A 1 

(J) = (G) + J Dq'Dq' exp(zAU) = (G) + J e . (7.111) 

Stated otherwise: the spin carried by the kaons 

J K = (J) - J e = (G) (7.112) 

is identical to the grand spin expectation value. This is generally expected based on the 1/Nc 
expansion for baryons, see e.g. ref. [129]. Since the eigenstates of the static Hamiltonian are 
also eigenstates of the grand spin operator, i.e. G^l/i) = M^lfi) the result for the kaonic 

(2) 

spin is easily obtainable. Repeating the calculation leading to A F , however, including the 
grand spin projection quantum number when taking matrix elements provides after a 
straightforward calculation [99] the third component Jk3, which in turn is proportional to 

a\(Lo)ai(io) — a|(cu)a 2 (cu) = a } ( T 3)i] a r Exploiting rotational invariance one finds 

Jk = ~\j ^d(Lo) ^ a](u)r %3 a 3 (u) j . (7.113) 

The spectral function d(u) is given in appendix E. Here it is interesting to discuss a conceptual 
difference in comparison with Skyrme type models. In these models the bound state approach 
involves classical static fields as well as kaon fluctuations. The former have vanishing grand 
spin while the latter carry grand spin 1/2. This yields the identity d(u) = 1 [30]. In the 
NJL model the situation is different because the functional trace involves quark spinors with 
arbitrary grand spin. These spinors get polarized by the soliton field as well as the kaon bound 
state causing d(u) to deviate from unity. 
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In the context of the collective quantization the main task now consists of determining the 
coupling of the kaon fluctuations to the angular velocities Q . For this purpose the fermion 
determinant is first expanded in J? 



A T = A?{f2 = 0) + X; n + 0(&) = Af + A$ + 0(&). (7.114) 

where the contributions of order J? 2 yield the non-strange moment of inertia a 2 (7.15). Here 
the interest is on the linear term which, in contrast to the two flavor model, does not vanish 
but rather provides the coupling between the fluctuations and the collective rotations. In the 
second step of the calculation, A^p is therefore expanded up to quadratic order in the kaon 
fields. To this end the J?-dependent terms of the Lagrange function may be extracted 

L n = \ a2Q2 -\\ % c{ ^ Q ' (.£ 4m^M j + • • , (7.H5) 

Again we refer to ref. [99] for details on this calculation while the explicit expression for the 
spectral function c(u) is presented in appendix E. 

As already mentioned hyperons with different strangeness are constructed by various oc- 
cupations of the kaon bound state i.e. the solution to the Bethe-Salpeter equation (7.96) with 
< rriK- Thus 

c := c(lo ) and d := d(io ) (7.116) 

are projected out from the spectral integrals (7.115,7.113). With this restriction to the bound 
state the collective Hamiltonian obtained from Lq reads 

H n = 2^[ Je + t £ a! Ttj(lj ) = h ( Je + xJk)2 (7 - 117) 

where the parameter x = —c/d has been introduced. The arguments of the annihilation 
(creation) operators have been omitted. According to the above discussion J@ and J k may 
be related to spin J and isospin I quantum numbers of the considered baryon. This yields 




H Q = ^ U J ( J +!) + (! " XW + 1) + ( d f) 2 S(S - 2) ) (7.118) 



because it can be shown that (j2i J=1 OjTijdj^j = S(S — 2) [30]. This term is already of fourth 

order in the bound state wave-function and has to be dropped for consistency. Each substitu- 
tion of a non-strange valence quark by a strange one changes the valence quark contribution 
to the energy from e va i to e va i — lo . The mass formula for physical baryons thus becomes 

M B = E d + Su + ( X J(J + 1) + (1 - xW + 1)) (7-119) 
2cr 

wherein E c \ is the classical energy (6.46). This expression is similar to the Skyrme model 
result [30], however, here the ratio x = —c/d had to introduced with d ^ 1. The term linear 
in lo should also come out by performing a calculation similar to that leading to eq (7.60) 
when the terms of odd powers in u are properly accounted for. 
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Table 7.8: Parameters for describing the hyperon spectrum as functions of the constituent 
mass m. Also listed are the empirical values which are obtained by the consideration of 
certain mass differences. (Data taken from ref. [99].) 



m(MeV) 



350 



400 



450 



500 



empir. 



cuo(MeV) 



207.1 

-0.20 

0.90 

0.22 

8.30 



182.6 

-0.36 

0.89 

0.40 

5.80 



163.6 

-0.46 

0.89 

0.52 

4.78 



148.8 

-0.53 

0.89 

0.60 

4.17 



-189.5 



c 



or 




0.62 
5.12 



Before going into the details of the numerical results™ it is important to note that bound 
states have only been obtained in the region —rriK < cu < and that these carry negative 
strangeness according to eq (7.99). No bound state has been observed in the interval < u < 
m K . 

As in the collective treatment only the mass differences with respect to the nucleon mass 
are considered. In this case eq (7.119) contains three quantities (a 2 ,cu ,x) which determine 
seven mass differences. Thus the baryon mass formula (7.119) may be inverted relating a 2 } io 
and x 1° hyperon mass differences. These relations may be employed to obtain empirical data 
for a 2 } io and \- We list these together with the calculated values for various constituent 
quark masses m in table 7.8. The empirical data are reasonably well accommodated in the 
region 400MeV < m < 450MeV although \ i s somewhat too small. 

Again, the only free parameter m is chosen to get a best fit to the experimental mass 
differences. The A-nucleon mass difference is the same as in the two flavor model, 3 /2a 2 . 
In order to reproduce this mass splitting kaon fluctuations need not be considered. This 
simplification has been used in ref. [99] to fix m = 430MeV such that a 2 = 5.12/GeV. Then the 
six other mass differences are predicted. The results are listed in table 7.9. There also the case 
is presented when fx is fixed to its experimental value rather than /„.. This yields a different 
value for the cut-off A and thus m needs to be readjusted in order to get a 2 = 5.12/GeV. 
This does not significantly alter the predictions on the mass differences since the ratio fx/fir 
is almost identical in both cases. Comparing with table 7.7 it is found that the bound state 
approach improves on the mass differences, especially for the spin 1/2 hyperons. 

The bound state approach may be extended to the investigation of kaonic channels others 
than the P-wave. A prominent case is represented by the S-wave since it in principle allows 
one to explore the odd parity A hyperon which is observed at 1405MeV experimentally. In the 
Skyrme model this state has been studied intensively, see e.g. refs.[29, 127]. The corresponding 
NJL model calculations are reported in ref. [130]. The kaon S-wave is described by the ansatz 



Although Ks possesses the same grand spin (G = 1/2) as the P-wave (7.94) these two channels 
nevertheless decouple due to the opposite parity. As for the P-wave, a Bethe-Salpeter equation 
for the radial function rj^(r) in the Fourier decomposition 




(7.120) 




and 





n 



The numerical method is explained in ref.[121]. 
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Table 7.9: The mass differences ol the low-lying | and | baryons with respect to the nucleon 
in the bound state approach to the NJL model. The up-quark constituent mass m is chosen 
such that the A-nucleon mass difference is reproduced. All data (from ref. [99]) are in MeV. 





A fixed 


fi< fixed 


Expt. 


A 


132 


137 


177 


E 


234 


247 


254 




341 


357 


379 


A 


293 


293 


293 


E* 


374 


375 


446 




481 


485 


591 





613 


622 


733 



can be derived. The solution of this equation determines the bound state energy, u>s • The 
quantization of spin and isospin proceeds as for the P-wave yielding cs and o?s, which denote 
the coupling of the S-wave bound state to the collective rotations and the contribution of the 
S-wave bound state to the spin, respectively. Following the quantization of the P-wave bound 
state it is then straightforward to derive the mass formula for the odd parity A 

M(odd parity A) = E d -u s + ^j (7.122) 

8 or 

wherein the ratio \s = ~cs/ds does not dependent on the normalization of f]^(r). 

Since the A(1405) is only about 40MeV below the kaon-nucleon threshold the bound state 
eigenvalue for the S-wave is expected to lie at cos ~ — 450MeV. Fluctuations possessing ener- 
gies as large as this value may raise problems because the NJL model is not a confining theory. 
In case the valence quark is only slightly bound, i.e. e va i — m is small, kaon energies of several 
hundred MeV may scatter the occupied valence quark state into the strange continuum. This 
situation corresponds to having an undesired nucleon-strange quark threshold, En s , below the 
kaon-nucleon threshold. In that case the second term on the RHS of eq (7.54) develops a pole, 
the position of which determines E^ s . Since for kaon fluctuations the perturbation carries 
unit strangeness, in (7.54) is the eigenvalue of a free Dirac-Hamiltonian in the strange 
sector. Therefore m s is the smallest eigenvalue which determines \En s \ ~ m s — e va i. It 

turns out that \En s \ is a monotonously rising function of the non-strange constituent quark 
mass m and that for m > 450MeV \En s \ lies above the kaon-nucleon threshold [130]. The 
predictions obtained in this region on the odd parity A hyperon are displayed in table 7.10. 

The resulting mass difference between the odd parity A hyperon and the nucleon is insensi- 
tive to the constituent mass m. As shown above, the A-nucleon mass difference is reproduced 
for m 430MeV. Then, unfortunately, the threshold E^ s lies slightly below the physical 
kaon-nucleon threshold. Nevertheless a reasonable description for both the relative position 
of this hyperon as well as the A can be obtained for m 450MeV. 
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Table 7.10: Parameters for the baryon mass formula (7.122) and the prediction for the mass 
of the odd parity A hyperon relative to the nucleon Ma s — Afjy as functions of the constituent 
mass m. Also the resulting values for the A nucleon mass splitting are presented. (Data taken 
from ref. [130].) 



m (MeV) 


a 2 (1/GeV) 


Lo s (MeV) 


X 


A s - N (MeV) 


A - N (MeV) 


450 


4.78 


-461 


0.46 


418 


314 


500 


4.19 


-461 


0.54 


419 


358 


550 


3.74 


-459 


0.57 


416 


401 


600 


3.42 


-456 


0.60 


411 


439 



8 Summary 

We have presented the description of baryons as chiral solitons within the NJL model, 
a microscopic model for the quark flavor dynamics. Within QCD we have stressed the role 
of chiral symmetry, its spontaneous breaking and its consequences in hadron physics as a 
basis for modeling the strong interaction. Furthermore, we have reviewed some important 
features of QCD in the limit of infinitely many color degrees of freedom. This provides the 
basic motivation for the picture that baryons emerge as soliton solutions of an effective meson 
theory. In the framework of QCD a realization of the soliton picture of baryons is not feasible. 
Therefore we have approximated the quark flavor dynamics by the NJL model which, like 
QCD, is invariant under global chiral transformations and breaks this symmetry dynamically, 
which is reflected by a non-vanishing quark condensate in the ground state. 

Using functional integral techniques the NJL model can be converted as an effective meson 
theory. This effective meson theory is highly non-local. At low energies one can, fortunately, 
resort to an gradient expansion, which describes the low energy (light flavor) meson dynamics 
reasonably well. In leading order the gradient expansion of the regular parity part of the 
effective meson theory yields the gauged linear cr-model. The associated massive gauge bosons 
are interpreted as vector and axial-vector mesons. In the limit of infinitely large (axial) vector 
meson masses this approximation to the effective theory reduces to the Skyrme model. The 
gradient expansion of the irregular parity part of the effective meson theory, which involves 
only odd numbers of spatial components of Lorentz-vectors, provides the Wess-Zumino action 
in leading order. The latter result furthermore allows one to identify the baryon current with 
the topological current of the Skyrme model. It should be emphasized that this gradient 
expansion refers to the vacuum (or sea) part of the action only. Hence the Skyrmion picture 
implies that the baryon number is carried by the polarized Dirac sea of the quarks. Stated 
otherwise, a soliton of baryon number B requires Nc • B valence quarks to be bound in the 
Dirac sea by the solitonic meson fields. This is an inherent assumption of all purely mesonic 
soliton models of baryons (Witten's conjecture). Restraining from the gradient expansion the 
NJL model allows one to test this conjectures. For this investigation we have considered the 
self-consistent soliton solutions of various field configurations in the NJL model. Although in 
order to describe a unit baryon number configuration explicit valence quarks have to be added 
to the chiral soliton of the pseudoscalar fields only, the baryon number is indeed carried by the 
Dirac sea once also the axial-vector meson fields are included. This result strongly supports 
the Skyrmion picture of the baryon. It furthermore provides a microscopic explanation for 
the fact that explicit valence quarks and (axial) vector meson degrees of freedom represent 
alternative approaches to describe baryons as solitons. 



100 



In order to simplify the investigation of static baryon properties we have resorted to the 
model containing only the pseudoscalar field while explicit valence quarks are present. In 
order to generate states with good spin and isospin quantum numbers we have adopted the 
semi-classical quantization scheme. In this scheme the time dependent coordinates, which 
parametrize the (iso) rotational zero modes, are quantized canonically. Since this treatment 
is analogous to the quantization of a rigid top it yields (iso) rotational bands of baryons. 
The energy splitting of the low-lying band members is of order 1 /Nc while the soliton mass 
is of order Nc- Like in Skyrme type models the absolute values for the baryon masses are 
generally predicted too large if 0(N C ) quantum corrections corrections are neglected. We have 
estimated these corrections and found that they may considerably reduce the baryon masses 
towards their experimental values. 

When extending this quantization prescription to the case of three flavors and including 
SU(3) symmetry breaking explicitly a reasonable description of the hyperons is obtained. 
Furthermore, an excellent description of hyperons within this approach is somewhat restricted 
because the NJL model predicts to small a ratio between the kaon and pion decay constants; a 
problem inherited from the meson sector. A alternative description of hyperons within soliton 
models relies on the existence of a strongly bound kaon meson in the background of the chiral 
soliton. As a matter of fact this meson becomes a zero mode if the SU(3) symmetry breaking 
is ignored. For the NJL model this treatment apparently provides a better description of 
the spectrum of the low-lying hyperons, especially for the J v = | + states. When applying 
this approach to excited states one is confronted with the existence of unphysical quark- 
(anti)quark thresholds. These reflect the non-confining character of the NJL model. Of 
course, any attempt to extend or refine the NJL model in such a way that this problem is 
avoided would be interesting. 

Nevertheless the NJL soliton model appears to be a well-established approach to study 
the low-energy properties of the nucleon. Especially it provides an apparent opportunity to 
compare the contributions of valence and sea quarks to various baryon observables. 
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Appendix A 

In this appendix we demonstrate the calculation of the Jacobi determinant (4.59) using 
Fujikawa's functional integral method [49]. In order to make the computation of the change 
of the integral measure, reflected by the relation J(0) 7^ 1, more transparent we will use a 
few simplifications. First, we will only consider infinitesimal chiral transformations, i.e. small 
0. Since the axial anomaly occurs only for the abelian subgroup £7a(1) we assume that 
is proportional to the unit matrix in flavor space, oc Up. This assumption will make the 
computation more feasible as commutes with other flavor matrices. The whole calculation 
may be done retaining these commutators with the result that the traceless parts of do not 
contribute. Second, we will work with a Dirac operator containing besides the usual kinetic 
and mass term only a vector field 

il} = + f) - m° =: if- m°. (A.l) 

In Euclidean space we are working with hermitian Dirac matrices and antihermitian vector 
fields the operator if is hermitian with respect to the ordinary scalar product. Therefore the 
eigenfunctions of the operator if can be chosen to form a complete and orthonormal set 

J d 4 xcpl(x)cp m (x) = 6 nm 

l>„(a)4(y) = S(x-y). (A.2) 

n 

We expand the dynamical quark field, i.e. the functional integration variable, in terms of the 
functions (p n 

<l( X ) = J2 anL Pn{x) 

n 

q(x) = J2 a n^n{ X ) (A. 3) 

n 

The expansion coefficients a n and a* are independent anticommuting Grassmann variables. 
The integration measure can now be written as 

DqDq = ~^\_da m da* m (A. 4) 

m 

The chirally transformed quark field (4.58) is given by 

X(x) = ^d m Lp m (x) = ^a m e l@i * x)lb Lp m (x) 

m m 

X(x) = £^R = EM^Re 80 ^. (A.5) 

m m 

where the expansion coefficients for the transformed quark field are related to the original ones 
by 

n 
n 

C mn = (m\e t@ ^\n) = [ d 4 x^Jx)e t@ ^'^ n (x) 



S mn + i I d 4 xQ{x) V l{x) 75Vn {x) + 0(0 2 ). (A.6) 
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Here we have used the orthonormality of the eigenfunctions <p n . As the sets a n , a* and a m , a* m 
are Grassmann variables we obtain 

J(0) = (DetC)" 2 . (A.7) 
For infinitesimally small the matrix C is given by 

C mn = S mn + (m\i l5 Q(x)\n) + 0(0 2 ). (A.8) 

Defining the matrix 

e mn = (m\i~f 5 Q(x)\n) (A. 9) 

and shorthandly writing C = 1 + e the functional determinant of this matrix is evaluated to 
be in first order in 



DetC = Det(l + e) = exp (Tr(e) + 0(e 2 ) 

= exp f d 4 xe{x)Y,ipl-f 5 i Pn {x) ] j + ... . (A.10) 

As it stands the operator appearing in the exponent is not well-defined. In order to attach a 
suitable definition we have to regularize (A. 10) and to remove the regularization at the end of 
the calculation 



B(x) := lim J2 ( fh^~ {Xn/A)2 ^n(x 



lim tr7 5 (x|e \x). 



A^oo 



(A.n; 



To apply the heat kernel expansion we introduce the "proper time" r = 1/A 2 as an expansion 
parameter and use A = —d^d^ + m 2 as unperturbed operator. Expanding now the matrix 
element 

J. oo 

(x\e- T V'\y) = (x\e~ TAo \y) £ h k (x, y)r k (A.12) 

k=0 

we obtain 

1 oo 

B ( x ) = 77^ lim . E r k - 2 tv{ l5 h k {x, x)). (A.13) 



(4tt 



k=0 



Note that only the heat coefficients h 0} hi and h 2 in the sum (A.13) contribute to B(x) because 
the coefficients of the higher order terms are suppressed by a factor r k ~ 2 and therefore vanish 
in the limit r — > (A — > oo). As 

ti^ 5 h (x } x) = tr7 5 = 

tr^h^x.x) = -^(757^7^)^ = ( A -14) 

only one term contributes 

tr l5 h 2 (x, x) = N c tT F F^F^. (A.15) 
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Here is the field strength tensor corresponding to the vector field and is the dual 
tensor. Note that a flavor trace has to be taken in (A. 15). Eqs. (A. 13) and (A. 15) can now 
be combined to the result 



N 

B(x) = -f-MFF^F, v (A.16) 



finally yielding the expression (4.61) for the Jacobian J(0): 

J(0) = (DetC) -2 = exp / d 4 xQ(x)tr F (F F)^ . (A. 17) 

Using the invariance of the classical action under (4.58) one can easily derive the anomalous 
Ward identity. Under infinitesimal chiral rotations the quark bilinear in the action transforms 
as 

q(i4-m°)q = xe"' eW75 (^-m O )e" ,0(a;)75 x 

= x(¥-^°)x + ^0(x)x7 A1 75X + 2zm o 0(x)x759 + C(0 2 ) (A.18) 

which leads to the following identities for the effective quark action 

Af = Trlog(zj$— m°) = log J DqVq exp j d 4 xq(ii — m°)q^ (A. 19) 

= iogJ(e) 

+ log j T>xT>xexp (i j d 4 xx{ii- m°)x + dv&WxYlsX + 2zm°0(:c)x75x) • 

As the first equation is independent of the last one has to be, too. Therefore, 
SA F 







SQ 

Defining the axial singlet current and density 



n n = -i^{tT F (FF)) - d.ixYlsx) + 2zm°(x75x)- (A.20) 



Js = (X7 A1 7sX) 

h = (X75X) (A.21) 



we finally arrive at the anomalous Ward identity (4.69) 



= 2im°j 5 - i^Ltv F ((FF)). (A.22) 
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Appendix B 

In this brief appendix we present the explicit representation of the eigenstates \v) of the 
Dirac Hamiltonian h (6.42,6.63) in coordinate space. As indicated in chapter 6 we first con- 
struct eigenstates to the grand spin operator 

G = Z + f + ^' (R1) 

which is the vector sum of orbital angular momentum Z, spin <r/2 and isospin r/2. G acts on 
the quark spinors. These are obtained by first coupling states with orbital angular momentum 
/ with 5 = 1/2 states to form states carrying total angular momentum j and projection j 3 

\m=Y,C>» h Jrn)\\a z ). (B.2) 

m,<r 3 z ^ 

Subsequently these states are coupled with isospin t = 1/2 states yielding the grand spin 
eigenstates 

\IjGM) = Y,C^Jjj 3 )\ 1 -t 3 ). (B.3) 

33 ,T3 

This coupling scheme obviously leads to the selection rules (6.45). Since the upper and lower 
components of Dirac spinors transform oppositely under the parity transformation the general 
form of Dirac spinors, which are eigenstates of G and parity II = ±1, are given by 

^ " \ff^ 1 \r)\G + lG+\GM)) + K-ff^)(r)\G-lG-\GM)j [aA) 

_ f^- ;1) (r)\G + 1G+ \GM) \ ( ig<f>-*\r)\G - 1G - \GM) \ . > 
^ -\ -f( G >-V(r)\GG+±GM) J + K fl G ~' 2 \r)\GG - \GM) J' ^- b) 

The second superscript labels the intrinsic parity Ili ntr and enters the parity eigenvalue via 
II = ( — 1) G X Ilintr. Ilintr represents a good quantum number since both grand spin and par- 
ity operators commute with h. Furthermore, the grand spin invariance is very helpful when 
choosing ansdtze for meson fields other than the chiral field, see section 6.2. The diagonaliza- 
tion (6.43,6.59) yields the radial functions g\^ ,+,1 \r), etc. as linear combinations of spherical 
Bessel functions (the solutions to the free problem). E.g. 



/ (G,+;i) (r) = '£V vk [Q]N kS gn(E kG )y/l - m/E kG j G+1 (k kG r) (B.6) 

k 

where the eigenvectors" K,fc[©] are obtained by diagonalizing h. In case that the Hamiltonian is 
not Hermitian (6.63) these eigenvectors are complex. N k are normalization constants. E kG = 



±ym 2 + k\ G denote the energy eigenvalues in the absence of the soliton. The momentum 
eigenvalues k kG are subject to the boundary condition j G (k kG D) = [81]. There are also 
other boundary conditions discussed in the literature [62, 98] and the pertinent boundary 
condition depends on the problem under consideration. 



Grand spin and parity indices are omitted. 
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Appendix C 

In this appendix the equations of motion for the meson fields as taken from ref. [73, 75] 
are listed. These equations result from the stationary condition SE/S(p for the Minkowski 
energy functional (6.39). The functional derivatives of the one-particle energy eigenvalues e v 
are extracted from the Euclidean Dirac Hamiltonian (6.63). 

It is appropriate to define densities" 

p{x,y) 



N c {p™\x, y) + p™(x, y)\ and b(x, y) = N c {b™\x, y) + b™ c (x, y 

Vv a \{^ va \{x)tp va l(y) + ^V'val(*)V'val(?/)} 

J2 {f R (e v /mM*)My) + fi(e,/A)%M*)My 



6 val ( 



x,y) 



??vai \ ^ va i(x)tpljy) - ^ij val (x)tpljy] 



*,y) = E\fi(^/ A )^M^ky)-fR(^J^M^ky 



(C.l) 



Real (3?) and imaginary (^s) parts refer to the expansion coefficients K,/^] of the free basis 
(cf. eq (B.6)). The regulator functions 



|sgn(e^)erfc 



rsgn e 



fi )erfc 



-sen e, 



erfc 



+ ^(4/A)exp( - (e?/A) 5 
Ai not regularized 

Ai regularized 



Ai not regularized, 

(C.2) 

Ai regularized 

(C.3) 



reflect the derivative of the energy functional with respect to the single quark energy eigen- 
values. The specific form of the equations of motion SE/S(p finally becomes 



$(r 
sin 0(r 
uj(r 
G(r 
F(r 
H(r 



m° 
m 
m 



9v 



cos0(r 
tr 



m 



ml fl 



-tr 



/JO 

/ — (cos0(r) + ijsT ■ r sin0(r)) p(x, x) 

J AtT 



Amy 



tr 



9y 
Am v 

9y 
Am\ 

9v 



jo 

Air 

dn 

Air 

dn 



sin0(r) — Z75T • r cos0(r)) p(x, x) 
b(x, x) 



Am v 



tr 



tr / ((7 x r) ■ r)p(x,x) 
r JO 

tr J -^H^-r)(T-r)-(a-T))p(x,x) 
r JO 

hrPtt (T ' f )( r • f ) - (°" • t))p(x,x). 

J Air 



(C.4) 
(C.5) 
(C.6) 
(C.7) 
(C.8) 
(C.9) 



Here the trace is taken with respect to Dirac and flavor indices only. Note that tr / JO b(r, r) 
represents the baryon charge density. In eq (C.6) isospin invariance has been assumed, i.e. 



m u 



m. 



m v . 



ipv{x) represents the coordinate space representation of \v) . 
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Appendix D 

In this appendix we present the explicit expressions lor the prelactors ol the collective 
operators which are relevant lor the description ol the hyperons as described in section 7.5.1. 

The explicit expression lor A ab may readily be found in rel. [62] where the corresponding 
calculation is described 

A^=^AM^f r (e„e v] A)(v\\ a \v)(v\Tp\ b T^ri (D.l) 
with the cut-off lunction 

sgn(e At )erlc ( f f J - sgn(e i/ )erlc 



friezes, A) 



D.2 



The contribution ol the Dirac sea to the coefficient 7 reads 

(HT/JAsTV). (D.3) 



2N C AM x . , 
7 = ^= — 2^ sgn(e At )erlc ( 



vac 

= 7= 7 y tin 1 e„ lenc 



A 



The second order terms in the expansion ol the vacuum contribution to the fermion de- 
terminant may be expressed as double mode sums over the eigenstates ol the static Dirac 
Hamiltonian (6.42) 

Nr , 

©r = fE M^, e„; A)(ii\\ a \v)(v\\ b \ii) (D.4) 

where the cut-off lunction /©(e^, t v \ A) is defined in (7.14). The symmetry breaking terms are 
found to be 

IT = %m - m s f J2 fr(e„ e„; N)^\T^\ a T^ \v) (v\T/3\ b T^ \fi) (D.5) 

with the cut-off lunction /r^, A) being defined in (D.2). 

The analogous expressions associated with the explicit occupation ol the valence quark 
level are found to be 



7 val = 


2 


m — 


P)val _ 


N c 


E 

/i^val 


aY = 

ab 


N C( 


— m< 


-pva\ 

1 ab — 


2 

-3AM 


m — 



(val|A a |//)(//|A 6 |val) 



(val|A a | / u)( / u|T/3A 6 Tt |val) 



(val|T/3A a Tt (fi\Tp\ b T^ |val) 



(D.6) 
(D.7) 

(D.8) 

(D.9) 



wherein |val) denotes the valence quark state \& va i- Finally the coefficients appearing in the 
mesonic part ol the collective action (7.85) are given by 



7 = y m ^ 



m° \ / m*\ ( m°\ /m s 

2 + — + 2 + HU — "1 



my \ m° 1 \ m 



j drr 2 (l - cos0),(D.lO) 
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17 = -Tf = —mi I — I - 1 ) ( — - 1 ) / Jrr^ (1 - cosO) , fD.ll] 

r * = f^g- 1 )^- 1 )/- 2 ! 1 — (!))• 

These various contributions finally add up to the collective Lagrangian which may be cast 
into the form 

1 3 f 7 f 

~^tJ2 D ^ D ^ ~ £ ^"AS" " -r 8 (f - AssAss) 

z=l a=4 
3 7 

+A T ^0 8J D 88 + A 5 ]T «A„. (D.13) 

i=l a=4 

The coefficients are sums of the quantities listed above 

a 2 = 0£ + 0£r, P = ©44 1 + ©™ + Pi, B = Vval + B™ c 

7 = 7 val + 7 vac + 2T™ 1 + 2TH c + 7 m , r T = r™ c + r™ c + 17, 

p pvac I pvac . pm p pvac pvac . pm 

i- S — 1 44 + 1 44 + 1 £, J-8 — — J-88 _ 1 88 + 1 8 5 

A T = A^ + A^, A 5 = A^ ac + A^ c . (D.14) 

Whenever necessary, the contribution of the trivial (£ = 1) field configuration has to be 
subtracted since the fermion determinant is normalized accordingly. 

In order to obtain the collective Hamiltonian from (D.13) it is important to note that the 
time-derivative of the collective rotations = fbcd^cDu only affects the right indices of the 
adjoint representation (see also eq (7.17)). Hence the momenta associated with dL/dQ a are 
the right generators R a of SU(3). These obey the algebra [R a} Rb] = —f a bcR c - Noting that 
[R a} R] = —RX a which also implies that [R a} Dib] = —if a bcDi c one easily verifies that 

[^,C(M)] = -%^. (D.15) 

Thus the quantization prescription reads 

r -(a 2 n a + A T D 8 a) = - J a , a=l,2,3 
R a = = ~(^a + A s D 8a ), a=4,..,7 (D.16) 

dUa { f B, a=8 

wherein J; (i = 1,2,3) denote the spin operators. This relation is the flavor SU(3) extension 
of (7.18). 

The Hamiltonian operator: 

H =-J2Ratta-L (D.17) 

(1=1 

may be diagonalized exactly by generalizing [131] the Yabu-Ando [28] approach to more com- 
plicated symmetry breaking terms. In the original Yabu-Ando approach only the 7(1 — Dgg) 
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type symmetry breaking term was considered. It should be noted that the constraint R 8 — 
y/3B/2 = commutes with the Hamiltonian and therefore is first class. Hence the baryon 
wave-functions live on the hypersphere R 8 = y/3/2 for B = 1. On this hypersphere only 
states with half-integer spin exist. Thus the Hamiltonian has the appreciated feature that 
its eigenstates are fermions [52, 26, 27]. The corresponding energy eigenvalue for baryon B 
is found to be given by the formula (7.90). The quantity csb appearing in that equation 
represents the eigenvalue of the SU(3) operator 

C 2 + /3 2 7 (1 - D 88 ) + /3 2 (A T /a 2 ) £? =1 D 8i (2Ri + A T D 8l ) + A s EL 4 D 8a (2R a + A s D 8a ) 

+/3 2 r 8 (l - D 88 D 88 ) + (3 2 T T £-=i D 8t D 8t + (3 2 T S EL 4 D 8a D 8a , (D.18) 

where C\ = J2a=i Ra denotes the quadratic Casimir operator of SU(3). Left SU(3) generators 
are constructed via L a = J2b=i D a bRb with I a = L a , (a = 1,2,3) and Y = 2Lg/y/3 being the 
isospin and hypercharge operators, respectively. A pertinent parametrization of the collective 
rotations in terms of eight "Euler angles" is given by [28] 

R = i?/(a,/3, 7 )e- 8i/A4 i? J (a / ,/3 / ,7 / )e-^ / ^ A8 . (D.19) 

Here Rj represents pure isospin transformations and, due to the hedgehog structure of the 
soliton, Rj corresponds to spatial rotations. For the parametrization (D.19) the generators 
R a are expressed as differential operators in the "Euler angles" a,f3,...,p. Employing the 
explicitly forms of the R a , as given in appendix A of ref.[131], the operator (D.18) can be 
formulated as a second order differential operator for the "Euler angle" v. This suggests the 
ansatz for the eigenfunctions [28] 

*Yii a ,jj a (R) = ("1) J " J3 E D[ a > L (a,M)fM L MMe^D^_ Ja (a' J', 7 '), (D.20) 

M L M R 

where the D's denote the SU(2) Wigner functions. The sum over the intrinsic projection 
numbers Ml^r is subject to the condition Ml — Mr = Y/2 — y/3 for B = 1. In order to 
finally obtain the eigenvalue csb (and thus the baryon spectrum) a system of coupled second- 
order differential equations for the isoscalar functions fM L M R (v) has to be integrated. The 
fM L M R (v) decouple for the various isospin channels and the integration is accomplished by 
standard (numerical) techniques. 
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Appendix E 

Here we briefly display the kernels $( 1,2 ) entering the Bethe-Salpeter equation (7.96) which 
represents the starting point lor the bound state description ol the hyperons. Also presented 
are the kernel $5 lor the strangeness charge (7.99) as well as quantities relevant lor the semi- 
classical quantization ol the kaon bound state. These expressions are computed in rel. [97] and 
the results are readily taken over. 

The local kernel $^)(r) acquires contributions Irom the meson part ol the action as well 
as those terms involving h( 2 ) 

*-.„, = -1^(1 + =,) (c„ s e + ^) 

-^??vai(m + m s ) j ^V'Jai( 7 ') M o(^)/3Mo(r)V'vai(^) (E.l) 
4 J1/A2 v Airs J 4tt <- ^ ns 



+ £e,e-'^J(r)/%(r)}. 



p=s 



The integral J(JO/47r) indicates that the average with regard to the internal degrees ol freedom 
has been taken. For convenience the unitary, sell-adjoint transformation matrix u 0} which 
represents a modified form ol the chiral rotation (7.42) has been introduced 

u (r) = if3~f 5 r ■ tT = f3 (sin® - i-( 5 r ■ rcos ^) • ( E - 2 ) 

Obviously u (r) acts as a unit matrix on the strange spinors. The bilocal kernel $( 2 )(cu; r, r') 
originates Irom the terms quadratic in and is symmetric in r and r' 

(h (2) { a N C( , 2 [dSl [dSl'f x , V)J al (r)u (r)V) p (r)V)|(r>o(r / )V) val (r / ) 

& ) {u;r,r) = — (m + m s ) /— /— <i]^i}^ 

4 J 4tt J 4tt [ ^ e va i - u - e p 

~ £ V'J(^)«o(r)V'p(r)V'J(r')«o(r')V',(^)^,p(^)} (E.3) 

li=ns ) 
p=s 

The regulator lunction appearing eq (E.3) has been obtained as 



2 2 

/CO I g ^Q~ St jJ, _|_ Q~ st p 
ds \j-{ ~ ^ t^ 2 ~~ + £ p ) 2 ]Ro(s;lo, e^, t p ) 

-2cue p i?i(s; u, e^, e p ) + 2ue tl R 1 {s; u, e p , e^)}. (E.4) 

The Feynman parameter integrals are defined in eq (7.52) 

<„, t.) = / *>dx exp (- S [(l - *)4 + «5 .M) . (E.5) 

Similarly the bilocal kernel $s(iu; r, r') for the strangeness charge (7.99) is given by 

~ ( >\ f ^ V- V'J a l(^)^o(^)V'p(^)V'p(^)^o(^)V'val(r / ) 

$5 {u;r,r ) = -rj val / — / — ^ 7 ^ E.6 

J 4tt J 4tt ^ (e va i - - e„) 2 
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+ / ^ / ^ /~ 2 j= £ ^{(1 - 2 a£ ;)^(r)t^(r) - ^„(r)t^(r) ^fr - r'^ 

dtt rdW ^ 



p=s 



X {-Ro(s; u;, e^, e p ) — + + e p ) ((e p — u)Ri(s; u, e^, e p 

To conclude this appendix we also list the spectral functions c(u) and d(u) which enter the 
relation between the spin and kaon spin expectation values (c/. eqs (7.113) and (7.115)). It 
is appropriate to list the contributions stemming from the explicit occupation of the valence 
quark level and the polarized Dirac sea separately. In what follows are understood as 
the perturbative parts of the Dirac Hamiltonian (7.43) with the properly normalized solution 
Vu( r ) t° the Bethe-Salpeter equation (7.96,7.100) substituted 

c(u) = c va i(u) + c vac (u) 

„ r (val|/t (2 )(cu,-cu)|/i)(/i|T 3 |val) 

CvalM = rjval 2^ [ Z + h - C - 

_ r (vaI|fe ( i)(a;)|/?)(/?|fe(i)(-a;)|//)(//|T3|vaI) 

T^val ( \( \ ' C ' 

M=«s l e val — e^Jl £ val — ^ — e pJ 



/•oo ^ ( „ g g se fl g g se J 

Cvac(^) = - J^^-J=\ (^l T 3k)(^l^(2)(^,-^)|^)^ ^ (E.7) 



fj,,i/=n.s. h> 

e 2 - e 2 



H,v=n.s. t v ) 

p=s. 

The relevant regulator function is given by a Feynman parameter integral 

k ,„m = /; i,^ + e , + e ,, exp (-.[(i - + «; - x ( i - ,V]) 

+ (2(1 - x )lo 3 + [e,, + t v ]u 2 - [2(1 - a;)(e J/ + t p ){t^ + e p ) - (^ - £,)> (E.8) 

2 2 \ 

+ (e A1 + e i/ )(e M + e p )(e i/ + e p )J 2 _ 2 exp (-s[(l - x)e^ - x(l - x)cu 2 ]^ 

Finally d(u) is found to be 

</.„m = -2 l , v „E^ <va " / ' i " (r r"' >< "' / '" ,( ,r"" Ya ' > 

t^s ( £ val - W - £p) 2 

/•oo A a •> r 

<W") = 2 -= }Z M,e-^{(l-2sel) (v\h {2) (r,u, 

-stf, (ii\h {1) (r,u)h {1) (r, -u) 

/•oo / 3 

"2 y i/A2 — Z M,(^(i)(r,cu)|z/)H/*(i)(r, + e„ - u) (E.9) 

x|i? (s;cu, e^, e v ) - s(uj + + e^) ((e^, - cu)i?i(s;cu, e^, e^) + ujR 2 (s;uj, e^))}. 
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